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FOREWORD 


This  report  contains  the  basic  unedited  Vu-graphs  of  the  presentations  at 
the  ’’Mechanics  of  Composites  Review”  sponsored  jointly  by  the  Non- 
metallic  Materials  Division  of  the  Air  Force  Materials  Laboratory  and 
the  Directorate  of  Aerospace  Sciences  of  the  Air  Force  Office  of  Scientific 
Research*  The  presentations  cover  current  in-house  and  contract  pro¬ 
grams  under  the  sponsorship  of  these  two  organizations. 

Since  this  is  a  review  of  on-going  programs,  much  of  the  information  in 
this  report  has  not  been  published  as  yet  and  is  subject  to  changes.  But 
timely  dissemination  of  the  rapidly  expanding  technology  of  advanced 
composites  is  deemed  highly  desirable.  Works  in  the  area  of  mechanics 
of  composites  have  long  been  typified  by  disciplined  approaches.  It  is 
hoped  that  such  a  high  standard  of  rigor  is  reflected  in  the  majority,  if 
not  all,  of  the  presentations  in  this  report. 

Feedback  and  open  critique  of  the  presentations  are  welcome.  Special 
thanks  are  due  to  James  M.  Whitney  of  the  Nonmetallic  Materials  Division 
for  his  effort  in  organizing  this  review.  Again,  suggestions  and  recom¬ 
mendations  from  all  participants  will  be  most  important  in  the  planning  of 
future  reviews*.^ — . 


Z.  M.  KELBLE,  Chief 
Nonmetallic  Materials  Division 
Air  Force  Materials  Laboratory 
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Moisture  Distribution 
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n  SINGLE  LAYER  -  Steady  State  Solution 
rn  EXPERIMENTAL  PROCEDURES 
ET  MULTILAYER  -  Unsteady. (numerical)  Solution 


SINGLE  LAYER  -  STEADY 


1.  T *  CONST 

a  a 

c  =  C  =  CONST 
a  a 

(STEADY) 

2.  INITIALLY  (t<0' 

C,  T  UNIFORMLY 
DISTR,BUTLD 


3.  K  »K(T)  ONLY 
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5.  D  =  D(T)  ONLY 
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2)  Edge  Corrections 
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GENERAL  PROBLEM 


(t) 
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RESULTS:  1)  Total  Moisture  Content  (WtGain)  as  a  Function  of  Time 

2)  Moisture  Content  (Wt  Gain )  of  each  Layer  as  a  Function  of  Time 

3)  Moisture  Distribution  in  Each  Layer 

4)  Temperature  Distribution  in  Each  Layer 
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WEIGHT  GAIN,  % 


SUMMARY  -  MOISTURE  PROBLfM 


1)  Exact  Steady  State  Solution 

2)  Experimental  Procedures 

3)  Test  Facility 

4)  Fully  Operational  Computer  Code  for 
General-Multilayered-Unsteady  Problems 
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2)  Temperature 
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MATERIAL  DEPENDENCE 
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2)  Fiber  Orientation 


NEEDED  :  COMPREHENSIVE  SET  OF  DATA 


1)  For  "Absolute" Values 
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4)  For  Guidance  in  Future  Tests 
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HT-S/(8183/137-NDA-BF?:  MEA) 

Hertz 

- 

- 

X 

JLJ 

HT-S/HYSOL  ADX-516 

Browning 

- 

X 

X 

- 

Hercules  HT-S/710 

Kerr  etal 

-  X 

-  X 

-  X 

HT-S/PDN 

Browning 

- 

-  X 

- 

Boron /AVCO  5505 

Hofer  etal 

X 

X 

X. 

X 

xj 

Boron /Narmco  5505 

Kaminski 

X 

X 

Browning 

- 

X 

X 

- 

GENERAL  CONCLUSIONS  -  ULTL? ATE  TENSILE  STRENGTH 


A.  Temperature 

1)  (P  and  7r/4 - Negligible  Effects 

2)  90^  - ►  Significant  Effects 

a)  Moisture  Dependent 

b)  60-90%  Max.  Loss  in  Strength 


B.  Moisture 


1)  0°  and  7r/4 


2)  90° 


Small  Effects 

a)  Negligible  Below  1% 

b)  Small  (20-30%  Max.  Strength  Loss)  Above  1% 

-►  Significant  Effects 

a)  Temperature  Dependent 

b)  60-90%  Max.  Loss  in  Strength 
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SUMMARY  -  ULTIMATE  TENSILE  STRENGTH 

1)  Extensive  Set  of  Data  for  Thornel  300  /  Fiberite  1034 

2)  Summary  of  Existing  Data 

3)  Interpretation  of  Existing  Data 

4)  General  Trends  Due  to  Moisture  and  Temperature 

5)  Guidelines  for  Future  Tests 


10 


ENVIRONMENTAL  EXPOSURE  EFFECTS  ON 
COMPOSITE  MATERIALS 

AFOSR  72-2214F 

Principal  Investigator:  Lawrence  J.  Broutman 
Illinois  Institute  of  Technology 


OBJECTIVES: 

1.  To  Study  the  Impact  Behavior  of  Advanced  Fiber  Composites. 

2.  To  Study  the  Influence  of  Fiber-Resin  Interface  Strength  on  Impact  Energy  Absorption. 

3.  To  Study  the  Effects  of  Moisture  Content  on  Performance  of  Carbon  Fiber  Composites. 
Determine  the  Relative  Roles  of  the  Interface  and  Matrix. 

3.1  Characterize  Moisture  Diffusion  as  a  Function  of  Applied  Stress  and  Temperature. 

3.2  Characterize  Interfacial  Moisture  Diffusion. 

3.3  Determine  Effect  of  Moisture  on  Strength  and  Impact  Behavior  by  Measuring 
Current  and  Residual  Properties. 

3.4  Analyze  the  Long  Term  Degradation  (11  Years)  of  Carbon  Fiber  Reinforced 
Epoxy  Resins. 
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Fobric  Style  I8[  FIN  112 
Z  -  6  0  3  0  Trea^ed 


Fabric  Style  181  FIN  112 
QZ  -  8  -  5069  Treoted 
(  Coupling  Agent  ) 


(  Coupling  Agent  ) 

Figure  30.  Effects  of  Envlromentol  Degradation  on  Impact  Strength  Of  Strongly-Bonded  Fiberglass-Polyester  Laminates 


20.68  □  dry 


No  S.T. 

No  S.T. 

XF-l-3563 

2-6030 

QZ-8-5069 

STYLE  7581 

STYLE  181 

STYLE  181 

STYLE  I8i 

STYLE  181 

No  FIN 

F1NII2 

FIN  112 

FIN  112 

FIN  112 

(Mold  Release)  (Coupling  Agent)  (Coupling  Agent) 

Figure  32.  Effects  of  Envlromentol  Degradation  on  the  Ductility  Index  of  Fiberglass-Polyester  Lamlnotes  In  Impact 
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INTERLAMINAR  FRACTURE 
SURFACE  WORK  (y),ft-lb/in2 


1  I 


o  WATER-SOAKED 
7  DAYS 


APPARENT  SHEAR  STRENGTH  (S„),  ksi 
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Interlaminar  Shear  Surfaces  For  (A)  Carbon  Fiber  and  (B)  Graphite 
Fiber  Specimens  Stored  in  Sea  Water  for  1 1  Years  (2000X) 


Table  7. 


Diffusion  Coefficients*  of  Ccrbon-Epoxy  Reexposed  to  Humid  Environments  After  11  Years 
Environmental  History  ^0  cm  ^/sec] 


Material 

Designation 

11  Yeors  History 

Re-exposure  Conditions 

Distilled  V/ater,  60°C 

98%R.H. 

,  60°C 

'>M 

^M'^C  ,nn 

- sr - X  too 

CD 

Dry 

89.04 

90.06 

-1.13% 

88.24 

87,^-') 

+.89% 

CWd 

Distilled  water 

114.81 

84.55 

+35.79% 

125.99 

82.11 

+53.44% 

CWj 

Sea  water 

122.36 

84.55 

+44.72% 

133.17 

82.11 

+62.18% 

**  measured  coefficient 
-  calculated  coefficient,  equation  (5) 
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Table  B.  Dlfljstcn  ^aotiiclcats  of  Carlnjn-L  paii)  Lorninofes 


Fiber/Rrsin 

Orlontotion 

Environmental 
Conditions  D[10 

-10  7,  1 

cm  /’.ccJ 

References 

T300/5208 

o^/t  “^5 

R.  H.  dO'^C 

32 

(2) 

T300/5208 

0/145 

SO^’oR.H.,  60°C 

9 

(3) 

T300/1031 

Q^/t  45 

80%R.H.,  60®C 

83,4 

(3) 

HMS/t034 

0/+  45 

80%R.H.,  60®C 

93.5 

(3) 

T300/1034 

+  45;  0 

water,  60®C 

11.3 

(1) 

Table  9.  Inferlominar  Shear  Strength  of  Corbon^Epoxy  Composites  Exposed 
for  n  Years  to  Environmental  Conditions  (x  10^  psl) 


Reference^ 

Test  Mode 

Environmental  History 

Value 

Dry 

Distilled  V/ater 

Sea  Water 

Current^ 

9.48+  .10 

6.21 +  . 20 

7.04+  .05 

10.2+  .1 

Partial 

drytng3 

9.58+  .20 

6.31  +  .20 

7.63+  .10 

Residue!^ 

9.81  +  .10 

7.24+  .10 

7.83+  .10 

Tested  by  Union  Carbide  in  1965 

tested  in  17/6,  tn.-r  u.uovin^  ii.e  specimens  from  the 

environmental  conditions 

tested  after  a  drying  period  of  10  days,  at  60°C  under  vacuum.  The 
water  content  in  each  case  is  shown  in  Figures  1-5. 

tested  after  complete  drying. 


Tolile  11.  Interlominar  Shear  Strength  of  Graphite -Epoxy  Specimens 
Exposed  for  11  Yecm  to  Environmental  Conditions  Dc  103  psi] 


Reference* 

Volue 

Test  Mode* 

Environmental  History 

Dry 

Disti lied  V/oter 

Sea  Woter 

Current 

2.59+  .20 

3.68+  .05 

3.60  +  .10 

3.3+  .1 

Residual 

2.61  +  .05 

3.34+  .05 

3.27+  .05 

*  See  notes  in  Toble  9. 
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Specimen 


Sea  Wotcf 
Hisfory 


Specimen 

"Dry"  History 

Distilled  Water 
History 

Sea  Wo1 
History 

Carbon 

Current  1 

— 

81.40 

73.67 

2 

— 

75.30 

76.35 

Average 

— 

78.35 

75.01 

Residual  ^ 

83.74 

84.90 

75.65 

2 

80.60 

78.56 

77.23 

Average 

82.17 

81.73 

76.44 

grqpbjil 


Current  1 

— 

66.98 

62.49 

2 

— 

68.83 

64.51 

Average 

— 

67.50 

63.50 

Residual  1 

60.24 

66.98 

66.76 

2 

64.79 

67.61 

64.56 

Average 

62.51 

67.30 

65.65 

Figure  3.  Moisture  Desorptlo-^  of  Grophite  Epoxy  Specimens  Exposed  to  11  Years 
Distilled  Water,  ord  Moisture  Acsorption  Upon  Reexposure  to  Humid 
Environments 


Figure  d.  Fi'ocrionol  V^ciu’f  Go'^’  of  Caroo'T£"OA  Sne.  •  F^po.ud  !o  11  v  ear'>  if  Dr.  E'"/iro-'m»r  *  Upon  Rue.*  pc  s *' 
•■O  Honr  ld  £ ’' ^if  U''meut , 
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FRACTIONAL 


PRINCIPAL  CONCLUSIONS 


1.  The  Initiation  Impact  Energy  (Unnotched  Charpy  Test)  Increases  with  Increasing  Values 
of  Interlaminar  Shear  Strength  for  Glass-Epoxy  and  Glass-Polyester  Laminates. 

Z  It  has  been  found  that  a  critical  value  of  interlaminar  shear  strength  (ILSS)  can  be 

defined.  Below  this  value  the  impact  propagation  energy  increases  with  decreasing  ILSS. 
Above  this  value,  the  propagation  energy  increases  with  increasing  ILSS.  Epoxy-glass 
laminates  are  always  above  this  critical  value  of  ILSS  while  polyester  laminates  can  be 
either  above  or  below  depending  upon  the  glass  surface  treatment.  Carbon-epoxy  laminates 
will  also  always  be  above  this  critical  value. 

3.  A  study  on  the  effects  of  1 1  years  of  exposure  to  dry  and  wet  (full  immersion  in  sea 

water  and  distilled  water)  environments  for  filament  wound  carbon  fiber  reinforced  epoxy 
NOL  rings  has  shown: 

3.1  The  11  years  immersion  in  wet  environments  considerably  affected  the  initial 
moisture  weight  gain  and  the  diffusion  coefficients  upon  reexposure  to  humid 
environments. 

3.2  A  considerable  shear  strength  reduction  (current  and  residual)  was  found  for 
carbon-epoxy  specimens  having  1 1  years  of  humid  exposure. 

3.3  The  shear  strength  degradation  has  been  qualitatively  correlated  and  interpreted 
with  the  moisture  absorption  characteristics. 

3.4  No  irreversible  damage  in  the  fiber  properties  was  observed  as  concluded  from 
flexural  strength  results. 

3.5  The  shear  strength  of  graphite-epoxy  specimens  exposed  to  humid  conditions 
was  greater  than  the  shear  strength  of  dry  specimens. 

4.  In  the  case  of  Thomel  300-epoxy  composites  the  following  can  be  concluded: 

4.1  The  equilibrium  content  for  complete  immersion  in  water  is  greater  than  for 
98%  relative  humidity. 

4.2  The  quality  of  the  exposed  surface  has  an  influence  on  moisture  diffusion  rate. 

4.3  The  exposure  temperature  has  an  effect  on  equilibrium  moisture  content. 

4.4  The  moisture  diffusion  coefficient  parallel  to  the  fiber  axes  is  at  least  3  times 
greater  than  the  coefficient  perpendicular  to  the  fibers. 

4.5  The  initial  moisture  absorption  for  a  unidirectional  lamina  subjected  to  a 
transverse  tensile  stress  (50%  of  ultimate  strength)  was  at  least  2  times  greater 
than  for  an  unstressed  specimen. 
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DURABILITY  OF  FIBER 


REINFORCED  RESIN  MATRIX  COMPOSITES 

JAMES  M.  WHITNEY 
AIR  FORCE  MATERIALS  LABORATORY 

I  n-House  Work  Unit 


DURABILITY  OF  COMPOSITES 

OBJECTIVE:  TO  DEVELOP  ANALYTICAL  AND  EXPERIMENTAL  METHODOLOGY 
FOR  PREDICTING  FAILURE  AND  FATIGUE  LIFE  IN  CONJUNCTION  WITH 
ENVIRONMENTAL  EFFECTS.  CONSIDERABLE  EMPHASIS  IN  THE  NEAR 
TERM  IS  PLACED  ON  ENVIRONMENTAL  EFFECTS. 

APPROACH:  ANALYTICAL  MODELS  BASED  ON  THEORETICAL  CONSIDERA¬ 
TIONS  AND  EXPERIMENTAL  OBSERVATIONS.  INTERDISCIPLINARY 
APPROACH  BETWEEN  MECHANICS  AND  MATERIALS  SCIENCE.  EMPHASIS 
ALSO  PLACED  ON  DEVELOPING  HIGH  QUALITY  EXPERIMENTAL  METHODS 
AND  PROCEDURES. 
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DURABILITY  OF  COMPOSITES 


CURRENT  ACTIVITIES 

I.  FAILURE  AND  FATIGUE 

A.  STRESS  CONCENTRATIONS 

B.  FATIGUE  OF  UNI  DIRECTIONAL  COMPOSITES 

IK  ENVIRONMENTAL  EFFECTS 

A.  DIFFUSION  STUDIES  FOR  ACCELERATED  TESTING,  INCLUDING  THERMAL  SPIKE 
EFFECTS 

B.  EXPANSIONAL  STRAINS  IN  COMPOSITE  LAMINATES 

C.  EFFECT  ON  MECHANICAL  PROPERTIES  OF  UNIDIRECTIONAL  AND  LAMINATED 
COMPOSITES 

D.  EFFECT  ON  COMPRESSION  LOADING 

E.  EFFECT  ON  NOTCH  STRENGTH 


TOPICS  OF  DISCUSSION 
STRESS  CONCENTRATIONS  -  ANGLE  CRACKS 

DIFFUSION  STUDIES  FOR  ACCELERATED  TESTING,  INCLUDING 
THERMAL  SPIKE  EFFECTS 

EXPANSIONAL  STRAINS  IN  COMPOSITE  LAMINATES 

EFFECT  OF  ENVIRONMENT  ON  MECHANICAL  PROPERTIES  OF 
UNIDIRECTIONAL  AND  LAMINATED  COMPOSITES 
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Theory  versus  experiment  -  average  stress  criterion. 
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i. 


STRESS  CONCENTRATIONS  -  ANGLE  CRACKS 


^ - - -  r  — - — 

r 

/\ 

Angle  crack  specimen,  W  ■  1".  =  0.25",  0.4",  6  •  15°,  30°,  45°,  60°,  75°,  90°, 

T300  /  5208. 
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mA 


Comparison  of  angle  crack  to  normal  crack,  average  stress  criterion. 


Least  squares  fit  for  c^ 
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(Oq  sin0 )/  W 

Comparison  of  angle  crack  to  normal  crack,  average  stress  criterion. 


Least  squares  fit  for  c^ 
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Angle  crack  compared  to  normal  crack  for  plexiglass. 


STRESS  -  CONCENTRATIONS  -  ANGLE  CRACKS 


CONCLUSION 

AN  ANGLE  CRACK  IN  A  LAMINATED  COMPOSITE  CAN  BE  TREATED 
AS  AN  EQUIVALENT  NORMAL  CRACK  HAVING  A  LENGTH  EQUAL  TO 
THE  PROJECTION  OF  THE  ANGLE  CRACK  ON  TO  THE  NORMAL 
DIRECTION. 
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DIFFUSION  STUDIES  FOR  ACCELERATED  TESTING 


Kelley-Bueche  Theory  for  glass  temperature  as  a  function  of  moisture  content. 


COMPOSITE 
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Edge  adjustment  factors  for  one-dimensional  diffusion  in  laminates  with 
quasi -isotropic  expansional  properties. 


TIME  -  DEPENDENT  RESPONSE  OF  A  RESIN  -  SLAB  SUBJECTED 
TO  A  THERMAL  SPIKE 


The  stress  response  at  the  boundary,  (z=h,t)  for  a  'Thermal-Spike"  and  the  Ensuing 
Fluctuation  in  Moisture  (case  1)  and  for  the  Same  Thermal  Spike  with  a  Fixed  Level  of 
Moisture  (case  2),  vs.  time. 


The  Stroee-Reaponee  at  the  Boundary,  »(*  Due  to  the 

•'Thermal  Spike"  and  with  w  *=0  for  aU  Timea. 


DIFFUSION  STUDIES  FOR  ACCELERATED  TESTING 


CONCLUSIONS 

1.  MOISTURE  LOWERS  THE  GLASS  TRANSITION  TEMPERATURE  OF  EPOXY  RESINS. 

THI S  ABSORPTI  ON  APPARENTLY  DOES  NOT  LEAD  TO  ANY  S 1  GNl  FI  CANT  CHEMI  CAL 
DEGRADATION,  ALTHOUGH  MECHANICAL  DAMAGE  CAN  OCCUR  IN  THE  FORM  OF 
CRAZING  OR  CRACKING. 

2.  FICK'S  LAW  ADEQUATELY  PREDICTS  TOTAL  MOISTURE  GAIN  IN  RESINS  AND 
COMPOSITES. 

3.  TWO  AND  THREE  DIMENSIONAL  DIFFUSION  CAN  BE  ESTIMATED  BY  USING  A 
MODIFIED  ONE  DIMENSIONAL  APPROXIMATION  FOR  A  WIDE  RANGE  OF  LAMINATE 
DIMENSIONS. 

4.  PRESENT  DATA  INDICATES  THAT  THE  DIFFUSION  PROCESS  IN  COMPOSITES  CAN 
INDUCE  MATRIX  CRAZING.  THUS  IT  IS  IMPORTANT  THAT  ACCELERATED  MOISTURE 
CONDITIONING  TECHNI QUES  ARE  NOT  USED  WHICH  ARE  MORE  SEVERE  THAN  REAL 
LIFE  CONDITIONS. 

5.  PRESENT  DATA  INDICATES  THAT  THERMAL  SPIKES  ABOVE  THE  GLASS 
TRANSITION  TEMPERATURE  IN  CONJUNCTION  WITH  RAPID  COOL  DOWN 
RATES  CAN  INDUCE  MATRIX  CRACKING. 


(%) 


III 


EXPANSIONAL  STRAINS  IN  COMPOSITE  LAMINATES 


T,  -c 


Transverse  coefficient  of  thermal  expansion. 
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STRAIN, 


EXPANSIONS  STRAINS  IN  COMPOSITE  LAMINATES 

CONCLUSIONS 

1.  MOISTURE  ABSORPTION  CAN  INDUCE  SIGNIFICANT  DIMENSIONAL  CHANGES 
IN  GRAPHITE  /  EPOXY  COMPOSITES. 

2,  CURRENT  DATA  INDICATES  THAT  DIMENSIONAL  CHANGES  INDUCED  IN  COMPOSITE 
LAMI  NATES  BYTEMPERATURE  AND  MOISTURE  CAN  BE  PREDICTED  BY  CLASSICAL 
LAMINATED  PLATE  THEORY. 
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V. 


EFFECT  OF  ENVIRONMENT  ON 
MECHANICAL  PROPERTIES 


Moisture  gradient  through-the-thickness  of  laminate  after  various  drying-out  times. 


2 


3  4 

STRAIN  {%) 


5 


6 


^trpcc-ctrain  curves  for  neat  resin. 


Stress-strain  curves  for  transverse  tension  of  a  unidirectional  composite. 
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INPLANE  SHEAR  STRESS  - 


Inplane  shear  stress-strain  curves  for  unidirectional  composite. 
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LOG 


T{°F) 

Flex  strength  for  unidirectional  composite  as  a  function  of  temperature  and  moisture. 


EFFECT  OF  ENVIRONMENT  ON  MECHANICAL  PROPERTIES 
OF  UNIDIRECTIONAL  AND  LAMINATED  COMPOSITES 

CONCLUSIONS 

1.  PARTICULAR  CARE  MUST  BE  EXERCISED  IN  EXPERIMENTAL  METHODOLOGY  WHEN 
PERFORMING  MECHANICAL  TESTS  ABOVE  THE  WET  GLASS  TRANSITION  TEMPERA¬ 
TURE  OF  THE  COMPOSITE. 

2.  MATRIX  DOMINATED  COMPOSITE  PROPERTIES  (E.G.  TRANSVERSE  UNIDIRECTIONAL 
TENSION.  UNI  DIRECTIONAL  SHEAR,  ETC.)  ARE  REDUCED  AT  ELEVATED  TEMPERATURE 
BY  THE  PRESENCE  OF  MOI STURE. 

3.  MOI  STURE  AT  ELEVATED  TEMPERATURES  CAN  I NDUCE  A  CHANGE  I N  FAI  LURE  MODE 
FROM  FILAMENT  DOMINATED  TO  MATRIX  DOMINATED.  CUSSIC  EXAMPLE  IS  THE 
ODEGRE  FLEX  TEST. 
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FUTURE  EMPHASIS 


1.  DIFFUSION  STUDIES  -  CONDITIONS  WHICH  PRODUCE  MATRIX  CRAZING  AND 
CRACKING. 

2.  EFFECT  OF  ENVI RONMENT  ON  LAMI  NATE  COMPRESSI  ON  STRENGTH. 

3.  DIMENSi  ONAL  CHANGES  DUE  TO  MOl  STURE  ABSORPTI  ON. 

4.  USE  OF  FLEX  TEST  FOR  EVALUATl  NG  MOl  STURE  EFFECTS  I N  COMPOS  ITES. 

5.  EFFECT  OF  ENVI  RONMENT  ON  LAMI  NATE  NOTCH  STRENGTH  ( HOLES  AND  CRACKS). 

6.  UNIDIRECTIONAL  FATIGUE  -  PURSUE  CONCEPT  OF  FATIGUE  LIMIT. 

7.  STRESS  CONCENTRATIONS  -  MATRIX  DOMINATED  LAMINATES. 
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ENVIRONMENTAL  EFFECTS  ON  THE  THERMOMECHANICAL  BEHAVIOR  OF  COMPOSITES 


T.  S.  COOKE 

SOUTHWEST  RESEARCH  INSTITUTE 


Overall  Objective 

Determine  the  influence  of  the  environment 
on  the  moisture  content  and  the  mechanical  proper¬ 
ties  of  a  resin  base  composite. 
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Specific  Goals 


Determine  the  effect  of  temperature  and  humidity 
on  water  absorption  by  the  composite . 

Investigate  the  interaction  of  the  environment 
with  the  composite  parameters  and  the  resulting  effect 
on  mechanical  properties. 


Test  Parameters 

To  carry  out  this  program,  the  following  parameters 
wili  be  investigated 

Composite  Parameters  Environmental  Parameters 

Ply  Conf  igu  ration  Temperatu  re 

Thickness  Relative  Humidity 

Volume  Fraction  Spiking 
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Moisture  Profiles  at  times  t  =  100,  200,  300  secs. 
Tor  the '*Ramp- Temperature*'  ease.  rf^or» 


ANALYSIS  AND  DESIGN  OF  COMPOSITE 
BONDED  JOINTS  UNDER  A  DYNAMIC  TYPE  LOAD 


JACK  R.  VINSON 
R,  BYRON  PIPES 


THE  CENTER  FOR  COMPOSITE  MATERIALS  AND 
THE  DEPARTMENT  OF  MECHANICAL  AND  AEROSPACE  ENGINEERING 
UNIVERSITY  OF  DELAWARE 


OBJECTIVE;  Phase  IV 

TO  SYSTEMATICALLY  CHARACTERIZE  PROMISING 
ADHESIVES  AT  BOTH  ROOM  AND  ELEVATED  TEMPERATURE 
USING  THE  UNIOUE  SHEAR  AND  TENSILE  TESTS  DEVELOPED 
AT  DELAWARE. 


OBJECTIVE-  Phase  V 

INVESTIGATING  THE  EFFECTS  OF  COMBINED  HIGH 
TEMPERATURE  AND  HIGH  HUMIDITY  (HYGROTHERMAL 
EFFECTS)  UPON  BONDED  JOINTS  BETWEEN  STRUCTURAL 
COMPONENTS  OF  COMPOSITE  MATERIALS. 
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0 


HOBAY  MONDUR  75 
HOBAY  MONDUR  MR 
MOBAY  MONDUR  PF 
CONAP  DPAD 
CONAP  AD-3 
3M  AF-30 
3M  W 
3M  m 
3M  AF-55 


D- MOBAY  MONDUR  MR 
I  -  3M  AF-30 


*  SHEAR  STRAIN  (IN./ IN.)  * 


9 


MAfiUFACTURER 

ULT.  SHEAR  STRESS (PS I) /STD.  DEV. 
1989  /  7.5 

2182  /  92 
2273  /  67 
W59  /  86 
2061  /  98 
2m  /  239 
2A05  /  107 
2171  /  37 


UNIV.  OF  DELAWARE 
ULT.  SHEAR  STRESS (PS I) /STD,  DEV. 

1622  /  43 
2512  /  189 
2379  /  74 


2080  /  125 
2560  /  63 
2175  /  282 
2239  /  196 


MANUFACTURER 

ULT.  SHEAR  STRAIN/STD.  DEV. 


MOBAY  MONDUR  75 

5.137  /  .932 

MOBAY  MONDUR  MR 

7.976  /  3.741 

MOBAY  MONDUR  PF 

7.037  /  1,484 

CONAP  DPAD 

1.911  /  .431 

CONAP  AD-3 

1.902  /  .3 

3M  AF-30 

>2.451  /  - 

3M  147 

1.858  /  .831 

3M  143 

1.681  /  .663 

3H  AF-55 

MANUFACTURER 

SHEAR  MODULUS  (PSD  /STD 

MOBAY  MONDUR  75 

813  /  152 

MOBAY  MONDUR  MR 

750  /  330 

MOBAY  MONDUR  PF 

768  /  210 

CONAP  DPAD 

992  /  209 

CONAP  AD-3 

2614  /  758 

3M  AF-30 

2351  /  213 

3M  147 

2566  /  869 

3M  143 

2469  /  737 

3M  AF-55 

UiHV.  OF  DELAWARE 
ULT.  SHEAR  STRAIN/STD.  DEV. 

6.651  /  2.931 
1.808  /  .868 
>2.229  /  - 


>2.918  /  - 
1.558  /  .378 
1.392  /  .633 
1.715  /  .536 


UN  IV.  OF  DELAWARE 
SHEAR  MODULUS  (PSD  /STD.  DEV. 

630/213 
362^*  /  1516 
2965  /  586 


1522  /  242 
3251  /  1083 
3717  /  1547 
2575  /  697 
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TEST  MATRIX- PHASE  V 


3  REPETITIVETESTS 

1  ADHEREND  SYSTEM  (Graphite-Epoxy) 

2  ADHESIVE  SYSTEMS 

5  TIMES  IN  SOAK:  0,1, 3, 5, 7  Months 

3  TESTS;  Static,  HC-LL  Fatigue;  LC-HL  Fatigue 

2  TEST  ENVIRONMENTS:  RoomTemp.,  Elevated  Temp. 

laO  TEST  PIECES 


OTHER  RESEARCH  COMPLETED 


,  BUCKLING  OF  COMPOSITE  PLATES  DUETO  IN-PLANE  LOADS 
FOR  MANY  BOUNDARY  CONDITIONS 


INTERLAMINAR  STRESSES  IN  COMPOSITE  CYLINDRICAL 
SHELLS 


•  HYGROTHERMAL  EFFECTS  IN  COMPOS ITE  LAMINATES 


PLANE  STRAIN  FINITE  ELEMENT  INCLUDING  HYGROTHERMAL 
•  EFFEQS 


•  PARAMETRIC  STUDIES  OF  STRUQURAL  ADHESIVE  JOINTS 
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FORCES  ON  A  TOO  LAMINA  COMPOSITE  SHELL 


BORON  EPOXY 

=  32,5  X  #  psi 

H=  0.25" 

E22  =  1-84  X  10^  PSI 

R  «  12.5" 

^^2  =  0.256 

L  =  50" 

^22  =  ^23  “  0,642  X  10^  PSI 

p  =  100  PSI 

G23  =  0.361  X  #  PSI 

BENDING  BOUNDARY  LAYER  IN  MERIDIONAL  DIRECTION 


ISOTROPIC  MATERIAL 

ORTHOTROPIC  MATERIAL 

LAMINATED  SHELL 
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— - 

tuA*AP60 

— 

S«Vi 

**  *  ^s* 

Single  Surface  Absorption  -  o,.  Stress  Disrrihurion. 


Single  Surface  Absorption  -  Tj^,.  Stress  Distribution. 


OTHER  RESEARCH  IN  PROGRESS 


ANALYSIS  OF  BONDED  JOINTS  OF  SEVERAL  CONFIGURATIONS 
•  USING  THE  ANALYTICAL  FINITE  ELEMENT  INCLUDING 
HYGROTHERMAL  EFFECTS 


ANALYSIS  OF  SINGLE  LAP  JOINTS  WITH  LINEARLY 
*  TAPERED  ADHERENDS 


NEW  TECHNIQUES  OF  SOLUTION  FOR  COMPOSITE  PLATES 
AND  SHELLS  USING  PERTURBATION  TECHNIQUES 
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AIR  FORCE  OFFICE  OF  SCIENTIFIC  RESEARCH 


AND  THE 

CENTER  FOR  COMPOSITE  MATERIALS 
UNIVERSITY  OF  DELAWARE 

WORKSHOP 

ON 

THE  EFFECTS  OF  RELATIVE  HUMIDITY  AND  ELEVATED 
TEMPERATURE  ON  COMPOSITE  STRUCTURES 
MARCH  30-31,  1976 
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EFFECT  OF  COMPRESSION  ON 
FATIGUE  OF  GRAPH ITE/ePOXY  COMPOSITES 


PRINCIPAL  INVESTIGATOR;  J.  T.  RYDER 
PROGRAM  MANAGER;  E.  K.  WALKER 


LOCKHEED-CALIFORNIA  COMPANY 


PURPOSE 

BUILD  STATISTICAL  DATA  BASE 

EVALUATE  WEAROUT  MODEL 

EVALUATE  EFFECT  OF  COMPRESS  ION  FATIGUE 
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PROGRAM  DEFINITION 


TWO  LAMINATES 

LAMINATE  I:  25%  0°  FIBERS 
16  PLIES 

(0IH5l90H3j90im)^ 
QUASI- ISOTROPIC 

LAMINATE  2;  66. 7%  0°  FIBERS 

24  PLIES 


PROGRAM  DEFINITION  (CONT'D.) 


TESTING 


STATIC  TENSION.  COMPRESSION 

S-N  CURVE,  T-T,  C-T 

EVALUATE  FATIGUE  SCATTER 

MINIMUM  TWENTY  TESTS  AT  EACH  STRESS  LEVEL 

RESIDUAL  STRENGTH 

FATIGUE  CYCLE  40  COUPONS  AT  STRESS  LEVELS  CHOSEN 
FAIL  HALF  TENS  ION  &  COMPRESSION 
ORIGINAL  STRENGTH  DISTRIBUTION 


PROGRAM  DEFINITION  (CONT'D) 

MATER  I  Ab  (934/T300) 

CURE  CYCLE:  350®F  CURE  FOR  2  HOURS 

CONDITIONING:  72  HOURS  AT  72  +  2°F,  40  +  I0%R.H. 

TEST  ENV I RONMENT:  72  +  2®F,  40  +  10%  R.  H. 
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COUPON  TYPE 


PARALLEL  SIDED 
5-1/2  IN.  CENTER  SECTION 
NON-TAPERED  END  TABS 


■“1 - 3.  OtO.  1 - -  5. 50+0. 1  - 1- 

-  8. 50+  0. 1  - - - 

<3 - 11.50+0.01  - - 1> 
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FABRICATION 


VARIATION  IN  AREA;  +  1.5%  WITHIN  COUPON 

VARIATION  IN  AREA;  +4%  FOR  ALL  COUPONS 

PERCENT  VOLUME  OF  FIBER; 

LAMINATE  I;  60.3-  66.4 

LAMINATE  2;  60.0-  63.1 

VOIDS;  «l/2% 


COMPRESS  ION  TEST  IMJ-601-23B 
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tensile  tests 


LAMINATE  1 


LAMINATE  2 


ULTIMATE 

ULTIMATE 

APPARENT  MODULUS 

STRESS, 

STRAIN 

OF  ELASTICITY, 

^ult’ 

®ult’ 

E* 

6 

ksi 

in. /in. 

psi  X  10 

69.2 

0.0096 

7.64 

+  6.2,  9i 

+0.0011,  Ufo 

+0.54, 

-  7.2, 

-0.0010,  10^ 

-0.69,  95J 

l4l.7  0.0092 

+20.1,  l4^  +0.0009,  10^ 

-20.3,  M  -0.0027,  29/a 


15.4 

+2.9,  19^ 
-1.9,  12fo 


TENSILE  TESTS  (CONT’D) 


LAMUTATE  1 


LAMINATE  2 


k 

e 

V 

■  R  &  FIT 

3  PARAMETER 

24.18 

-0.0753 

70.50 

EQUAL 

2  PARAMETER 

24.12 

0 

70.57 

3  PARAMETER 

18.87 

-0.4304 

144.4 

EQUAL 

2  PARAMETER 

18.76 

0 

144.8 

P  (x)  -  exp 

L  -  ^  J  ^ 

FOR  2  PARAMBIER,  or  =  k,  F  = 
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COMPRESSION  TEST  IMJ-5SO-1A 
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CCMPRESSION  TESTS  (COMT’D) 


ULTIMATE 

ULTIMATE 

APPARENT  MODULUS 

STRESS, 

STRAIN, 

OF  elastichy, 

°ult» 

®ult’ 

'a. 

ksi 

In. /in. 

psi  X  10"" 

69.7 

0.0110 

6.95 

LAMItTATE  1 

+11.8,  17^ 

+0.0022,  20^ 

•K).29, 

H 

-12.1,  17^5 

-0.0025. 

-o;83. 

-izi 

114.2 

■■■■ 

LW'IINATE  2 

±15.0,  I3i 

+0.8, 

li 

-10.2,  Si 

■■ 

-0.5, 

hi 

COMPRESSIOT  TESTS  (COKT'D) 


k 

e 

V 

R  &  FIT 

LAMINATE  1 

3  PARAMETER 

12.09 

-0.4572 

72.19 

EQUAL 

11.94 

0 

72.53 

LAMINATE  2 

3  PARAMETER 

25.01 

-0.2592 

115.82 

EQUAL 

2  PARAMETER 

24.91 

0 

116.07 
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MAXIMUM  STRESS,  o-.  KSI 
max 


S-N  CURVE  SCAN 


THREE  COUPONS  AT  EACH  STRESS  LEVEL 
CHOSE  CONSTANT,  VARIED 

ALLOWED  COMPARISON  BY  ^  AND 
TEST  FREQUENCY,  10  Hz 


LAM  I  NATE  1 


CYCLES  TO  FAILURE,  Np 
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A  ® 

G  Q 

AO  B 

O  Q 

A  G>  D 

A  O 

A  O 


CYCLES  TO  FAILURE, 


LAMINATE  1 
FATIGUE  TEST  RESULTS 

O 


O 


LAMINATE  1 
FATIGUE  TEST  RESULTS 

A  A  o;^in“-16ksl 

O 

O 

A 

O 

X2  A 

O 

Q  A 
© 
a 

o 


CYCLES  TO  FAILURE,  N. 
82  *■ 


Maximum  Stress,  ksi 


180  r- 


TEKSlorr  -  TEH3I0N 
FATIGUE  TEST  RESULTS 
LAMTOTE  2 


min 


0  kGi 


G 


© 


CG 

0  G . » 


(E& — ► 

3 

0-> 


10' 


10- 


10 


10' 


10 


10' 


Cj'-cles  to  Failure,  N 


FATIGUE  RESULTS 

CHOSE  SPECIFIC  LEVELS 

rndA 

TWENTY  COUPON /  C  LEVEL 

PICKED  SPECIFIC  STRESS  LEVELS  FOR  RESIDUAL 
STRENGTH  STUDY 
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MAX  I  MuX  STRESS,  max,  ksi  MAX  lAVJM  STRESS,  CTmsx,  ksi 
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5%  ProUhility 
of  Survival 


95%  Probability' 
of  Survival 
? 

Tension-Tension 
FatioueTest  Results 
Laminate  1 

O'  •.  '^Oksi 
min 


IS 

I  I 

20 


Note:  Numbers  refer  to 
number  of  Test  Points 


INDICATES  NO  FAILURE 


CYCLES  TO  FAILURE,  N, 


-V  X  5%  Probability 

\  of  Survival 

< 


95%  Probability^ 
of  Survival 


1\ 

\ 

20 


Compression-Tension 

Fatigue  Test  Results 

Laminate  1 

CT  .  - -16  ksi 
^  min 

f  “  10 II 


Note;  Numbers  refer 
to  Number  of 
Test  Points 


CYCLES  TO  FAILURE,  Np 
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'ollure,  N,,,  Million 


Maximum  Stress,  a  ,  ksi 


80 


70 


lAI'IIKATE  I 

caMPKi‘:ss]:oN  -  aiiKsiav 

FATIGUi':  Tt::JT  RESULTS 

o  .  =  -  Id  ksi 

min 

f  -  10  Hz 


6o 


6  20 


^0(^20 


40 


30 


20 


(5  20 
2)0  3 
^  3  00 

O  GE)  O 


00025  <S)  00  0  19  © 


O  O 
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NOTE:  Numbers  indicate  Total 
Number  of  Coupons  Tested 
at  that  Stress  Level. 


0 

0 


2  4  6 

Cycles  to  Failure,  Np,  Millions 


WEI  BULL  PARAMETERS 


LAMINATE  1 


Stress 

ksi 

k 

e 

V 

max 

min 

58 

0 

1.988 

-  2,030 

4,899 

50 

0 

4564 

-  20,711 

64,621 

42 

0 

1.846 

-280,781 

1,900,898 

58 

-16 

2.409 

-  1,015 

1,979 

50 

-16 

2.796 

-  4,973 

10,964 

42 

-16 

2.915 

-  3,0^0 

33,806 

26 

-16 

1.214 

+  74,246 

2,892,580 

Maxiaum  Stress,  cr  ,  ksi 


Probability,  P;  Survival 


.0001 

.0005 

.001 

.002 

.005 

.01 

.02 

.05 

.10 

.20 

.30 

.40 

.50 

.60 

.70 

.80 

.90 

.95 

.98 

.99 

.998 

.999 

.9999 


Failure  Stress,  a,  ksi 


RESIDUAL  STRENGTH  STUDY  -  LAMINATE  1 

COUPONS  FATIGUED  TO  P  AT  EACH  (T  LEVEL 

FAILED  COUPONS  REPLACED 

40  COUPONS/  c::r  LEVEL,  1/2  T  AND  1/2  C 

FAILURE  MODE  SAME  AS  PREVIOUS  T  AND  C, 
DISSIMILAR  TO  FATIGUE 


RESIDUAL  STRENGTH  STUDY 
LAMINATE  1 

FATIGUE  STRESS  LOWEST  TENSILE  FAILURE  MODE,  AVERAGE  TENSILE  STRENGTH 
LEVEL,  KSI  STRESS,  KSl  KSI  KSI 


OTO  42 

49.2 

64.1 

62.4 

OTO  50 

54.6 

66.7 

65.1 

-16  TO  42 

59.0 

68.2 

67.2 

-16  TO  50 

58.3 

69.0 

67.7 

STATIC  TENS  ILf  TEST 

62.0 

70.3 

69.2 

COUPONS  FAILED  IN  FATIGUE  NOT  INCLUDED  -  FAILURE  MODES  DIFFER 
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SUMMARY  OF  RESIDUAL  STRENGTH  TESTS  -  LAMINATE  1 
TENSILE 


Fatigue  Stress  Level, 
ksi 

Weibull  Parameters 

Average 

Ultimate 

Stress, 

Average 

Ultimate 

Strain, 

®avg’ 
in. /in. 

Average 

Initial 

Slope, 

E  , 
avg' 

k 

e 

V 

a  ,ksi 

avg' 

.  _ -6 
psi  X  10 

Tensile  Test  Results 
Before  Fatigue 

24.18 
24 . 12 

-0.0753 

0 

70.50 

70.57 

69.2 

0.0096 

7.64 

0  to  1(2 

12.12 

11.88 

-0.6276 

0 

64.52 

65.04 

62.4 

•0.0090 

6,84 

0  to  50 

12^85 

12.78 

-0.1688 

0 

67.31 

64.45 

65.1 

0.0093 

7.01 

-16  to  42 

18.51 

18.48 

-0.0608 

0 

68.43 

68.49 

67.2 

0.0094 

7.17 

-l6  to  50 

18.34 

18.27 

-0.1326 

0 

_ 

69.26 

69.38 

67.7 

0.0095 

7.13 

SUMMARY  OF  RESIDUAL  STRENGTH  TESTS  -  LAMINATE  1 
COMPRESSIVE 


Fatigue  Stress  Levels, 
ksi 

Weibull  Parameters 

Average 

Ultimate 

Stress, 

a  ,ksi 

avg 

Average 

Ultimate 

Strain, 

^avg’ 
in. /in. 

Average 

Initial 

Slope , 

E  , 
avg’ 

psi  X  10^ 

k 

e 

V 

Compressive  Test  Results 
Before  Fatigue 

12.09 

11.94 

-0.4572 

0 

72.19 

72.53 

69.7 

0.0110 

6.95 

0  to  42 

7.22 

7.07 

-0.6574 

0 

60.15 

60.52 

57.2 

0.0098 

6.37 

0  to  50 

13.02 

12.94 

-0.2053 

0 

65.62 

65.79 

63.7 

0.0111 

6.64 

-16  to  42 

12.54 

12.46 

-0.2387 

0 

67.25 

67.45 

65.3 

0.0108 

6.66 

-16  to  50 

14.74 

14.52 

-0.5418 

0 

moo 

d  d 

68.6 

0,0110 

6.66 
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RESIDUAL  STRENGTH  STUDY  -  LAMINATE  2 

T-T  COUPONS  CYCLED  TO  10^  AT  0  TO  100  ksl 
C-T  COUPONS  CYCLED  TO  !0^  AT  -30  TO  70  ksi 
40  COUPONS  AT  EACH  CONDITION,  1/2  T  &  1/2  C 
FAILED  COUPONS  REPLACED 
FAILURE  MODE  SAME  AS  PREVIOUS  T  &  C 


RESIDUAL  STRENGTH  STUDY  -  LAMINATE  2 


C 


TYPE  OF 
STATIC  TEST 
T  “  TENSILE 
-  COI-tPRESSICN 


T 


FATIGUE 

STRESS 

WEIBUL  PARAMETERS 

AVERAGE 

ULTIMATE 

STRESS 

AVERAGE 

ULTIMATE 

STRAIN 

e 

LEVEL 

ksi 

k 

e 

V 

^avg,  ksi 

avg, 
in. /in. 

TENSILE  TEST 
RESULTS  BEFORE 
FATIGUE 

18.87 

-0.4304 

l44.,4l 

i4i.7 

.0092 

0  to  ICO 

17.19 

-0.3926 

146.81 

143.5 

.0092 

-30  to  70 

21.57 

-0.1409 

146.42 

143.7 

.0094 

COMPRESSION 

TEST  RESULTS 
BEFORE  FATIGUE 

0  to  100 

-30  to  70 

25.02 

-0.2592 

115.62 

114.2 

.0098 

AVERAGE 

INITIAL 

SLOPE 


avg,  g 
psi  X  ,10- 


15.7 

15.3 

11.5 


LAMINATE  1 


COMPARISON  OF  ACTUAL  AND  PREDICTED  FAILURE 
RATES  AT  P  gg 

(THREE  PARAMETER  WEI  BULL  ANALYSIS) 


TOTAL  SAMPLE  SIZE  240 

FAILURES  24 

PREDICTED  FAILURES  23 


LAMINATE  1 

COMPARISON  OF  ACTUAL  AND  PREDICTED  FAILURE 
RATES  AT  P  o. 


TOTAL  SAMPLE  SIZE 
FAILURES 

PREDICTED  FAILURES 
PREDICTED  FAILURES 


240 

12 

12  (3  PARAMETER  FIT) 
19  (2  PARAMETER  FIT) 
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PROCABiLiTY  OF  FAILURE 


90 

80 

70 

60 

30 

40 

30 

20 

15 

10 

5 

4 

2 


CYCLES  TO  FAILURE,  N, 

LAMINATE  1  ^ 

PROBABILITY  OF  FAI  LURE  FOR  TENS  ION-TENS  ION 

FATIGUE  TESTS  AT  M2i<SI 
max 


CALCULATED 

SURVIVABILITY 


SAMPLE  SIZE 

EFFECT  OF  SAMPIT.  SIZE  ON  PROBABILITY 
OF  SURVIVAL  (0TO‘12  KSl) 
LAMINATE  1 
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SUMMARY 


1.  STATIC  COMP  RES  SI  ON  SCATTER  LARGER 
THAN  TENS  ION 

2.  LARGE  FATl GUE  DATA  SCATTER 

3.  COMPRESSION  LOWERS  FATIGUE  LIFE 

4  STATIC  AND  FATIGUE  FAILURE  MODES  DIFFER 

5.  EFFECTS  OF  cr  AND  /icT  UNCLEAR,  NEITHER 
CORRELATE  WELL 

6.  LAMINATES  1  AND  2  DIFFER  IN  STATIC  AND  FATIGUE 
RESPONSE 

7..  FIRST  CYCLE  FATIGUE  FAILURES  OCCUR  FOR  LAMINATE  2; 
INDICATES  LARGE  STATIC  STRENGTH  SCATTER? 

8.  FOR  LAMINATE  1,  AT  (T  LEVELS  WHERE  FATIGUE  FAILURE 
OCCURS,  WEAROUT  OCCURS 

9.  FOR  LAMINATE  2,  NO  WEAROUT  AT  ^  LEVELS  WHERE 
FATIGUE  FAILURE  OCCURS 

10„  THREE  PARAMETER  FITS  SHOULD  BE  USED  FOR  FATIGUE 

11.  SAMPLE  S IZE  TWENTY  OR  LARGER 

12,  CAN  RESULTS  BE  EXTENDED  TO  LARGE  COUPONS  OR 
PANELS? 
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COMPOSITE  STRUCTURE  OPTIMAL 
DESIGN  BY  ENERGY  METHODS 

M.  Murphy 
Norwich  University 


OBJECTIVES:  To  formulate  and  test  appropriate  optimization 

criteria  for  the  design  of  fail-safe  aerospace  com¬ 
posite  structures.  Since  design  criteria  are  based 
on  failure  modes  and  many  composites  of  interest 
in  aerospace  applications  fail  by  fracture,  a  specific 
research  objective  is  the  development  of  criteria 
predicting  fracture  onset,  fracture  path  and 
fracture  arrest. 
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ACCOMPLISH?tnNTS:  A  theory  of  fracture  has  been  developed  for  a 

quasistatic,  adiabatic  fracture  process  based  on  the  principle 
of  naximum  rate  of  entropy  production.  New  minimum  principles 
have  been  developed  from  the  Calculus  of  Variations  for 
predicting  the  crack  path  and  crack  initiation.  Specifically, 
these  developments  lead  to  a  theory  of  crack  movement,  a 
numerical  procedure  for  crack  initiation,  conditions  for 
crack  arrest,  and  conditions  for  instability  in  the  running 
crack.  The  numerical  results  for  fracture  in  isotropic  media 
are  in  agreement  with  results  of  other  investigators. 
Presently,  numerical  methods  applying  these  minimum 
principles  to  simple  composite  fracture  are  being  developed 
and  evaluated.  A  now  test  method  for  measuring  the 
environmentally  sensitive  fiber-matrix  interface  fracture 
parameters  is  under  evaluation  also. 


JM  IMCHMMT  or  CRACK  LEHCTH 
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da  da 


U(a,a)  ^  U(Gi,a)  ♦  /J(C(a,C)  -  G(a,C)>ilC 
U(o,a)  <  lJ(a,a) 
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SUMMARY  OP  CONCEPTS 


‘  t  t  f  ' 


n  T  i  T  ' 


Crack  In  a  Teat  Spaclaen 


Crack  Path  Description. 
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Initial  Position  of  Crack 


Final  Position  of  Crack 


Initial  Position  of  Crack 
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«  c 

Mesh  Used  and  Crack  Movement  in  the  Test*  Specimen. - 
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crack 


PBHHH 

HHBiH 

smaBmmssa 


asaBBBaBHB 

HHBHH 

HHHm 


Initial  Position  of  Crack 


BBBBBBBBBB 


Final  Position  of  Crack 


^Extension  Percent  error 

Kodel 

Results  5. 7x10^  -5% 

Griffith 

Solution  6x10^^ 

Energy 

Difference 

Solution  8, 1x10^^  4-35Z 


Griffith  Solution  Results 
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MATERIALS  SCIENCES  CORPORATION 


Effects  of  Environment,  and  Damping  and  Coupiing  Prqpertifs 
OF  Composite  Laminates  on  Panel  Fiuttfr 
(AFOSR  contract  F44620-76-C-0080) 


objectives 

I.  develop  epoxy-matrix  composite  laminate  constitutive  relations 

WHICH  ACCOUNT  FOR  THE  EFFECTS  OF  LAMINAE  STACKING  SEQUENCE^ 
transverse  SHEAR  DEFORMATION^  AND  TIME-DEPENDENT  BEHAVIOR  IN 
A  FORM  THAT  IS  APPLICABLE  TO  STRUCTURAL  DYNAMICS. 

II.  EXTEND  CONSTITUTIVE  RELATIONS  UNDER  "l."  TO  CONSIDER  EFFECTS 

OF  MOISTURE  ABSORPTION  AND  ELEVATED  TEMPERATURE  ON  EPOXY-MATRIX 
ADVANCED  COMPOSITES. 

III.  EVALUATE  DAMPING  CHARACTERISTICS  OF  COMPOSITE  PANELS  USING 
CONSTITUTIVE  THEORY  DEVELOPED. 

IV.  PERFORM  PANEL  FLUTTER  ANALYSIS  USING  LINEAR  PISTON  THEORY 

AERODYNAMICS  AND  THE  CONSTITUTIVE  THEORY  DEVELOPED.  QUANTI¬ 
TATIVELY  EVALUATE  COMPOSITE  MATERIAL  AND  DESIGN  FACTORS 
WHICH  HAVE  MAJOR  INFLUENCE  ON  PANEL  FLUTTER. 
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amplitude  „  (in.) 


I.  LAMINATE  CONSTITUTIVE  RELATIONS 


OBSlfiyAIJONS 

•  DAMPING  IS  AN  INHERENT  PROPERTY  IN  FIBER-REINFORCED 
EPOXY-MATRIX  COMPOSITES. 

•  BENEFICIAL  EFFECT  IS  THE  ATTENUATION  OF  PEAKS  AT  RESONANCE 
(particularly  at  HIGHER  FREQUENCIES). 


AMPLITUDE  OF  LATERAL  DEFLECTION  OF  A  VISCOELASTIC  PLATE  AS  A 
FUNCTION  OF  FREQUENCY 


approach 
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4K 

w*  4  K  A  ^  a  * 

K  +  (1  +  -— ) 


_  n 

"  /»/*  "  ^ 

2=t 


where  V™  and  V^  are  the  volume  fractions  of  matrix  and  fibers  in 
the  fiber  bundle,  respectively, 


:  =  ^  bV)  (1  t  (£1^ 

^  Y  -  1  1  + 


iV^)h 


-  V*)  (1  +  ~  (V^)^) 

^  ^  1  +  y6^ 


-  3V*(V™)^  (B"')^ 


Y  ’ 


s'"  =  - ^;;^-  and 

2g” 

1  +  ^ 
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2g; 

1  -4- 

K 
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OBTAINED  FROM  LAMINA  EFFECTIVE  COMPLEX  MODULI  AND 
ASSUMED  DISPLACEMENT  FIELD 

h/2  ^ 

{A .  .  / B .  . , D .  , )  =  P  -i  (^#2,2  )  dz  (i#j 

ij  ij  ij  </  ij 

-h/2 
C . 

°ij  “  ^ij  ■  ^3a  (plane  STRESS  REDUCED  STIFFNESSES) 

h/2 

^ij  '  ■'ij  /  ° 

-h/2 


C^j  =  LAMINA  COMPLEX  STIFFNESS  MATRIX 

k,  .  =  SHEAR  CORRECTION  FACTORS  FOR  TRANSVERSE 
■’  SHEAR  (real) 

CAN  BE  COMPLEX? 

A?,  (oj)  +  iA?  .  (u)  =  A?,  (id) 


1  +  i 


A^.lo.) 


LOSS  FACTOR 


TRANSVERSE  SHEAR  DEFORMATION  EFFECTS 


RATIOS  OF  OR  FOR  COMPOSITES  HIGHER  THAN 

ISOTROPIC  MATERIALS. 

"softer"  material  in  SHEAR  INCREASES  THE  IMPORTANCE 
OF  SHEAR  DEFORMATION. 

SINCE  STRUCTURAL  DAMPING  IN  VIBRATING  LAMINATED  BEAM/ 
PLATES  ARISES  PRIMARILY  DUE  TO  INTERLAMINAR  SHEAR^ 

THE  COUPLED  EFFECT  OF  SHEAR  DEFORMATION  AND  DAMPING 
MAY  BE  IMPORTANT. 

UTILIZE  A  SHEAR  DEFORMATION  THEORY  WHICH  CONSIDERS 
SHEAR  CORRECTION  FACTORS  FOR  OBTAINING  LAMINATE  DYNAMIC 
VISCOELASTIC  CONSTITUTIVE  EQUATIONS. 

DETERMINE  THE  SHEAR  CORRECTION  FACTORS  AS  A  FUNCTION  OF 
STACKING  SEQUENCE  AND  LAMINATE  CONSTRUCTION  FROM 
"exact"  ELASTIC  ANALYSIS. 


SHEAR  CORRECTION  FACTOR 
VS.  NO.  OF  LAYERS 


n.  moisture/temperature  effects 


OBSERVATIONS 

•  EFFECT  OF  MOISTURE/tEMPERATURE  IS  TO  ALTER  LAMINA/lAMINATE 
COMPLEX  MODULI  (PARTICULARLY  THOSE  WHICH  ARE  MATRIX  CONTROLLED). 

’  PLASTICIZING  EFFECT  OF  MOISTURE/TEMPERATURE  REDUCES  THE  STORAGE 
MODULUS  AND  INCREASES  THE  LOSS  FACTOR  (FOR  EXAMPLE^  MECHANICAL 
DAMPING  INDEX  PEAKS  WHEN  THE  RELATIVE  RIGIDITY  DROPS  AT 
T  “  Tg  FOR  AN  EPOXY;  HOWEVER^  DATA  FOR  COMPLEX  DYNAMIC  PROPER¬ 
TIES  OF  EPOXIES/lAMINA  FOR  DIFFERENT  MOISTURE  CONTENTS  AND 
TEMPERATURES  ARE  LACKING). 


APPFDACH 

•  DETERMINE  LAMINA  DIFFUSIVITIES  FROM  CONSTITUENT  DIFFUSIVITIES 
(using  conductivity  analogy  and  composite  CYLINDER  ASSEMBLAGE 

results) . 

.  DETERMINE  THE  THROUGH-THE-THICKNESS  TRANSIENT/sTEADY-STATE 
MOISTURE  DISTRIBUTION  IN  A  LAMINATE  FROM  THE  1-D  DIFFUSION 
EQUATION. 

•  FOR  THE  transient/steady  STATES^  DETERMINE  EFFECT  OF  MOISTURE 
DISTRIBUTION  ON  LAMINATE  PROPERTIES  (USE  APPROPRIATE  KNOCK¬ 
DOWN  FACTORS  FOR  STORAGE  AND  LOSS  MODULI;  ALSO,  CONSIDER  THE 
SITUATION  WHEN  SYMMETRIC  LAMINATES  MAY  BEHAVE  AS  ASYMf-IETRIC 
LAMINATES  WHEN  MOISTURE  DISTRIBUTION  IS  NOT  SYMMETRIC). 
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III.  DAMPING  CHARACTERISTICS  OF  COMPOSITE  PANELS 


AREfiQAQl 

•  INTEGRATE  THE  ANALYSES  FOR  THE  DETERMINATION  OF  (a)  DYNAMIC 
VISCOELASTIC  PROPERTIES  OF  LAMINATES  FROM  CONSTITUENT 
PROPERTIES^  (b)  shear  (ELASTIC  CASE)  CORRECTION  FACTORS 
FOR  DIFFERENT  LAMINATES.  AND  (c)  MOISTURE/TEMPERATURE 
EFFECTS  INTO  A  SINGLE  PREDICTIVE  TOOL  (COMPUTER  PROGRAM). 

•  COMPARE  THE  LAMINATE  DYNAMIC  COMPLEX  MODULI  AS  A  FUNCTION 
OF  FREQUENCY  OBTAINED  FROM  THE  COMPUTER  PROGRAM  WITH 
EXPERIMENTAL  DATA. 

•  DETERMINE  THE  DYNAMIC  RESPONSE  CHARACTERISTICS  (PEAK 
AMPLITUDE  VS.  FREQUENCY;  NO  AERODYNAMIC  LOADING)  FOR 
LAMINATES. 


Ill 


IV,  PANEL  FLUTTER 


APPROACH 

•  SOLVE  THE  CYLINDRICAL  BENDING  AND  A  FINITE  SIZE  PUTE 
FLUTTER  PROBLEM  USING  GALERKIN  PROCEDURE. 

•  INVESTIGATE  THE  EFFECTS  OF; 

.  STRUCTURAL  DAMPING  (DYNAMIC  VISCOELASTIC  LAMINATE 

properties) 

•  AERODYNAMIC  DAMPING 

•  PRESTRESS 

•  SHEAR  DEFORf'V\TION  WITH  SHEAR  CORRECTION  FACTOR  OBTAINED 
AS  A  FUNCTION  OF  STACKING  SEQUENCE  AND  LAMINATE 
CONSTRUCTION. 

•  BENDING/EXTENSION^  bending/twisting  COUPLING 

•  moisture/temperature 

•  QUANTITATIVELY  EVALUATE  FACTORS  WHICH  HAVE  MAJOR  INFLUENCE 
ON  PANEL  FLUTTER. 


PANEL  CONFIGURATION  AND  COORDINATE  SYSTEM 

VARIOUS  PLATE  GEOMETRIES 
y  S.S,  (x  =  0,  a) 
y  -►<»,  clamped  (x  =  0)  ;  free  (x  =  a) 
S.S.  (x  =  0,a) ,  (y  =  0,b) 
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•  SHEAR  DEFORMATION  THEORY  FOR  LAMINATED  PLATES  (AFTER 
WHITNEY  S  PAGANO)  WITH  COMPLEX,  PLANE  STRESS  REDUCED 
STIFFNESSES  AND  THE  SHEAR  CORRECTION  FACTORS  INTRODUCED 
DIRECTLY  INTO  THE  GOVERNING  EQUATIONS. 

•  CALCULATED  FOR  EACH  LAMINATE  FROM  THE  THICKNESS  SHEAR 
MODE  CUT-OFF  FREQUENCIES. 

.  FROM  PISTON-THEORY  AERODYNAMICS, 

AERODYNAMIC  LOAD  =  "  ^  (w,„  cos  6  +  w,„  sin  8)  -pcw,^ 

D  *  y  t 

<3  =  FREE-STREAM  DYNAMIC  PRESSURE^  j  PU^ 

U  -  FREE-STREAM  AIR  VELOCITY 

B  =  COMPRESSIBILITY  FACTOR,  i  (m  IS  THE  MACH  NUMBER) 
0  =  FLOW  ANGLE  wrt  X-AXIS 
c  =  FREE-STREAM  SPEED  OF  SOUND 
P  “  FREE-STREAM  AIR  DENSITY 

SO 

p(x,y,t)  *  “Ppanel'^'tt  AERODYNAMIC  LOAD 
"FREQUENCY  LOOPS"  IN  THE  PLANE  (FOR  FIRST  FIVE  MODEs) 
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PLANE 


1/2 

"frequency  loops"  in  the  (w) 


200 


160 


"frequency  loops"  in  the  plane 

[0/90] g  laminate 

loss  factor  =0.05  (for  ^22*^1 

~  real 

— —  IMAGINARY 


FLUTTER 


FREQUENCY  LOOPS"  IN  THE  w  PLANE 


PRELIMINARY  OBSERVATIONS 

•  CRITICAL  VALUE  OF  FLUTTER  PARAMETER  ^  IS  REDUCED  WHEN 
DAMPING  IS  CONSIDERED  (TENDS  TO  CONFIRM  EARLIER 
OBSERVATIONS  THAT  INSTABILITY  LOAD  IS  LOWERED  FOR  NON¬ 
CONSERVATIVE  SYSTEMS  WITH  DAMPING). 

.  FURTHER  STUDY  OF  THIS  PHENOMENON  REQUIRED  BY  INVESTIGATING 
BEHAVIOR  OF  THE  COMPLEX  FREQUENCIES  IN  THE  VICINITY  OF 
FLUTTER  POINTS  (V/HICH  CORRESPOND  TO  REAL  FREQUENCIES)  AND 
INCLUSION  OF  THE  FREQUENCY  DEPENDENCE  OF  COMPLEX  STIFFNESSES. 
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FRACTURE  AND  FATIGUE  OF  Bl-MATERlALS 


J.  W.  Mar 

Massachusetts  1  nstitute  of  Technology 


•  FRACTURE 

•  FATIGUE 

•  MECHANISMS 

•  CORRELATION  WITH 
FRACTURE 
MECHANICS 
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-  STR£SS  CONCBhJT/iATIOhl 


BIAXIAL  TESTING  OF  GRAPHITE/EPOXY 
COMPOSITES  CONTAINING  STRESS  CONCENTRATIONS 


BY 


I.M.  DANIEL 

I IT  RESEARCH  INSTITUTE 
CHICAGO^  ILLINOIS  60616 


OBJECTIVE 


DETERMINE  EFFECT  OF  BIAXIAL  LOADING  ON  DEFORMATION 
AND  FAILURE  OF  FLAT  GRAPHITE/EPOXY  LAMINATES  WITH  CIRCULAR 
HOLES  AND  THROUGH  THE  THICKNESS  CRACKS. 

INVESTIGATE  ALSO  NOTCH  SIZE  EFFECT. 
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SCOPE 


MATERIAL:  SP-286T300  (3M  COMPANY) 

LAMINATE  ORIENTATIONS:  [0/+45/9013,  [02/±45]g 

TESTS 


UNNOTCHED  UNIAXIAL 
NOTCHED  UNIAXIAL 
UNNOTCHED  BIAXIAL 


NOTCHED  BIAXIAL 

NOTCHED  UNIAXIAL  WITH  INCLINED  CRACKS 
(EFFECT  OF  NORMAL  STRESS  IN  DIRECTION  OF  CRACK) 


PROCEDURE 


•  USE  STRAIN  GAGES,  PHOTOELASTIC  COATINGS  AND 
MOIRE  TO  MASURE  STRAINS  IN  VICINITY  OF 
DISCONTINUITY  AND  IN  FAR  FIELD 


•  DETERMINE  STRAIN  CONCENTRATIONS,  STRAIN 
DISTRIBUTIONS,  FAILURE  MODES  AND  STRENGTH 
REDUCTION  RATIOS 


•  CORRELATE  EXPERIMENTAL  RESULTS  WITH  EXISTING 
THEORY 
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Determine : 


0  Moduli:  ^22*  ^12 

•  Poisson* s  ratios:  ^21 

•  Strengths:  ^22T, 

^lic*  ^22C’  ®12 


u  u 

o  Ultimate  strains:  ^xiT*  ^22T 


u  u  .u 

^llC’  ^^220’  12 


PROPERTIES  OF  UNIDIRECTIONAL  GRAPHITE/EPOXY 
SP-286T300 


Property 

Value 

Ply  Thickness 

0.130  mm  (0.0051  in.) 

Longitudinal  Modulus, 

151  GPa  (21.9  X  10®  psi; 

Transverse  Modulus ,  E22 

10.6  GPa  (1.53  x  10®  psi) 

Shear  Modulus ,  G^2 

6.6  GPa  (0.96  x  10®  psi) 

Major  Poisson* s  Ratio,  v^2 

0.31 

Minor  Poisson* s  Ratio,  V2^ 

0.014 

Longitudinal  Tensile  Strength, 

1401  MPa  (203  ksi) 

Ultimate  Longitudinal  Tensile 

Strain,  e^ix 

0.00925 

Longitudinal  Compressive 

Strength, 

1132  MPa  (164  ksi) 

Ultimate  Longitudinal  Compressive 

Strain,  e^lC 

0.01040 

Transverse  Tensile  Strength, 

^22T 

54  MPa  (7.8  ksi) 

Ultimate  Transverse  Tensile 

Strain,  e^2T 

0.00541 

Transverse  Compressive  Strength, 

®22C 

211  MPa  (30.6  ksi) 

Ultimate  Transverse  Compressive 

Strain, 

0.02565 

In-plane  Shear  Strength,  S^2 

72  MPa  (10.5  ksi) 

Ultimate  Shear  Strain,  e^2 

0.0110 
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UNNOTCHED  UNIAXIAL  TESTS 


SPECIMENS:  STANDARD  2.54  cm  x  23  cm  (1  in  x  9  in) 

COUPONS  OF  [0/+45/90]g  LAMINATE 


TESTS:  TWO  TESTS  EACH  IN  O-DEG.  AND  90-DEG. 

DIRECTIONS 


RESULTS : 


Layup 

Modulus 

^xx 

GPa  (lO^psi) 

Poisson’s 

Ratio 

V 

xy 

Strength 

^xxT 

MPa  (ksi) 

Ultimate 

Strain 

^xxT 

(10^£) 

[0/+45/90]g 

57 

(8.3) 

0.31 

502 

(73) 

10.0 

[90/+45/0]g 

55 

(7.9) 

0.30 

541 

(78) 

10.3 

STATIC  TENSILE  PROPERTIES  OF  UNNOTCHED 
[02/+45]g  GRAPHITE/EPOXY  LAMINATES 


Direction 

Loading 

Modulus 

*^xx 

GPa  (10^  psi) 

Poisson's 
Ratio,  v^y 

Strength 

^xxT 

MPa  (ksi) 

Ultimate 

Strain 

0-deg 

82 

(11.8) 

n 

805 

(116) 

10.0 

22.5-deg 

73 

(10.6) 

290 

(  42) 

4.0 

45 -deg 

49 

(  7.1) 

0.11 

356 

(  52) 

>9.0 

90-deg 

24 

(  3.46) 

0.20 

139 

(  20) 

7.5 
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Strength  Reduction 


UNIAXIAL  [0/t45/90]  LAMINATES  WITH  CIRCULAR  HOLES 


ISOCHROMATIC  FRINGE  PATTERNS  IN  PIIOTOELASTIC  COATING  AROUND  CRACK  OF 
SPECIl^N  NO.  4-6  AT  APPLIED  STRESS  LEVELS:  (a)  121  MPa  (17.6  ksl) , 

(b)  152  MPa  (22.0  ksi)  ,  (c)  172  IlPa  (24.9  ksi)  ,  (d)  192  MPa  (27.9  ksi) 
(e)  212  MPa  (30.8  ksi)  and  (f)  223  MPa  (32.3  ksi). 
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UNIAXIAL  [0/+45/90]g  LAMINATES  WITH  CRACKS 


cm 


2.54  (1.00) 

2.54  (1.00) 

1.91  (0.75) 

1.91  (0.75) 

1.91  (0.75) 

1.27  (0.50) 

1.27  (0.50) 

0.64  (0.25) 

0.64  (0.25) 


Modulus 
^xx  ^ 

GPa  (10°  psi) 

52 

55 

53 

(7.7) 

55 

(7.9) 

57 

(8.3) 

55 

(7.9) 

58 

(8.3) 

56 

(8.1) 

Strength 

Sxx 

MPa  (ksi) 

202 

(29) 

233 

(34) 

180 

(26) 

231 

(33) 

224 

(32) 

293 

(42) 

261 

(38) 

320 

(46) 

304 

(44) 

Experimental 

—  Point  Stress  Criterion  (d^  *  1  mm) 

—  Average  Stress  Criterion  (a^  =  5  mm) 


(in) 

H 

15  (mm) 


Half  Crack  Length,  a 
STRENGTH  REDUCTION  AS  A  FUNCTION  OF  CRACK  LENGTH  FOR  [0/+45/90]g  GRAPHITE/ 
EPOXY  PLATES  WITH  HORIZONTAL  CRACKS  UNDER  UNIAXIAL  TENSILE  LOADING 


Strength  Reduction, 


EPOXY  PLATES  WITH  CIRCULAR  HOLES  AND  CRACKS  UNDER  UNIAXIAL  TENSILE  LOADING 


CONCLUSIONS 

UNIAXIAL  TESTS  >  STRAINS 

1.  IN  THE  ELASTIC  RANGE,  THE  MEASURED  STRESS  CONCENTRATION 
IS  CLOSE  TO  THE  THEORETICAL  VALUE  OF  3  FOR  THE  UNIAXIAL 
PLATES  WITH  CIRCULAR  HOLES. 

2.  STRAINS  NEAR  THE  DISCONTINUITY  BECOME  NONLINEAR  AT  A 
STRAIN  LEVEL  OF  APPROXIMATELY  0.006e  CORRESPONDING  TO 
FIRST-PLY  FAILURE  OF  THE  90-DEG.  PLIES.  THE  ULTIMATE 
TRANSVERSE  TENSILE  STRAIN  IS  0.0054e 

3.  THE  STRAINS  NEAR  THE  CRACK  TIP  SHOW  AN  ADDITIONAL  POINT 

OF  RATE  CHANGE  OCCURRING  AT  A  STRAIN  LEVEL  OF  APPROXIMATELY 

0.002e. 

4.  MAXIMUM  STRAINS  AT  FAILURE  ON  THE  HOLE  BOUNDARY  OR  NEAR 
THE  CRACK  TIP  HAVE  EXCEEDED  TWICE  THE  ULTIMATE  STRAIN  OF 
THE  UNNOTCHED  MATERIAL. 
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CONCLUSIONS 


UNIAXIAL  TESTS  -  STRENGTH 


1.  THE  STRENGTH  REDUCTION  RATIO,  RATIO  OF  NOTCHED  TO 
UNNOTCHED  STRENGTH,  RANGES  BETWEEN  0.40  AND  0.64 
FOR  NOTCHES  BETWEEN  2.5  cm  and  0.6  cm  IN  SIZE. 

2.  STRENGTH  REDUCTION  RESULTS  FOR  UNIAXIAL  SPECIMENS  ARE 
IN  AGREEMENT  WITH  OTHER  EXPERIMENTAL  DATA  ON  GRAPHITE/ 
EPOXY,  BORON/EPOXY  AND  GLASS /EPOXY. 

3.  THE  EFFECT  OF  NOTCH  SIZE  CAN  BE  SATISFACTORILY  DESCRIBED 
BY  USING  THE  STRESS,  POINT  OR  AVERAGE,  NEAR  THE  BOUNDARY 
OF  THE  DISCONTINUITY  AS  A  CRITERION.  THE  AVERAGE  STRESS 
OVER  A  DISTANCE  OF  3.8  TO  5  mm  SEEMS  TO  GIVE  THE  BEST  FIT 

4.  THE  STRENGTH  REDUCTION  IS  INDEPENDENT  OF  NOTCH  GEOMETRY. 
SPECIMENS  WITH  HOLES  AND  CRACKS  OF  THE  SAME  SIZE  GAVE 
THE  SAME  STRENGTH  REDUCTION. 

5.  THERE  IS  A  CRITICAL  NOTCH  SIZE  BELOW  WHICH  FRACTURE  IS  AS 
LIKELY  TO  OCCUR  THROUGH  THE  UNNOTCHED  SECTION  AS  THROUGH 
THE  NOTCHED  SECTION. 

CONCLUSIONS 

UNIAXIAL  TESTS  -  FAILURE  MODES 

1.  IN  THE  SPECIMENS  WITH  CIRCULAR  HOLES,  FAILURE  INITIATES 
OFF  THE  HORIZONTAL  AXIS  AT  POINTS  WHERE  THE  INTERLAMINAR 
AND  MEMBRANE  SHEAR  STRESSES  REACH  MAXIMUM  VALUES. 

2.  NO  CRACK  EXTENSION  WAS  VISIBLE.  FAILURE  IN  THE  FORM  OF 
CRACKING  AND  DELAMINATION  INITIATED  OFF  THE  VERY  TIP  OF 
THE  CRACK  AND  PROPAGATED  AT  AN  ANGLE  TO  THE  CRACK. 

3.  FINAL  FAILURE  OCCURS  WHEN  THE  DAMAGED  AREA  AROUND  THE  TIP 
OF  THE  CRACK  REACHES  SOME  CRITICAL  SIZE. 
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406  (16) 
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SKETCH  OF  BIAXIAL  COMPOSITE  SPECIMEN 
(Dimensions  are  in  millimeters  and  inches) 


.<A.,  -i  >'-ri7T*'“«L^»'  ,MFi'*t«(5iii.>-S‘:'Pii^  ,  £-As 

«ri8iiit:«lfi 

■/i  ■:  :\*Ti‘^S:fr~'i:-C^M.^ri 

>Ms'S??S;'l;#''Sf^5-^i§ 


BIAXIAL  SPECIMEN  WITH  HOLE  MOWTTED  IN  LOADIIIO  FRA?® 
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TING  OF  FLAT  LAMINATES  AND  ASSOCIATED 


♦4» 


RECORDING  FRINGE  PATTERNS  IN  PIIOTOELASTIC  COATING  OF 
BIAXIAL  SPECIMEN 


STRESS  •  V •  ( ) 


50i 


(1  in.)  DIAMETER  HOLE  UNDER  EQUAL  BIAXIAL  LOADING  (Spec.  No.  6-7). 


50| 


-1300 


-H200 


100 


I  2  3  4  5  6  7  8 

STRAIN  .lO^c 

HORIZONTAL  STRAINS  IN  t0/+45/90lg  GRAPHITE/EPOXY  SPECIMEN  WITH  2.54  cm  (1  In.) 
DIAMETER  HOLE  UNDER  EQUAL  BIAXIAL  LOADING  (Spec.  No.  6-7). 


STRESS, ^VVf(MPa) 


ISOCHROllATIC  FRINGE  PATTERNS  IN  PHOTOELASTIC  COATING  OF  [0/+45/90] 
GRAPHITE /EPOXY  SPECIMEN  UNDER  EQUAL  BIAXIAL  LOADING  (Spec.  No.  6-6 
far-field  biaxial  stress  marked). 
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ISOCHROMATIC  FRINGE  PATTERNS  IN  PHOTOELASTIC  COATING  OF 
SPECIMEN  NO.  6-7  (Far-Field  Biaxial  Stress  Marked). 


STRAIN,  1 

FRINGE  ORDER  AND  CIRCUMFERENTIAL  STRAIN  AT  TWO  LOCATIONS  ON  THE  HOLE  BOUNDARY 
FOR  [0/+45/90]g  GRAPHITE/EPOXY  SPECIMEN  WITH  2.54  cm  (1  in.)  DIAMETER  HOLE 
UNDER  EQUAL  BIAXIAL  LOADING  (Spec.  No.  6-7). 


FAILURE  PATTERN  IN  [0/+45/90]s  GRAPHITE /EPOXY  SPECIMEN  WITH  1.91  cm 
(0.75  in.)  DIAMETER  HOLE  UNDER  EQUAL  BIAXIAL  LOADING  (Spec.  No.  6-6) 


FAILURE  PATTERN  IN  [07+45/90]^  GRAPHITE/EPOXY  SPECIMEN  WITH 
2.54  cm  (1  in.)  DIA>IETER  HOLE  UNDER  EQUAL  BIAXIAL  LOADING 
(Spec.  No.  6-7). 
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Strength  Reduction, 


Table  1 


BIAXIAL  [02/145/90] g  LAMINATES  WITH  CIRCULAR  HOLES 


Spec. 

No. 

Hole  Diam. 

2a 

cm  (in.) 

Failure 

Stress 

^xxT“®yyT 

MPa  (ksi) 

Maximum 
Strain  at 
Failure 

_ (lO'^e) _ 

Strength 

Reduction 

Ratio 

S  /S 
xx'  0 

9=0" 

0=22.5" 

6-7 

2.54 

(1.00) 

293 

(42) 

17.5 

0.582 

6-8 

2.54 

(1.00) 

258 

(37) 

iB 

10.7 

0.512 

6-6 

1.91 

(0.75) 

299 

(43) 

10.0 

19.5 

0.593 

6-9 

1.91 

(0.75) 

261 

(38) 

10.3 

11.8 

0.518 

6-10 

1.27 

(0.50) 

312 

(45) 

- 

0.619 

6-11 

1.27 

(0.50) 

328 

(48) 

ByQ 

14.9 

0.651 

6-15 

0.64 

(0.25) 

383 

(56) 

15.0 

18.5 

0.761 

6-22 

0.64 

(0.25) 

349 

(51) 

14.3 

17.3 

0.693 

i.oa^ 

\ 


0.8 


1-  ^  \ 


0.6 


0.4 


\ 


\ 


\ 


N; 


i- 


$ 


Biaxial  (0^^^  -  a^) 


Uniaxial 


,  0.25  0.50 

'I - ^ - - \ - 

®  ^  Hole  Radius .  a  ® 

STRENGTH  REDUCTION  AS  A  FUNCTION  OF  HOLE  RADIUS  FOR  (0/+45/90J,  GRAPHITE/ 
EPOXY  PLATES  WITH  CIRCULAR  HOLES  UNDER  UNIAXIAL  AND  1 :  l“ BIAXIAL  LOADING 


Strength  Reduction, 


CIRCULAR  HOLE  IN  BIAXIAL  STRESS  FIELD 


J  f  I  i  M  L 

Stresses : 


■“I  I  t  I  I  I  r 

Average  Stresses  (In  annulus  r  *»  a  to  r  =*  r^^) 


®rr  ' 

»ee  '  % 


a 


Failure  Criterion 

•  Failure  occurs  when  a  ^  and  fall  on  failure  envelope 

rr  V  w 

for  laminate 

•  (r^^  -  a)  =  characteristic  dimension 
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CONCLUSIONS 


[0/+45/90]  LAMINATES  WITH  HOLES  UNDER  EQUAL  BIAXIAL 

'*^ _ S  - 

TENSION 


STRAINS  AND  FAILURE 


1.  STRAINS  ON  THE  BOUNDARY  OF  THE  HOLE  BECOME  NONLINEAR 

AT  APPLIED  STRESSES  OF  »  120  MPa  (17.5  ksi) 

CORRESPONDING  TO  A  STRAIN  LEVEL  OF  APPROXIMATELY  4  x 
lO’^e.  THIS  IS  THE  STRAIN  LEVEL  AT  WHICH  THE  RESPONSE 
OF  THE  UNIAXIAL  UNNOTCHED  SPECIMENS  BECOMES  NONLINEAR. 

2.  INITIALLY,  THE  CIRCUMFERENTIAL  STRAIN  IS  UNIFORM  AROUND 
THE  BOUNDARY  OF  THE  HOLE.  SUBSEQUENTLY,  WITH  INCREASING 
LOAD,  REGIONS  OF  HIGH  STRAIN  CONCENTRATION  WITH  NON¬ 
LINEAR  RESPONSE  DEVELOP  AT  EIGHT  CHARACTERISTIC  LOCATIONS 
22.5-DEG.  OFF  THE  FIBER  AXES. 

3.  MAXIMUM  STRAINS  AT  FAILURE  ON  THE  HOLE  BOUNDARY  REACH 
VALUES  UP  TO  19.5  x  10'^  COMPARED  TO  THE  ULTIMATE 
STRAIN  10"^  OF  THE  UNNOTCHED  LAMINATE. 

4.  FAILURE  INITIATES  ON  THE  HOLE  BOUNDARY  AT  POINTS  22.5 
DEG.  OFF  THE  HORIZONTAL  AND  VERTICAL  AXES. 
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CONCLUSIONS 


[0/+45/90]g  LAMINATES  WITH  HOLES  UNDER  EQUAL  BIAXIAL 
TENSION 


STRENGTH 

1.  THE  STRENGTH  REDUCTION  RATIO  OF  NOTCHED  TO  UNNOTCHED 
STRENGTH  RANGES  BETWEEN  0.55  AND  0.73  FOR  HOLES 
BETWEEN  .^.54  cm  and  0.64  cm  IN  DIAMETER. 

2.  THE  STRENGTH  REDUCTION  RATIOS  ARE  HIGHER  THAN 
CORRESPONDING  VALUES  FOR  UNIAXIAL  LOADING  BY  26  TO 
33  PERCENT. 

3.  THE  VARIATION  OF  STRENGTH  REDUCTION  RATIO  WITH  HOLE 
DIAMETER  CAN  BE  SATISFACTORILY  DESCRIBED  BY  USING  AN 
AVERAGE  BIAXIAL  STRESS  CRITERION  OVER  A  RADIAL 
DISTANCE  OF  3  mm. 

4.  THE  STRENGTH  REDUCTION  RATIOS  FOR  UNIAXIAL  AND  EQUAL 
BIAXIAL  LOADING  REPRESENT  LOWER  AND  UPPER  BOUNDS.  THE 
STRENGTH  REDUCTION  RATIO  FOR  ANY  OTHER  BIAXIALITY 
RATIO  WOULD  FALL  BETWEEN  THESE  TWO  EXTREMES. 


DESIGN  RECOMMENDATION 


QDASI-ISOTROPIC  COMPOSITE  PANELS  WITH  CIRCULAR  CUTOUTS 
CAN  BE  DESIGNED  CONSERVATIVELY  BY  ASSUMING  STRENGTH 
REDUCTION  RATIOS  CORRESPONDING  TO  UNIAXIAL  LOADING. 
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OBJECTIVE : 


SCOPE; 


RESPONSE  OF  GRAPHITE  COMPOSITES 
TO  LASER  RADIATION 
(F44620-76-C-0030) 


K.  G.  KIBLER 

Materials  Research  Laboratory 
GENERAL  DYNAMICS 
Ft.  Worth,  TX  76101 


DEVELOP  A  BASIC  UNDERSTANDING  OF  THE  PHYSICAL  BEHAVIOR 
OF  GRAPHITE-EPOXY  COMPOSITE  MATERIALS  EXPOSED  TO 
CW  CO2  LASER  RADIATION 

1.  DETERMINE  PENETRATION  RATE,  THERMAL  RESPONSE, 

AND  STRENGTH  DEGRADATION  OF  COMPOSITES 

2.  IDENTIFY  THERMOPHYSICAL  MECHANISMS  OF 
DEGRADATION,  OTHER  THAN  COMPLETE  PENETRATION 

3.  FORMULATE  SIMPLE  MODELS  TO  PREDICT  PENETRATION 
RATE  AND  STRENGTH  DEGRADATION  IN  COMPOSITES 

4.  COMPARE  EFFECTS  ON  COMPOSITES  WITH  EFFECTS 
ON  ALUMINUM 
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BACKGROUND:  J.  COMPOSITE  MATERIALS  9,  28  (1975) 


CWARACTERIZATIOK)  Of 

Strew^th  De^radat/ow 


NbE 


MEASltRE’ 

THERMAL 

Response 


INDUCE 

Damage 
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MEASaftE 

STREWQTH 


CoAReLATE'  i 
PREDICT 


TEST  SPECIMEN)  ScchhARV 


Effects 


Use«.  resr  CoK>recu 

SpeciHcMS  speciHEMS 

-(  0/i:^s')s - U.,  4'8"^t'>/A“ct'> 

-fo/±45r)« - u.  vMX(T):i''a\'ov 

z  (o/tAs)^f - U-,  ''q"Ct),  'A^IT) 

LdWii.vs.  T^Awsv. 

St«cio<;th 

I 

■-Cio/±4s)zs - U, 

-(±As/‘}o)2s - LL 

asA - u,  V‘a),  '4"(T),  K(V 


LAse« 

ExpePinewTi 

DME’ 

P9W£R 

Sadw  dePrn 

THeRhac  Response’ 

PRCLoAO 
Moist  aer 
(spot  Siir  ") 


MflSS  Loss  vs.  \NC.10eMT 


INCIDEMT  ENERGY 


Am  (e) 

5208/T 300 


LASER  DAMAGE  IN  5208 /T300  (0/+45)2s 

INCIDENT  BEAM  750  w  /  0.5  cm^ 


2 


25,6  sec. 
(Burnth rough) 


5 


12.8  sec. 


LASER  DAMAGE  IN  5208/T300  (0/+45)2g 

INCIDENT  BEAM  750  w  /  0,5  ctn^ 


X-RAY  WITH  TBE 


5  lo  15. 

BEAM  EMERQY 


LASER-lMtiU.Ce’D  STREM<,th 
ReoutTiOM  IN  (o/±4S)^^  5‘2£iS/TJoo 


A^ho\  • 


D 


o  a 


zo  e:  (o) 


2 


57 


30 


oA  Am  (□) 

40  L%Os  (•) 


PRINCIPAL  RESULTS  TO  DATE  (  1  SEPT  1976  ) 


O  REQUIRED  BURNTHROUGH  ENERGY  FOR  ALUMINUM  ABOUT  1/3 
THAT  FOR  GRAPHITE-EPOXY 

O  ROUGH  CALCULATION  FOR  EFFECTIVE  HEAT  OF  ABLATION  GIVES 
VALUE  UNREALISTICALLY  HIGH  -  -  SIGNIFICANT  CONDUCTION 
LOSS  IN  G/E  ? 

o  CONSIDERABLE  DATA  SCATTER  FOR  A1  --  NOT  SO  FOR  G/E 

o  NO  EFFECT  OF  STACKING  SEQUENCE  ON  MASS  LOSS  OR 
PENETRATION  DEPTH  -  - 
SLIGHT  THICKNESS  EFFECT 

O  STRENGTH  REDUCTION  IN  G/E  APPEARS  PREDICTABLE 
USING  "SIMPLE"  PARAMETERS 
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DEFECT  -  PROPERTY  RELATIONSHIPS  IN  CONPOSITE  I^TERIALS 


Investigators 

E.  6.  HENNEKE  IL  K.  L.  REIFSNIDER.  W.  W.  STINCHCOMB 
Virginia  Polytechnic  Institute  S  State  University 
BlacksburGv  Virginia 


OBJECTIVES; 

1.  To  IDENTIFY  THE  PRECISE  NATURE  OF  DAMAGE  DEVELOPMENT  IN  QUASI -ISOTROP IC 
GRAPHITE-EPOXY  LAMINATES  UNDER  VARIOUS  LOAD  HISTORIES. 

2.  To  DETERMINE  THE  PHYSICAL  PARAMETERS  WHICH  LEAD  TO  A  LOSS  OF  STRENGTH 

and/or  life. 

3.  To  ESTABLISH  THE  MECHANICS  OF  THE  INDIVIDUAL  AND  COMBINED  ACTION  OF  THESE 
PARAMETERS  AS  THEY  INFLUENCE  MECHANICAL  RESPONSE. 

M.  To  ADDRESS  THE  QUESTION  OF  HOW  THESE  FINDINGS  CAN  BEST  BE  DESCRIBED  BY 
ANALYSIS. 


Figure  1 


INVESTIGATIVE  PROGRAM: 

Type  I  material  -  [O.iAB.OOlg.  AS3501  Graphite-Epoxy 
Instrumented  tensile  tests 

Laminate  theory  calculations  -  failure  predictions 
SEM  and  light  microscope  studies 
Acoustic  emission  and  ultrasonic  attenuation  studies 
Fatigue  loading 

Video  and  thermographic  studies 
Sectioning  studies 

Type  II  materials  -  [0j90^±A5]g^  AS3501  Graphite-Epoxy 
Instrumented  tensile  tests 

Laminate  analysis  and  3-D  FEii  -  failure  predictions 

SEM  AND  LIGHT  MICROSCOPE  STUDIES 

Acoustic  emission  and  ultrasonic  attenuation  studies 
Instrumented  fatigue  loading 
VlDEO^  CINEj  and  thermographic  STUDIES 
Sectioning  studies 

General 

Load  history  studies 

Transverse  crack  and  delamination  investigation  -  initiation,  growth  and  fracture 
Technique  development  -  ultrasonic  attenuation,  thermography,  replication 

Analysis  evaluation  and  development 
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TABLE  III 


SUMMARY  OF  PROPERTIES 
FOR  AS-3501 

Elastic 

Stiffness 

Tensile 

Properties 

Value 

El 

20.2x10®  psi 

E2 

1.4x10®  psi 

'’12 

0.28 

®12 

0.65x10®  psi 

Strength  Properties 

X 

235,000  psi 

Y 

8,200  psi 

S 

17,900  psi 

Other  Properties 

Thermal  expansion  coefficients; 

a.j  (fiber  direction)  =  -0.2x10®  in/in/^F 
og  (transverse  direction)  =  13.0x10  ®  in/in/°F 
Fiber  volume  fraction  0.62 
Ply  thickness  0.0052  +  0.0004  in. 

Void  content  Z% 

Density  0.057  Ib/in^ 

Stress-free  temperature  278°F 

*  Data  supplied  by  Hercules,  Inc. 

Figure  3 


Compression 

Value 

17x10®  psi 

1.6x10®  psi 

0.28 

0.65x10®  psi 

180,000  psi 
25,000  psi 
17,900  psi 
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5TRE:55  [  K5Z  ] 


Figi  ^  “  Transverse  stress  distribution  through  the  thickness  fop  Type  I  material 


Z  [ PLY  thickness] 

Fi6.  5  -  Interlaminar  normal  stress  distribution  through  the  thickness  for  Type  1  material 
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APPENDIX 


For  a  stress  free  temperature  which  differs  from  room  temperature  the  midplane  strains 
ARE  given  by 


[e°]  =[a]"^N  +  (a) 

WHERE  [e]  IS  THE  MATRIX  OF  MIDPLAME  STRAINS^  [a]”^  IS  THE  INVERSE  OF  THE  PLATE  STIFFNESS 
MATRIX  IN  THE  ABSENCE  OF  BENDING^  [N]  IS  THE  COLUMN  MATRIX  OF  APPLIED  IN-PLANE  RESULTANTS^ 
AND  [n"'']  IS  THE  AMBIENT  TEMPERATURE  RESULTANTS  DUE  TO  THERMAL  RESIDUAL  STRESSES  GIVEN  BV 

[N^]  =  JtQ][a]ATdz  (b) 

WHERE  t“]  ARE  THE  THERMAL  EXPANSION  COEFFICIENTS^  LISTED  IN  TABLE  III,  ThE  STRESSES  IN 
A  GIVEN  LAYER  ARE  THEN 


[a]'^  =[Q]'^[[^'’l  -  l«]AT]'^ 
The  PERTINENT  MAXIMUM  STRAIN  CONDITION  IS^ 


- 2 - 2“ 

SIN  e-''2TC0S  e 


(C) 


(d) 


Fig.  6  •  Summary  of  maximum  strain  criterian  including  thermal  stresses 
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Delamination  and  transverse  cracking  in  Type  I  material  at  2700  lb. 


Fig,  9  -  Typical  observed  stress-strain  curves  for  Type  I  material 


VIDEO-THERflOGRAPHY  STUDY 

1.  Cycled  type  one  specimen  at  30  HZ  at  82%  of  static  failure  load.  Specimen 

HAD  1/4  IN.  CENTER  HOLE.  BrOKE  AFTER  A  FEW  HUNDRED  CYCLES.  VeRY  LITTLE 
HEAT  GENERATION  OR  CHANGE. 

2,  Similar  to  1.  At  70%  of  failure  load.  No  substantial  pattern  change. 

3  Ran  some  BEp  and  BAl  for  comparison.  BAl  is  similar  in  heat  output. 

4.  Introduced  variable  frequency  excitation  -  "Vibro-thermography",  and 
identified  damage  more  readily  IN: 

i)  0°  GEp  after  static  loading 

II)  NOTCHED  UNIDIRECTIONAL  GEf  AFTER  STATIC  LOADING 
III)  QUASI-ISOTROPIC  LAMIANTE  WITH  HOLE 

iv)  ±60°  GEp 

v)  90°  GEp  unidirectional 

VI )  TYPE  I  AND  II  GEp  without  HOLE 

5.  Identified  grip  region  as  a  special  problem. 

Fig.  19  -  Summary  of  video-thermography  study 
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(0,^45,90)3 


0  -45  90  ?45  0 


Fig.  IJ  -  Schematic  diagram  of  damage  development  in  Type  I  material 


TYPE  I 

Transverse  cracking  of  90°  plies: 

1.  Commonly  observed  at  1500-1800  lb  under  load. 

2.  Cracks  propagate  through  thickness^  first  through  90°  plies^ 
then  into  ±i|5°  plies. 

3.  Branching  in  ^5°  plies  is  common^  usually  at  2500-2700  lb, 

^1.  Delamination  occurs  beginning  at  transverse  cracks  and  relaxes  the 

BRANCHES  OF  THE  PRIMARY  CRACK. 


Fig,  12  -  Summary  of  damage  process  in  Type  I  material 
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5TRE55  [  K5I  ] 


Fie.  13  -  TRAHSVEftSE  STRESS  DISTRIBUTION  THROUGH  THE  THICKNESS  FOR  TyPE  II  MATERIAL 


12  3  4- 

Z  [  PUY  THICKNESS  ] 


Fig.  W  -  Interlaminar  normal  stress  distribution  through  the  thickness  for  Type  II  materiau 
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Fig.  17  Edge  replicas  of  Type  II  speci.'iens  at  200  lb.  and  2700  lb. 


FigurelS-  Ultrasonic  transducer  clumped  to  a  composite  specimen 
mounted  in  testing  grips. 
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5TRHIN  C 10 


Fig.  19  -  Real-time  attenuation  and  strain  verses  load  for  a  test  of  Type  II  material 


OSCILLOSCOPE 


PULSED  ATTENUATION 


Fig.  20  -  Schematic  of  real-time  attenuation  apparatus 
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DB/EII-HP 


KIPB 


STRAIN  GAGE  STUDY: 

1.  Observed  as  much  as  111  difference  in  readings  from  edge  to  centerline. 

2.  Differences  can  be  caused  by  variations  in  tab  thicknesseSj  imprecise 
alignment  of  grip  faces^  or  nonuniform  grip  purchase. 

3.  Fracture  initiated  from  high  strain  region. 

4.  Obtained  data  wherein  fracture  occurred  through  the  gages  -  strains  were 

ESSENTIALLY  IDENTICAL  TO  FRACTURE  “  NO  EVIDENCE  OF  PREMATURE  RELAXATION 
AT  EDGES  OR  CENTER. 

5.  Temperature  was  found  to  influence  strain  gage  readings. 


Fig.  23  -  Summary  of  results  of  strain  gage  study 


CENTERLINE  SECTIONING  SERIES: 

1.  Type  I  specimen  loaded  to  1500  lb 

a.  essentially  same  number  of  major  transverse  cracks  on  interior 
&  exterior  edges  in  90°  layers 

B.  more  cracks  in  45°  layers  in  interior 

C.  MANY  more  angular  AND  MINOP  CRACKS  IN  INTERIOR  -  LARGEST  SINGLE 
CRACK  FOUND  IN  INTERIOR 

D.  90°  CRACKS  MORE  VISIBLE  AT  lOOX  ON  LIGHT  MICROSCOPE  UNDER  LOAD  THAN 
AT  400X  ON  SEM^  UNLOADED 

E.  CRACKS  IN  45°  LAYERS  MORE  VISIBLE  IN  SEM  (UNLOADEDI  THAN  UNDER  LOAD 
WITH  LIGHT  MICROSCOPE 

F.  No  CRACKS  ON  OUTSIDE  EDGE  OBSERVED  ON  SE'I 

2.  Type  I  specimen  loaded  to  1000  lb 

A.  TRANSVERSE  CRACKS  WERE  FOUND  IN  INTERIOR  WITH  SEK 

B.  NO  CRACKS  WERE  FOUND  ON  EXTERIOR  EDGE  IWTH  SEM 

C.  CRACKS  IN  45°  LAYERS  IN  INTERIOR  ARE  FREQUENTLY  ANGULAR 

D.  TRANSDUCER  REGION  DOES  NOT  SHOW  ANY  IRREGULARITIES 

E.  AT-LOAD  LIGHT  MICROSCOPE  STUDY  SHOWED  LARGER  NUMBER  OF  CRACKS  ON 
INTERIOR  THAN  EXTERIOR  EDGE.  RESULTS  ON  TWO  HALVES  SHOWED  GOOD 
REPRODUCIBILITY 

Fig.  24  -  Summary  of  results  of  centerline  sectioning  series 
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3.  Type  II  specimen  loaded  to  1000  lb  stress  amplitude  for  850^000  cycles 

A.  LIGHT  MICROSCOPE  EXAMINATION  AT-LOAD  SHOWED  LARGER  NUMBER  OF 
TRANSVERSE  CRACKS  THAN  FOR  THE  SPECIMEN  LOADED  STATICALLY 

B.  MANY  MORE  TRANSVERSE  CRACKS  WERE  FOUND  IN  THE  INTERIOR  EDGE  THAN 
ON  THE  EXTERIOR  EDGE 

C.  MORE  INTERFACIAL  AND  IRREGULAR  CRACK  FORMATION  BEHAVIOR  VIAS  NOTED 
THAN  FOR  THE  STATIC  LOADING  CASES 

D.  THE  LARGEST  CRACKS,  ESPECIALLY  IN  THE  ANGLE  PLIES,  VIERE  FOUND  IN 
THE  INTERIOR  EDGE 

4.  Type  I  specimen  loaded  at  1000  lb  load  amplitude  for  850  thousand  cycles 

A.  TRANSVERSE  CRACK  COUNT  WAS  36  ON  ONE  EXTERIOR  EDGE  (WHICH  DELAMINATED) 
AND  20  ON  THE  OTHER,  INTERIOR  EDGE  COUNT  WAS  24. 

B.  DELAMINATION  OCCURRED  AT  ONE  END  AND  NEAR  THE  CENTER  ON  ONE  EDGE  ONLY 

C.  STIFFNESS  APPEARS  TO  DECREASE  WITH  DAMAGE  BUT  INCREASE  WITH  CYCLES  WHEN 

DAMAGE  IS  NOT  DEVELOPING 

5.  Type  I  specimen  in  as-received  condition 

A.  observed  NINE  OR  TEN  INITIAL  CRACKS  WHICH  HAVE  THE  APPEARANCE  OF  SMALLER 
versions  of  WHAT  WE  FIND  IN  LOADED  SPECIMENS 

B.  Interior  and  exterior  edges  appear  to  be  essentially  identical 

C.  THE  INITIAL  CRACKS  THAT  ARE  SEEN  DO  NOT  ALWAYS  APPEAR  WITH  ORIENTATIONS 
transverse  to  the  specimen  LONGITUDINAL  AXIS 

D.  A  PREDOMINANT  PROPORTION  OF  THOSE  INITIAL  CRACKS  WHICH  WERE  OBSERVED 

WERE  FOUND  TO  BE  LOCATED  IN  A  REGION  WITHIN  1  CM  OR  SO  FROM  THE  TABS 
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0^  +  -  t|5°  90°  -  +  i|5°  0°  0°  +  -  i|5°  90°  -  +  45°  0° 


INTERIOR  EDGE 

Fig.  25  -  Schematic  diagram  of  crack  branching  on  interior  and  exterior  surfaces 


FIG.  26-through-the-thickness  damage  sequence  in  type  i 

SPECIMENS  (edge  VIEW) . 


TYPE  II 


ESI 


FAILURE  LOAD  (lBS) 


Fig,  27  -  Distribution  of  failure  loads  for  Type  I  and  Type  II  specimens 


TYPE  II 

1.  "knee"  in  stress-strain  curve  nearly  imperceptible  -  MORE  like  a  gentle 
curve. 

2.  Transverse  cracking  more  difficult  to  see  -  first  observe  at  1800-2300  lb. 

3.  The  growth  of  transverse  cracks  is  slower  than  in  type  i.  Spread  into 
45°  PLIES  occurs  between  3000-3200  lb. 

4.  Delamination  does  not  stop  the  transverse  propagation^  but  a  complete 
transverse  crack  through-the-thickness  forms  only  in  a  few  instances. 

5.  Formation  of  a  through-the-thickness  transverse  crack  occurs  quite  close 

TO  THE  FINAL  FAILURE  LOAD.  USUALLY  IMMEDIATELY  PRECEEDING  IT. 


Fig.  28  -  Summary  of  damage  process  in  Type  II  material 


DAMAGE  CHARACTERIZATION: 

1.  Ultrasonic  attenuation  and  acoustic  emission  indicate  damage  development 
BEFORE  visual  CRACK  OBSERVATION  OF  "kNEE"  IN  CURVE>  AT  LEVEL  QUITE 
CLOSE  TO  PREDICTED  ONE. 

2.  Cracks  do  not  appear  on  both  edges  at  the  same  time  =>  in  plane  growth, 

3.  No  evidence  that  cracks  grow  from  edge  to  center. 

4.  Need  to  understand  growth  process  to  characterize  damage  correctly^  and 

THEREBY  ESTABLISH  CONTROLLING  PARAMETERS  FOR  STRENGTH  AND  LIFE. 

Fig.  29  -  General  damage  characterization 


PRINCIPAL  FINDINGS  AND  CONCLUSIONS 

1.  In  quasi-isotropic  laminates  it  would  appear  that  "first-ply  failure"  is  more 
correctly  interpreted  as  the  initiation  of  a  process  than  the  discrete  occurrence  of  a 

SINGLE  EVENT. 

2.  If  thermal  residual  stresses  are  included  in  the  calculations^  laminate  theory 

PREDICTIONS  ARE  WITHIN  ABOUT  5%  OF  THE  LEVEL  AT  WHICH  OUR  NDT&E  TECHNIQUES  INDICATE 
DAMAGE  INITIATION. 

3.  The  "knee"  in  the  stress-strain  curve  is  more  or  less  well  defined  depending  upon 

THE  CONSTRAINT  IN  THE  REGION  OF  CRACK  FORMATION  IN  THE  90°  PLIES  AND  IS^  THEREFORE^  A  FUNC¬ 
TION  OF  STACKING  SEQUENCE.  OUR  OBSERVATIONS  OF  THE  KNEE  GENERALLY  OCCUR  AT  HIGHER  STRESS 
LEVELS  THAN  DO  THE  FIRST  INDICATIONS  AND  OBSERVATIONS  OF  MICRO-DAMAGE^  SOMETIMES  AS  MUCH 
AS  50-75%  HIGHER. 

4.  Branching  of  transverse  cracks  from  the  90°  plies  into  the  adjacent  plies, 

GENERAL  CRACKING  OF  THE  ANGLE  PLIES  APPEARS  TO  BE  MORE  EXTENSIVE  IN  THE  INTERIOR  OF  THE 
SPECIMEN  THAN  AT  THE  EDGE.  TRANSVERSE  CRACKS  COMPLETELY  THROUGH  THE  THICKNESS  APPEAR  TO 
FORM  AT  THE  EDGE. 

5.  The  strength  of  the  [0°j90°^±45°]g  laminates  is  somewhat  greater  than  the 
[0°^±45°,90°]5  laminates. 

6.  These  materials  do  exhibit  memory  effects  in  some  cases,  time  delayed  fracturf  at 
constant  elongation  has  been  observed. 

7.  Ultrasonic  attenuation  change  appears  to  be  a  clear,  precise,  consistent  and 

REPRODUCIBLE  INDICATION  OF  DAMAGE  INITIATION  AND  GROWTH.  It  CORRELATES  WELL  WITH  MICRO¬ 
SCOPE  observations. 


Figure  30 
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FATIGUE  IN  ADHESIVELY  BONDED  LAMINATES 


J.  A.  ALIC 

DEPARTMENT  OF  MECHANICAL  ENGINEERING 
WICHITA  STATE  UNIVERSITY 

OBJECTIVES 

CRACK  INITIATION  {7075-T6) 

ISOLATION  OF  FATIGUE  DAMAGE 

S-N  PERFORMANCE,  WITH  AND 
WITHOUT  STRESS  CONCENTRATION 

CRACK  PROPAGATION 

EFFECT  OF  LAMINA  THICKNESS 
(7075-T6) 

7075-T6(51)  BONDED  TO  2024-T3(51) 


LAMINA  TYPF.cj 
THIN  LAYERS 

8  LAYERS.  0.032  INCH  THICK 
22  LAYERS,  0.011  INCH  THICK 
THICK  LAYERS 

2  LAVERS.  ’A  INCH  THICK 
4  LAYERS.  Vs  INCH  THICK 
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CRACK  PROPAGATION 
COMPACT  SPECIMEN.  R=0.1 

MONOLITHIC  7075-T6(51)  AND 
2024-T3(51) 

THIN  LAMINAE  {7075-T6) 

8  LAYER 
22  LAYER 

THICK  LAMINAE  (7075-T651) 

2  LAYER 
4  LAYER 

COMBINED  7075/2024 
2  LAYER 
4  LAYER 


CRACK  INITIATION  AND  S-N 

CYCLIC  TENSION,  R=0.1 

UNNOTCHED  AND  NOTCHED  (CENTRAL 
HOLE.  1Tt=2.5) 

MONOLITHIC,  ’A  INCH  THICK 
8  LAYER 
22  LAYER 
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BACKGROUND 


CRACK  DIVIDER 

hECHANISM 

CRACK  ARRESTOR 

THICKNESS  EFFECT 

DISCONTINUITY 

ON  CRACK  GROWTH 

AT  INTERFACE 

REAL  STRUCTURE  ^WY  BE  INTERMEDIATE  SITUATION 


PRELIHINAfTl'  RESILTS 


THIN  LAMirWE 
;8  OR  22  LAYERS) 


ADIESIVE:  HYSOL  EA  9410 


TENSILE  PROPERTIES 


YS 

(Ksn 

TS 

(KSI) 

EIDNG. 

(Z) 

(KSl/IR) 

8-U\YER 

74 

82 

15 

-55 

22-LAYER 

72 

79 

16 

52 
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DELTA  K  (MPA  (M)  ••1/2)  DELTA  K.  (MPA  (M)  ••1/5) 


CONCLUSIONS 


SOTEWHAT  LaCR  GRa^H  RATES  IN  22-LAYER  LW1INATES  At 
INTERMEDIATE  LOAD  LEVELS. 

PROPAGATim  ERRATIC  IN  LAMINATES  AT  LOW  GROVIFH  RATES. 
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FAILURE  MECHANISMS  AND  CRITERIA  I N  COMPOSITE  MATERIALS 

N.  J.  PAGANO 

AIR  FORCE  MATERIALS  LABORATORY 
In-House  Work  Unit 
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VARIATIONAL  PRINCIPLE  FOR  LAMINATES 


where 


/  FdV  -  /  P.U.dS 

V  S' 


Traction  boundary  surface 


W  “  Complementary  energy  =  W(Ojj,ej.) 


ej|  «  expansional  strain  tensor 
P|  *=  prescribed  tractions 


For  a  body  v/ith  interfaces,  we  get 
N  f/U:  :  +  U:  : 


8q  ‘  s  / 

K'l  V,. 


I.J  J,l 


5))  Hi  •  \j«"i  ] 


.  N-I  .  /(!<)  (K)  (K+1)(K+1)\ 

+  /Ti8u.ds+  S  y  Ir  8u.  +r  Su  /di^  -  o 

s" '  '  K=i  1/  '  ’  '  '  ' 


Curing  Stresses 


Curing  Stresses  caused  by 


a)  Non-uniform  inelastic  deformation  of  resin 
(microstresses) 

b)  Thermal  deformations  of  "hard"  layers 
{ Effective  Modulus  theory  applies) 


For  the  latter,  we  assume 


crj(T)=  Q.j(T)[e;(T)-ej(T)] 


inelastic  ' 


e.  =  free  thermal  exp. 


and  e.  (T)  “  e,  (T)  -  e.  (T  )  (T<T  ) 

j  j  j  0  0 


cri(T^)  =0 


As  a  first  approx.,  we  may  pi;t  (MT)  ’  a.IT)  (T-T^).  Ihcii  only 
R.T.  mear.nreiiimits  are  r<.''|i/irMl. 


MAXIMUM  STRAIN  CRITERION 


TRANSVERSE  LAYER 


€  =  -  +  C 

f  T 


LAMINATE  STRENGTH 


'  ^22^^ 


(:k-  ^ 


TABLE  3  PREDICTED  APPLIED  STRESS  AT  FIRST 
PLY  FAILURE  (psi) 


WITH  CURING  STRESS  WITHOUT  CURING  STRESS 


LAMINATE/CRITERION 

MAX.  0- 

MAX .  e 

MAX .  cr 

MAX.  e 

[90,0.-45]^ 

20, 530 

19, 150 

41, 480 

37, 680 

(0,90]^ 

25, 605 

26, 100 

51,730 

51 , 370 

[90, ±30]^ 

20, 530 

19, 150 

41,480 

37,  680 

[90,i45]^ 

11.470 

7,890 

20, 390 

18, 530 

. 
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SWELLING  STRAINS 


=  -f  msc 

FOR  T300/5208: 

=  -2167  X  10"® 
AND  Sc  =  1.6 

THUS  6^"^’=  /ef>/ 

WHEN  m  =  .27% 


TABLE  III.-  CURING  STRAINS  AND  CURING  STRESSES 


-e°T.  % 

.e°T 

90 

.  % 

%5  * 

% 

< 

EXP. 

CAL. 

EXP. 

CAL. 

EXP. 

CAL. 

T 

^LT 

B/Ep 

0.  123 

0.109 

0.234 

0.259 

0.171 

0.184 

0.99 

0.  1  3 

B/PI 

0.  100 

0.075 

0.214 

0.186 

0.163 

0.131 

4.09 

0.33 

Gr/En 

-0.018 

(-0.009) 

-0.019 

0. 104 
(0.068) 

0. 122 

■1 

0.052 

1.09 

0.11 

Gr/PI 

0.002 

-0.004 

0.051 

0.017 

0.028 

1.20 

0.11 

Gl/Ep 

0.  105 

0.117 

0.260 

0.361 

0.  182 

0.239 

1.39 

0.  65 

Remarks:  Numbers  inside  the  parentheses  for  Gr/Ep  are 
curing  and  the  heating  stage  of  postcuring. 

measured  during  the 

:ooUng  stage  of 
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Figure  12.  Transverse  Strain  During  Thermal  Cycling 


TIME,  min. 

Figure  13.  Moisture  Desorption  During  Thermal  Cycling 
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STRAIN 
IN  90° 
% 


STRAIN 

% 


STRESS -FREE  TEMPERATURE,  °K 

Figure  14  Effect  of  Stress-Free  Temperature  on  Laminate  Thermal 
Strain  Between  Room  Temperature  and  450°K 


ANGLE,  DEG. 

Figure  15.  Curing  Strains  of  [02^±6^s  Laminates 


STRAIN 
IN  0° 
% 
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0 


20 


40  60 

ANGLE,  DEG. 


80 


100 


Figure  16.  Curing  Stresses  Within  0°  Plies  of  [02^^^  Gr/Pl  Laminates 


Figure  17.  Effect  of  Preloading  on  Total  Number  of  A.  E.  Counts 
at  Two  Different  Temperatures 


FAILURE  SURFACE  IN  (r„-(r,  _  PLANE 
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PLY  FAILURE  UNDER  COMBINED  LOADING 

•  UNI  DIRECTIONAL  TUBES 

•  Sj’j  ■  EXCEPT  I N  ONE  CASE  (0"’) 

• 

^12  ^12 
•  *16  • 

•  APPLICATION  OF  INVARIANTS 

•  STRENGTH  DATA 

•  ANGLE-PLY  TUBES  ([±0]^^) 

•  FIRST-PLY  FAILURE 

•  BEHAVIOR  ARER  FPF 

•  FINAL  FAILURE 

•  EFFECT  OF  PLY  FAILURE  ON  I NITI AL  MODULI 
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Figurell.  Laminate  Damage 


a  ^12  ^13  14  ^15  16 


/  14^15  16  \ 

^22  ^23  ^'24  ^25  ^26 
^33  *^34  ^33  ^36 

I  I  t 

symm. 

^55  ^56 


'66 


Yn  "  ®rz 


I  Yze-^zs 
\  Yre-  ^re 


For  axisymmetric  loadings 


u  *  U(rl 

V  *  V  (r)  +  Arz 

€g-U/r 

Xrz'W.r 

if  =  U.r 

Xze  “  Ar 

/r-^‘=V,r  -  V/r 

Equilibrium  Equations 

(r  cr^),  ^  -cTe  =  0 

^re  ■ 

This  yields  3  D.E.  with  variable  coefficients  in  U,V,VV. 
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Governing  Equations 
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ULTRASONIC  NON-DESTRUCTIVE  TESTING 
OF  COMPOSITE  MATERIALS 


W.  Sachse  and  Y.H.  Pao 

DEPARTMENT  OF  THEORETICAL  AND  APPLIED  MECHANICS 

Cornell  university 
ITHACA.  N.Y.  1A853 


(April  I,  1976  -  March  31,  1977) 

ACCOMPLISHMENTS  (as  of  October,  1976) 

1,  PROCURED  composite  SPECIMENS  AND  FABRICATED  TEST 
specimens  of  various  FIBER  ORIENTATIONS, 

GRAPHITE-EPOXY  (0/  +  45)  ,  BORON-EPOXY  (UNIDIRECTIONAL) 

2,  SET  UP  PULSED  ULTRASONIC  EQUIPMENT  AND  ANALOG  AND 
DIGITAL  SIGNAL  PROCESSING  SYSTEMS. 

3,  DEVELOPED  AND  TESTED  TECHNIQUES  FOR  THE  EXPERIMENTAL 
DETERMINATION  OF  THE  DISPERSIVE  BEHAVIOR  OF  COMPOSITE 
MATERIALS. 

A.  COMPLETED  PRELIMINARY  ULTRASONIC  DISPERSION  MEASURE¬ 
MENTS  FOR  ARBITRARY  ANGLES  OF  INCIDENCE  TO  THE 
REINFORCING  FIBERS  IN  UNFLAWED  COMPOSITE  SPECIMENS. 

5.  FORMULATED  A  METHOD  BASED  ON  THE  SCATTERING  MATRIX 
(T-MATRIX)  to  study  the  EFFECTIVE  PROPERTIES  OF  A 
FIBER-REINFORCED  COMPOSITE.  THE  STUDY  INVOLVES 

THE  MULTIPLE  SCATTERING  OF  SH-  WAVES  INCIDENT  NORMAL 
TO  THE  FIBERS  WHICH  MAY  BE  OF  ANY  ARBITRARY  CROSS 
SECTION. 

6.  EXTENDED  THE  T-MATRIX  FORMULATION  OF  ELASTIC  WAVES 
(P  AND  S)  TO  THE  MULTIPLE  SCATTERING  SITUATION  OF 
NORMAL  INCIDENT  WAVES  ON  PARALLEL  CYLINDERS  OF 
ARBITRARY  CROSS  SECTION. 
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FRLQUEKCY-DtPi-K'DENT  VELOCITY  fCASUREFiENTS 

A.  ANALOG  TECHNIQUES 

1.  Harmonic  Wave:  tt  -  phase  difference  measurements 


Display 

Oscilloscope 


Composi te 

_ _  Specimen 

Swept  r.  f. 

Osci  1  lator  - - i  x 


X-Y  Plotter 


r~ 


Wideband 
Ampl i f ier 


Frequency 

Counter 


amplitude,  phase 


Vector 

Voltmeter 


2.  Pulse:  group,  del  ay- time  meas'.'i-ements 


Pulser/ 
Receiver 
Sync  Source 

R  S 

Tronp.duccrr 

/_  f  X  1  - 


•®  Dual 
Verficof 
^  Amplifier 


[  Delaying/ 
Sync  L-p-c  Deloyed 
I  Time  Bose 


Specinic  n 


Oscilloscope 

Mainframe 


Digital 
^  Deloying 
Time  Base 


Display  J 


Delayed 
Time  Bose 


Delay  Interval  Oui 


Time-  i 

Interval  i 

Counter  t' 
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FREQUENCY-DEPtNDENT  VELOCITY  MEASUREMENTS 
B.  analog/digital  techniques 


-  ^1(5) 


0(ff)  =  ±  In 

c(f)  = 

fcomp  1  +  Rff) 
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rai  Forrni/fa  for  ibe  Scaiiering  of  Elasiic  Waves  TransHion  Matrix  for  Scafferfng  by  a  CyiirJer 


BIAXIAL  FATIGUE  OF  COMPOSITES 


P.  H.  Francis 


Southwest  Research  Institute 


PROGRAM  OBJECTIVES 


Determine  the  characteristics  of  first  ply  failure  (FPF)  of  laminated 
6r/Ep  composites  under  combined  states  of  plane  stress  (static  & 
-fatigue). 


.SUMMARY  OF  PRINCIPAL  FINDINGS 

1.  FPF"locallzed,  uniformly  distributed  matrix  failures  of  a  ply—is  a 
real  phenomenon.  Both  matrix  tension  and  shear  FPF  mechanisms 
have  been  observed  using  two  techniques. 

2.  FPF  predictions  can  be  very  sensitive  to  input  property  data. 

3.  FPF  generally  occurs  outside  predicted  envelopes,  suggesting 
laminate  theory  is  conservative  in  predicting  FPF. 

4.  FPF  is  a  precursor  to  laminate  failure,  i.e.,  measured  ultimate 
strength  always  exceeds  predicted  FPF. 

5.  FPF  data  so  far  suggest  only  weak  interaction  between  matrix  shear, 
tensile  modes. 


ORGANIZATION  OF  TALK 

1.  Tubular  specimens  for  characterizing  mechanical  performance 

of  laminated  composites. 

2.  Analytical  results  -  FPF. 

3.  Experimental  results  -  FPF. 

4.  Conclusions. 


209 


Exposed  View  of  Fabrication  Assembly 
for  Thin  and  Thick  Walled  Tubes 
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TUBES  VIS-A-VIS  FLAT  COUPONS  FOR  CHARACTERIZING  COMPOSITES 

RELATIVE  ADVANTAGES 


COUPON 

TUBE 

Geometric  Aspects 

_  _  . 

"thin”  (R/t  >  10) 

Edge  Effects 

influence  strength 

absent 

Grip  Aspects 

conventional  tabs 

rigid  tabs 
compensation 

Data  Value/Test 

limited 

full  orthotropic 
good  shear  specimen 

CosttoMfg.  +Tab 

inexpensive: 

cost  of  tube 

costly 

A  DEFINITION  OF 

TUBE  "THINNESS" 

REQUIREMENT 

R/t  for  given  %  stress  gradient  In  wall 

Loading  Mode 

5% 

10% 

Axial 

4 

2 

Hoop 

23 

13 

Torsion 

20 

10 

•Based  on  a  16-piy  (0/+  45)  Layup. 
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COMPARISON  OF  AVERAGE  LAMINATE  PROPERTIES 

MEASURED  USING  TUBES  AND  COUPONS 

PRO  PERTY  (0)^  TU  BE  DATA  COU  PON  DATA 

Exx 

19.3 

23.0 

jUt 

X 

— 

214 

.  (0)j 

"xy 

0.35 

0.34 

4 

^yy 

1.44 

1.6 

y 

3.5 

6.8 

•  ™8 

^yx 

0.03 

0.03 

0.94 

0.87 

.  {^(+451), 

s'" 

xy 

11.1 

10.8 

Fiberite  hy-E  1034C  graphite/epoxy  (T-300  fiber) 
Strength  values  In  ksi 
Modulus  values  in  10^  ksi 
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AXIAL  STRESS 


TORSIONAL  STRESS  (  ksi ) 

Ay 


Figure  Predicted  First  Ply  Failure  Surfaces  of  ( 0/90  Tubes  by  Laminate  Analysis 


TORSIONAL  STRESS  T^y  (  ksi ) 


Figure  Predicted  First  Ply  Failure  Surfaces  of  {  ±45  Tubes  by  Laminate  Analysis 
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TORSIONAL  STRESS  (  ksi ) 

Figure  FPF  Prediction  input  Data  Variation  {  0190)^ 


PLANE  SHEAR 
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^^6-  PfinctTarvI 
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P  P P  RfiLpl  \ c.o^  M  icrosco] 


AXIAL  STRESS 


^113  ksi 


TORSIONAL  STRESS  (  ksi ) 


Figure  Experimental  FPF  Measurements  on  (  0/90 Tubes 


TORSIONAL  STRESS  (  ksi } 


Experimental  FPF  Measurements  on  (  ±45  Tubes 
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Figure  . 


AXIAL.  STRESS^  ks.’ 


depend  ENT  FE^CTU^E  or  CoNPOS/TES 

P.  /?.  ^cha^e/y 

J^ercr/JTLY  C0MPLFrB£>  SruO/Es 

o  l^scoe/as/Y<^  Oyd^E  Ee /vy  Eiy/'Y^oEopE 

SL/ici  F/es’/rdi}y?eJ  Ee. 

0''  A  Thecry  0/  -  Depey c/^nr  l//scosAE/c  A^eJ/a 

/fe^/Ay  <y  /Ai^s  in  sd^/U  p/c^p>E/e>y'i  ayd 

^c/h^sd^f  coy')payei^'f‘  of  s^Yd/h^  /x/cfanj 

CO^R.P/Jr  PeSBA£cH 

o  ApproxiP?^-fe,  Method  of  Aoa/ysh  c/  Tir/je,  -  D^pendeyd  OacPp^ 
and  deJfy  for  prol/dns  f  wJ)icl  Phe  ExfeydeJ  C^rrespondeyce. 

PrAc/p/e.  do^as  joyd  App/y,  c^-pj  d^yyamic  prol/arnS'  f 

tyvohiyy  fef^p<2f^aii4/-e,  ^  rnyJs-/^/^  jiaJjbyPs) 

^ //(jynaricc^l  (A^Yery^AaMy  of  SEosJ  iPh/er>s}fy  EzLcYors  c/s/hy  Foda/ 

Pciy'Y  Y  f/ara  F E.  f  Eksi  Foor/er  Fyayj'fhr/>r  SoEd/oyiS^ 


Figure  1. 


Cruss-section  of  .'in  idoalizt^d  crack. 
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(psi-i) 


c/eAc^s  /AJ  o/srNoT^op/c  mg-p/a 


CrjL<k  S^eeck^'  P. 

*  K.:(,f)Ki  =4r 

C^fco  fic/ioyi  cf  oi  : 

Kr,  =o 


I^ni/tdh  'm  Ti/ue,  ti 

K.I  f  f  \1‘  It,,’  H-.-Vj^hr 


Figure  2.  Viscoelastic  functions  and  for  a  glass 
flbcr-roinforccd  plastic, 
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Spe&e/j 


fi  I  ^v  \ia..} 


f; 


C^/cu/i-^j^y\  o-^  ^ 


n 


^  -  32Tr  ^  Dt 

^oTo;' 


^/V4 


er/a. 


0.^=0 


Gj,=  CO 


If  K.  a  r<s'J 

'f  i<.sfs’rj/c'^-9j'^fc^v'V<,‘-"3'^ 

\»lne‘^fiv5K  Ki  IS  te-fwee^  +>'ei<2.  liVni'^’/ 
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EXPLORATORY  DEVELOPTOfT  OF  PROBABILISTIC  CHARACTERIZATION 
Am  TESTING  TECHNIQUES  FOR  FATIGUE  BEHAVIOR 
OF  COMPOSITE  LAMINATES 


Contract  No.  F33615-75-C-51J2 

BY 

Robert  A.  Heller 
Jann  M.  Yang 
John  J.  Mcgowan 

Department  of  Engineering  Science  and  Mechanics 
Virginia  Polytechnic  Institute  and  State  University 
Blacksburg,  Virginia 


PURPOSES: 


Examination  of  the  statistical  distributions  of  static  stre,ngth  and  stiffness 
FOR  GrAPHITE/EpoX'Y  composites  and  their  effects  ON: 


1.  Constant  Amplitude  Fatigue  Life 


2.  Random  Fatigue  Life 

3.  Residual  Strength  Reduction 
Stiffness  Changes 

5.  Creep  and  Creep  Rupture 

6.  Periodic  Proof  Tests 


Fig.  1 


Purposes 
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Autocorrolat 


Frequency,  jj  in  llz. 


18.  Output  spectral  density  (vnde  hand). 


rjE  0,  'jO,  ■- 

Panel  1 

«  Static  Tests  -  39.9  kf-.l 
o  .yroun  1  "i  '  ~  ‘tean  P  3  KMl:.  —  .^5  +  3(.1._) 

*  Til 

w  Croup  5  =  Mean  +  3  JO-IS  *=  .^5  +  3(,13) 

•  Croup  6  0^  -  Mean  +  3  Kt-tS  =  .40  i-  3 (.13) 

F  I'atieuc  Failure 

o  ♦ 

o  o 


o 

o 

F 


of  cvcle.s,  N  in  lU 


Moxinun  Stress  or  Residife.1  Strenqth,  o/o 


Group  10 

Static  Tests,  a  =  65.65  ksi 
Group  3,  0^  =  hiean  +  6  RMS  = 

Group  2.  Off,  =  Mean  +  3  RMS  " 

Group  4a ,  =  Mean  +  3  RMS  • 

Group  4,  Ofji  =  Moan  +  4  RMS  = 

Group  5.0^=  Mean  +  3  RMS  = 

Fatigue  Failure 


.41  +  5(.ll) 
.45  +  3(.12) 
.41  +  3(.09; 
.41  4  3(.09) 
.41  +  3(.12) 


C7E  0,  90,  +  4S/S 
Panel  2 


iUTtber  ui  cycle ' 


Fig.  20  Randoiii  Fatigue  Test  Results 


C/]i  0,  90  1  45, 


t  -  0  on  U/U 
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G/F.  0,  90,  i  Ub!y. 
Panel  1 

t  0  on  11/11/75 


♦  Statir 

*  Group  1 

•  Group  0 
o  Group  5 

•  Conrt.  Ampl 


1.0  4 


n/h  0,  i  45,  0/r. 
*  Static  Tests 
Constant  An:pl  . 
=  0  on  5/] 5/73 


Table  1  .  Creep  Tests  on 


G/E  0,+45,  0/5  Specimens,  0^  =  98.00  KSl 


cn  CU 

c  s~ 


D1  -r-  -r-  -r-  rS 

C  CO  -O  I/:  4J 


DATE 

SPECIMEN 

NO. 

STRESS 

KSI 

STRESS 
%  OF 

•r-  0. 

■0 

rrj  UD 

0  0 

UJ  X 

r3  D- 

0 

r—  VO 

C  0 

3  i~ 

LU  X 

Temoera 

REMARKS 

8/12/75 

G3I117 

80 

8.35 

200 

y^.ns 

G3I1C9 

85 

10.26 

200 

10/27/75 

G5I6S 

78.00 

80 

10.86 

80 

10/27/75 

G5I33 

78.00 

80 

10.00 

80 

11/7/75 

G5I12B 

8.15 

80 

Broke  durinq  loadinq 

11/10/75 

G5I133 

88.20 

90 

10.00 

80 

11/10/75 

G5I150 

88.20 

90 

10.00 

30 

n/12/75 

G5I71 

83.30 

85 

10.16 

80 

n/19/75 

G5I6] 

78.40 

80 

9.34 

200 

Temperature  controls  incorrect 

12/16/75 

G5I57 

78.40 

80 

10.74 

200 

1/14/76 

G5I21 

83.30 

85 

10.19 

10.19 

ion 

zf^ne 

G5I21 

88.20 

90 

10.00 

10.00 

200 

Stress  increased  to  90% 

3/23/76 

65121 

93.10 

95 

10.00 

10.00 

200 

Stress  increased  to  0^% 

2/4/76 

65122 

83.30 

85 

9.52 

2C0 

Temperature  controls  incorrect 

2/25/76 

65 146 

83.30 

85 

10.00 

200 

Tempera turc  controls  incorrect 

3/30/75 

65115 

83.30 

85 

10.13 

10.13 

.200 

4/7/76 

Gbiij 

88,20 

90 

9.33 

9.80 

Stress  increased  to  90% 

4/5/76 

65163 

83.30 

85 

11.54 

11  .18 

200 

4/19/75 

G5I63 

38.20 

90 

n  .33 

11  .32 

200 

Str'ess  increased  to  90% 

^/27/ie 

G5I53 

93.10 

95 

11  .15 

200 

Stress  Increased  to  95% 

DATE 


Table  2  ,  Creep  Tests  on 

G/E  0, +45, 90/5  Specimens,  ==  68.38  KSI 
Proof  Stress  =  63  KSI 


SPECIMEN 

STRESS 

STRESS 
%  OF 

NO. 

KSI 

o> 

c: 


ZJ 


UJ  X 


a» 

i- 

3 

+j 


o 

QJ 

O- 

E 


<u 

r- 


REMARKS 


1/26/76 

2340^7 

58.13 

85 

7.5 

7.5 

80 

2/5/76 

2340-77 

61.55 

90 

7.5 

7.5 

80 

Stress  increased  to  90% 

2/10/7C 

2340.77 

64.97 

95 

7.5 

7.5 

80 

Stress  increase.'i  to  95% 

1/27/76 

2336=10 

80 

Broke  during  loading 

1/29/76 

2374^2 

80  - 

Broke  during  loading 

zf7ne 

2340=17 

58.13 

85 

7.35 

7.35 

80 

2/24/76 

2340=17 

61 .97 

90 

7,35 

7.35 

80 

Stress  increased  to  90% 

3/23/76 

2340?! 7 

64.96 

95 

7.31 

7.31 

80 

Stress  increased  to  95% 

2/10/76 

2.’61=5 

54.71 

8C 

6.10 

80 

Broke  after  one  day 

3/29/76 

2337=22 

58.13 

85 

6.54 

6.36 

80 

4/18/76 

2337=20 

61.55 

90 

6.27 

6.10 

80 

5/5/76 

2337=20 

64.97 

95 

6.20 

80 

Stress  increased  to  95%. 

5/3/76 

2336=7 

5.88 

80 

Broke  during  proof  loading 

5/14/76 

2374=3 

6.71 

80 

Broke  during  proof  loading 

5/17/76 

2303+17 

58.13 

85 

5,63 

5.58 

200 

Proof  loaded  at  80®F.  E=5.69 

5/19/V6 

2326=16 

58.13 

85 

5.27 

6.27 

80 

6/10/76 

?335;7’9 

58.13 

85 

6.25 

200 

Proof  loaded  at  80'’ F.  E=6.98 

6/8/76 

2375+17 

58.13 

85 

6.92 

6.90 

.80 

7/27/76 

2375=17 

61  .55 

90 

6.90 

6.74 

80 

Stress  increased  to  90% 

6/10/76 

2325fl 

61.55 

90 

7.00 

6.85 

80 

7/18/76 

2325=1 

64.97 

95 

80 

Broke  during  loa-J  increase 

6/10/76 

2340=17 

58.13 

85 

7.40 

7.15 

200 

Proof  loaded  at  80'’ F.  E=7.:'’ 

7/27/76 

2325--20 

58.13 

85 

6.4? 

80 

RrekP  during  loadinq 

7'27/76 

2375=18 

58.13 

85 

6.0 

80 

‘'’roi'.'.  du^Mng  loadi''  ) 
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Table  3  .  Creep  Test^  on 


DATE  SPECIMEN.  STRESS 
NO.  KSI 


G/E  0, 

90.+45/s 

Specimens . 

Panel  1 

Qjj  =  59.57  KSI 

STRESS 
%  OF 

^u 

C71  -r- 

C  c/) 
•r-  a. 

-o 

fO 

o  o 

UJ  X 

E  unloading 

X  10^  psi 

Temperature 

op 

REMARKS 

12/29/75 

5-11 

50.63 

1/15/76 

5-11 

53.61 

1/21/76 

12/22/75 

4-14 

50.63 

1/16/76 

4-14 

53.61 

85 

6.11 

6.11 

80 

90 

6.11 

6.11 

80 

6,20 

6.20 

80 

35 

7.06 

7.06 

80 

90 

6.92 

6.92 

80 

Stress  increased  to 
Specimen  broke  during  load 
increase 

Stress  increased  to  90;;! 
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Date 

Specimen 

»/■> 

r3 

3 

-o  tn 
.r-  c> 

ro 

o 

X 

•O 

O)  V. 

•>-  0) 

'a.  Tj 

No. 

o  +-► 

CL 

c;  oo 

2-11-7C 

4-n 

58.48 

526.30 

2-10-76 

3-8 

. 

2-3-76 

5-6 

- 

73.00 

2-5-76 

2-1 

• 

12.40 

2-4-76 

1-10 

59.43 

427. 5D 

2-3-76 

6-13 

2.00 

1  -18-76 

2-3 

64.96 

787.00 

1-23-76 

5-5 

- 

.640 

1  -22-76 

6-12 

- 

80.59 

1-21-76 

6-11 

51.11 

500.00 

1-19-76 

1-13 

- 

11  .34 

1-16-76 

1-14 

353.70 

1-15-76 

5-8 

- 

184.70 

1-14-76 

2-3 

65.03 

423.10 

1-13-76 

6-lC 

•• 

4.07 

11-13-/6 

4-3 

1.12 

TABLE  4  RAMDOM  FATIGUE  TEST  RESULTS 

G/E  0,90, +45/  Panel  No.  1  Group  No.  1 
Mearr  Stress  =  26.80  ksi 


RMS  Stress 

=  6.95 

ksi 

Clipping  at 

Mean 

+3  RMS 

4-> 

A-> 

•M 

fd 

rt3 

ra 

fO 

•r~  t/1 

•r-  lA 

■*-  V) 

•r-  10 

l/S  (U 

t/1  <U 

<A  CJ 

to  QJ 

fO 

D-  I— 

CL  r- 

CL  r-- 

a. 

♦->  tA 

VJD  VD 

4D  >, 

to  >, 

•1-  O 

O  O 

O  u 

o  u 

O  (.J 

►  4  VO 

n 

fO 

n 

ro 

O 

X  o 

X  o 

X  c> 

X  o 

o 

r—  * 

CSJ 

ro 

'S’ 

UJ  X 

UJ  X 

I.LJ  X 

UJ  X 

UJ  X 

6.21 

6.00 

5.89 

5.79 

5.79 

at 

at 

at 

at 

1.00 

5.00 

10.00 

56.49 

5.89 

6.32 

6.18 

5.97 

5.75 

at 

at 

at 

2.00 

8.62 

50.00 

5.89 

5.89 

5.89 

5.80 

at 

at 

at 

1 .00 

5.00 

10.00 

6.22 

6.20 

6.04 

5.8b 

5.84 

at 

at 

at 

at 

1 .00 

5.00 

io.ee 

34.90 

6.02 

6.22  5.90 

5.90 

5.73 

5.72 

at 

a  t 

at 

at 

1.00 

5.95 

5.00 

50.00 

787.00 

6.09 

6.01 

6.01 

5.82 

at 

at 

at 

1.00 

5.00 

10.00 

6.14 

5.89 

5.80 

5.76 

5.54 

at 

at 

at 

at 

10.00 

50.00 

100.00 

600.00 

6.31 

5.98 

5,98 

at 

at 

1 .00 

5.07 

5.81 

5.66 

5.65 

5.38 

5.28 

at 

at 

at 

at 

3.00 

20.00 

50.00 

100.00 

5. S3 

5.61 

5.47 

5.30 

at 

at 

at 

6.02 

50.00 

100.00 

6.37 

6.07 

6.01 

5.88 

5.87 

at 

at 

at 

at 

1,00 

6.36 

70.93 

423.00 

6.10  5.88 

at 

2.00 


■M 

•r-  \A 
iTt  OI 

n.  f- 
o 

lO  >, 

o  u 
ro 

Remarks  (FAILURE  MODE) 

UJ  X 

5.79  5.69  X  10^  at  526.3  x  10^ 

at  cycles.  Broke  in  middle 

100.00  during  tension  test 

Broke  while  proof  tested 

Broke  in  middle  during 
fatigue 


Broke  in  lower  tab 
du»'in9  fatigue 


5.77  Broke  in  middle  during 

at  tension  test 

427.50 

Broke  in  middle  ( doubt f'j* 
data)  during  fatigue 

Broke  in  middle  during 
tension  test 


Broke  near  lower  tab 
during  fatigue 

Broke  in  top  grip 
during  fatigue 


Broke  in  top  grip 
during  tension  test 


Broke  in  top  grip 
during  fatigue 


Broke  in  top  grip 
during  fatigue 


Broke  in  middle 
during  fatigue 


Broke  in  lower  grip 
during  tension  test 


Brokt.'  in  to;.-  grip 
dur'ing  fatigue 


Broke  in  top  grip 
during  fatigue 
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TABLE  5  RANDOM  FATIGUE  TEST  RESULTS 

G/E  0,90, +45/  Panel  No.  2  Group  No,  2 
Mean  Stress  =  29.54  ksi 
RMS  Stress  =  7.66  ksi 
Clipping  at  Mean  +^3  RMS 


4J  +J  4J  4->  4-> 

(0  fl3  «tJ  »0  *0 


Date 

Specimen 

CO 

o 

CO 

-O  cO 

CJ 

CO  I_ 

o 

r— 

X 

•o 

O  Cl 
•I-  <3J 

*Q.  Tj 

<a 

*.>  to 
•f  CL 

c 

i-r  to 

O 

•r*  CO 

CO  O 

CL  p- 
U 

to  >> 

O  U 

m 

X  o 

cri 

CO  at 

Cl,  p— 

u 

VO  >> 

o  u 

ro 

X  o 

♦r-  CO 

t/>  OJ 

CL  P- 
<J 

VO  >, 

O  'w) 

CO 

X  o 

•p-  Cl 

Cl  at 

D.  »— 

Cl 

to  >. 
o  o 

fO 

X  o 

•P-  Cl 

Cl  a> 

CL  p— 

u 

to  >, 
o  o 

ro 

X  o 

un 

No. 

■>-* 
c:  to 

CL  >, 

C  O 

LU^  X 

Li7~  X 

X 

LJ^  X 

Ut  X 

UJ  X 

3-30-76 

9-5 

84.37 

100.00 

7.m 

6.88 

at 

6.80 

at 

6.50 

at 

G.50 

at 

6.50 

at 

1.00 

5.00 

10.00 

50.00 

100.00 

2-24-76 

12-1 

76.88 

100.00 

6.74 

6.74 

at 

6.74 

at 

6.74 

at 

6.74 

at 

6.74 

at 

1.00 

5.00 

10.00 

50.00 

100.00 

2-19-76 

12-3 

64.45 

100.00 

6.69 

6.52 

at 

6.52 

at 

6.52 

at 

6.52 

at 

6.52 

at 

1.00 

5.00 

10.00 

50.00 

100.00 

2-17-76 

1-4 

11.64 

6.43 

6.27 

6.03 

G.03 

at 

at 

at 

1.00 

5.00 

10.00 

2-18-76 

13-5 

69.50 

330.00 

6.49 

6.33 

at 

6.30 

at 

6.30 

at 

6.30 

at 

6.30 
at  ’ 

1.00 

5.00 

10. on 

50.00 

330.00 

2-17-76 

14-1 

- 

.85 

6.30 

2-16-76 

7-6 

65.41 

100.00 

6,85 

6.66 

at 

1.00 

6.44 

at 

5.00 

5.44 

at 

10.00 

6.34 

at 

50.00 

6.27 

at 

100.00 

2-12-76 

4-2 

63.47 

100.00 

6.51 

6.39 

at 

1.00 

6.8C 

at 

5.00 

6.42 

at 

10.00 

6.27 

at 

50.00 

6.27 

at 

100.00 

TABLE  6  RANDOM  FATIGUE  TEST  RESULTS 

G/E  0,90 ,+45/  Panel  No.  2  Group  No.  3 
Mean  Stress  =  22.93  ksi 
RMS  Stress  7.02  ksi 
Clipping  at  Mean  +_5  RMS 


4^  4-> 

lO  ^  ^ 


CO 

a 

Cl 

Cl 

<u 

•P-  o 

o  ai 

o 

•/I  a» 

(/» 

Cl 

OJ 

Cl 

v> 

Q> 

in 

PH 

a. 

*— 

CL  p— 

O  1 — 

p. 

Ti 

fL 

>C 

VO 

to  ^ 

VO  ^ 

VO 

>, 

VO 

*>1 

u 

O. 

o 

o 

o  o 

o  o 

o 

u 

O 

o 

a 

Cl 

a» 

c 

t£> 

O 

p— 

ro 

CO 

r  — 

. 

Date 

.Specin'M-ti 

T3 

Cl 

0) 

t. 

CL 

QJ 

U 

27. 

X 

o 

X  o 

X  o 

X 

‘  o 

X 

o 

No. 

S 

CH 

CO. 

o. 

c 

F? 

o 

UJ  X 

UJ~  X 

CVJ 

UJ  X 

CO 

UJ  X 

UJ  X 

VO 

lO 

X 

4-l-7i 

1-6 

53 

o 

o 

6 

.47 

6, 

.28 

6.28 

6.28 

6. 

28 

at 

at 

at 

at 

1, 

.00 

5.00 

10.00 

50. 

00 

Remarks  (FAILURE  MODE) 


br'oke  in  middle 
during  tension  test 


Broke  in  top  grip 
during  tension  test 


Broke  in  lov/er  grip 
during  tension  test 


Broke  in  top  grip 
during  fatigue 


Broke  in  middle 
during  tension  test 


Broke  in  middle 

Specimen  damaged  in  loading 

Broke  in  middle 
during  tension  test 


Broke  in  middle 
during  tension  test 


Remarks  (FAILURE  MODE) 

Broke  in  middle 
during  fatigue 
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TABLE  7 

RAriDO:-:  i 

FATIGUE 

lES'i  RESULTS 

G/E  0,90, +45/  Panel  No. 

2  Group  No.  4 

Kean  Stress 

=  26.80  ksi 

RMS 

Stress 

=  5.96 

ksi 

Clipping  at 

Mean  +4  RMS 

4-> 

4^ 

+j 

4-> 

4J 

<o 

to 

fO 

<0 

fO 

ro 

•  r-  trt 

•r-  KA 

,r-  (/I 

♦r-  t/l 

V)  QJ 

t/) 

o 

ro 

CL 

CL  f— 

CL  r— 

Q.-  ^ 

(j 

a.  r- 
o 

X 

Ul  VD  >,  LD  >> 

•43  >. 

O 

lO  >, 

o  o 

LD  >, 

o  a 

Date 

Spacimen 

■3  trt 

XJ  «/' 

•r-  OJ 
t/)  ^ 

“O 

a;  to 

■T-  <D 

c 

w_  VO 

o 

ro 

X  o 

fO 

X  o 

CO 

X  o 

X  o 

ro 

X  o 

Remarks  (FAILURE  MODE) 

No. 

CJ  4-» 

ca  to 

CL  >. 

<c  o 

o 

Ul  X 

uT”  X 

CVI 

LU  X 

ro  ^ 

UJ  X 

UJ  X 

un 

UJ  X 

5-27-76 

13-1 

76.88 

677.00 

7.33 

7.13 

7.13 

7.09 

6.94 

7.17 

Broke  in  middle 

at 

at 

at 

at 

at 

during  tension  test 

1.00 

5.00 

10.00 

50.00 

677.00 

5-20-76 

n-5 

75.15 

508.70 

6.39 

7.22 

7.05 

7.05 

7.05 

7.05 

7.05x10^  at  503.7x103  cycles 

at 

at 

at 

at 

at 

Broke  In  top  tab 

1 ,00 

5.00 

10.00 

50.00 

96.50 

during  tension  test 

5-20-76 

2-5 

18.00 

6.88 

6.8y 

6.88 

6.88 

Broke  near  top  tab 

at 

at 

at 

during  fatigue 

1.00 

5.00 

10. CO 

5-13-76 

11-3 

78.10 

511 .00 

7.13 

6.90 

6.81 

6.90 

6.81 

6.97 

6.97x10^  at  511x103  cycles 

at 

at 

at 

at 

at 

Broke  in  lov/er  tab 

1  .00 

5.0C 

10. cc 

50.  CO 

100.00 

during  tension  test 

4-29-76 

5-6 

75.45 

517.50 

7.43 

7.10 

7.10 

7.31 

6.99 

6.96 

6.91x10^  at  517.5x103  cycles 

at 

at 

at 

at 

at 

Broke  in  middle 

1.00 

5.00 

10.00 

50  00 

81  .00 

during  tension  test 

5-8-76 

9-2 

75.76 

501 .00 

7.20 

7.16 

7.21 

7.01 

7.01 

7.06 

7.06x10®  at  501x10'^  cycles 

at 

at 

at 

at 

at 

Broke  in  top  tab 

1 .00 

5.00 

10.00 

50.00 

100.00 

during  tension  test 

4-22-76 

5-5 

500,00 

6.76 

6.76 

6.76 

6.76 

6.76 

6.76’ 

6.49x10^  at  500x10-^  cyclr-; 

at 

at 

at 

at 

at 

Broke  accidentally 

1,00 

5.00 

10.00 

50.00 

100.00 

Doubtful ,data 

4-16-76 

12-4 

63.4/ 

500,00 

7.03 

7.03 

7,33 

7.18 

7.18 

7.03 

7.03x10®  at  500x10^  cycles 

at 

at 

at 

at 

at 

Broke  in  lov/er  tab 

1.00 

5,00 

10.00 

50.00 

100,00 

during  tension  test 

4-15-76 

10-3 

- 

86.50 

7.22 

6.99 

6.99 

G.99 

6.99 

Possible  machine 

at 

at 

at 

at 

fiial  function 

1.00 

5.00 

10.00 

50.00 

4-15-76 

10-3 

86.50 

7.22 

6.99 

6.9y 

6.99 

6.99 

Broke  in  middle 

at 

at 

at 

at 

Machine  malfunction 

1 .00 

5.00 

10.00 

50.00 

during  fatigue 

4-13-76 

9-] 

- 

486.00 

7.20 

6.90 

6.90 

6.90 

6.90 

6.90 

Broke  at  bottom  grip 

at 

at 

at 

at 

at 

during  fatigue 

1 .00 

5.00 

1C.  00 

50.00 

lOO.Oo' 

4-8-76 

13-2 

- 

54.70 

6.59 

6.59 

6.47 

6.40 

6.30 

Broke  in  lower  tab 

at 

at 

at 

at 

during  fatigue 

1 .00 

5  ,  T'; 

10.00 

50.00 
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TABLE  8  RANDOM  FATIGUE  TEST  RESULTS 

G/E  0,90, +45/  Panel  No.  2  Group  No.  4A 
Mean  Stress  =  26.80  ksi 
RMS  Stress  =  5.96  ksi 
Clipping  at  mean  +3  RMS 


+-» 

fO 

4J 

nj 

to 

4-» 

to 

n 

o 

in  Qj 

in 

m  0) 

•1-  in 
in  o 

♦r-  to 

in  Oi 

to 

to  <y 

X 

Tj 

<c 

4-* 

•r~  ci. 

CL  1— 

u 

uo  >, 

o  u 

O- 

u 

>> 

O  u 

D.  r— 

O 

>» 

O  O 

Ol  *~ 
o 

to  >, 

o  o 

O.  r— 

o 

VO  >, 

o  u 

Date 

Specimen 

3  </) 

TJ 

■r-  0) 

t/1 

a»  «/> 

•v  <u 

o 

c 

)  \£> 

^  O 

X  o 

ro 

X  o 

CO 

X  o 

ro 

X  o 

CO 

X  o 

Remarks  (FAILURE  MODE) 

No. 

CJ  4-> 

cc  lo 

D.  >> 

<X. 

o 

Ul  X 

uj”  X 

e\» 

LU  X 

fO 

UJ  X 

uT  X 

ulf'  X 

6-9-76 

8-2 

73.71 

500.00 

7.16 

7.45 

7.31 

7.18 

7.36 

7.37 

Broke  in  top  grip 

at 

at 

at 

at 

at 

during  tension  test 

1.00 

5.00 

10.00 

50.00 

500.00 

6-7-76 

2-4 

64.90 

6.33 

6.27 

6.43 

6.38 

Broke  in  lower  tab 

at 

at 

at 

due  to  twisted  torque 

1.20 

5,00 

13.80 

cell  during  fatigue 

6-8-76 

3-5 

69.65 

503.00 

6.87 

6.93 

6.71 

6.60 

6.64 

6.08 

6.64x10®  at  503x10^  cycles 

at 

at 

at 

at 

at 

Broke  in  top  tab 

1.00 

5.20 

10.00 

50.00 

100.00 

during  tension  test 

6-1 5-76 

11-1 

73.54 

614.00 

7.00 

6.81 

6.81 

5.62 

6.62 

6.62 

Broke  in  middle 

at 

at 

at 

at 

at 

during  tension  test 

2.60 

5.00 

10.00 

329.00 

614.00 

6-17-76 

11-6 

7Z.n2 

500.00 

7.27 

6.70 

6.82 

6.72 

6.92 

6.90 

6.90x10^  at  500x10^  cycles 

at 

at 

at 

at 

at 

,  Broke  in  lov/er  tab 

1 ,00 

5.00 

10.00 

50.  OU 

100.00 

during  tension  test 

6-21-76 

6-3 

500.00 

7.55 

7.17 

7.17 

6.79 

6.79 

6.79 

6.85x10^  at  500x103  cycles 

at 

at 

at 

at 

at 

Broke  accidentally 

1.00 

5.00 

14.70 

53.80 

102. CU 

Doubtful  data 

6-22-70 

9-4 

72.47 

452.00 

6.91 

7.11 

7.10 

7.10 

7.02 

Broke  in  top  tab 

at 

at 

at 

at 

during  tension  test 

1.00 

5.00 

14.90 

452.00 

6-25-76 

4-3 

83.22 

555.00 

7.53 

7.53 

7.21 

7.13 

7.33 

7.45 

7.13x10®  ut  555x10^  cycles 

at 

at 

at 

at 

at 

Broke  near  lower  tab 

1.00 

5.00 

10.00 

66.90 

109.00 

during  tension  test 
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TABLE  9  RANDOM  FATIGUE  TEST  RESULTS 

G/E  0,90, +45/  Panel  No.  1  Group  No.  5 
Mean  Stress  =  26.80  ksi 
RMS  Stress  =  7.94  ksi 
Clipping  at  Mean  j\3  RMS 


«/» 

■U. 


=}  t/l 

-a  k/> 


Date 

Specimen 

No, 

QJ 

</)  S- 
0)  .M 

a:  m 

1-15-76 

5-4 

- 

n-12-75 

2-10 

- 

n-11-75 

3-2 

0.0b 

2.80 

1.68 


6.13 

6.48 

6.30 


Remarks  {FAILURE  MODE) 


Broke  in  top  tab 
during  fatigue 

Broke  in  lov/er  tab 
during  fatigue 

Broke  in  top  tab 
during  fatigue 


11-11-75 

3-13 

“ 

3.18 

6.12 

5.94 

at 

1.68 

Broke  in  top  tab 
during  fatigue 

12-5-75 

4-4 

- 

0.10 

6.11 

Broke  in  middle 
during  fatigue 

12-10-75 

4-5 

57.79 

803.00 

5.30 

at 

803.00 

Broke  in  top  tab 
during  tension  test 

12-12-/E 

6-3 

65.41 

693.00 

6.27 

6.00  5.75 

at  at 

1.11  32. CO 

5.73 

at 

693. CG 

Broke  in  middle 
during  tension  test 

11-10-75 

6-7 

- 

70.00 

6.08 

Broke  in  top  tab 
during  fatigue 

11-11-75 

4-6 

- 

85.00 

6.3' 

6.04 

at 

1.00 

Broke  in  lower  tab 
during  fatigue 

TABLE  10  RANDOM.  FATIGUE  TEST  RESULTS 

6/E  0,90.+45/  Panel  No.  1  Group  No.  6 
Mean  Stress  =  23.83  ksi 
RMS  Stress.  =  6.95  ksi 
Clipping  at  Mean  +3  RMS 


+.»  +j  +j 

<o  «  OJ 


Date 

Specimen 

No. 

■O  m 

.>•  o 

l/l 

QJ  +J 

a:  in 

•r-  01 

"q-  Tj 

CL  >> 

o 

—1  VO 
w.  O 

o 

LU  X 

ro 

X  o 

Lu""  X 

1-8-76 

2-7 

55.93 

116.00 

8.03 

8.03 

t/l 

t/»  o 


8,03 


•f-  in 
tn  Oi 
O. 

u 

lO  >, 
O 


at  at 
2.40  116.00 


•r-  in 

in  <b 

Ol  r— 

u 


ro 


o. 

o 

L  >> 


Remarks  (FAILURE  MODE) 

Broke  in  top  tab 
during  tension  test 
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TABLE  11  RANUOM  FATIGUE  TEST  RESULTS 

G/E  0,90, +45/  Panel  No.  2  Group  No.  10 
Mean  Stress  -  26.80  ksi 
RMS  Stress  =  7.4  ksi 
Clipping  at  Mean  +^3  RMS 


■M 

<o 

4-J 

TO 

TO 

TO 

TO 

_ _ _ 

CO 

•r-  tA 

■r- 

•*“  iA 

•r-  «/ 

‘I-  tA 

•r- 

o 

r— 

(/>  a» 

V)  <U 

iA  01 

tA  a 

tA  QJ 

r- 

.|2  .r- 

u 

O 

*0 

*0 

to 

^  a 

O  l/ 

o  <y 

o  u' 

o  o' 

o 

3  tfl 

(U  g 

c 

(O 

»■—  ^ 

^  m 

Date 

Specimen 

•r-  a 
iA  ^ 

'q.  SS 

-  o 

X  o 

X  o 

X  o 

X  o 

X  o 

Remarks  (FAILURE  MODE 

No . 

cc  in 

LJ°  X 

uj”  X 

UJ^  X 

up  X 

up  X 

up  X 

7-20-76 

10-4 

72.05 

424.00 

6.73 

6.73 

6.73 

6.73 

6.73 

Broke  in  middle 

at 

at 

at 

at 

during  tension  test 

1.00 

5.00 

10.00 

423,59 

7-19-76 

7-1 

70.21 

759.00 

6.62 

6.43 

6.62 

7.00 

7.00 

Broke  in  middle 

at 

at 

at 

at 

during  tension  test 

1.00 

5.00 

10.00 

759.00 

7-21-76 

7-2 

66.22 

500.00 

6.3/ 

6.37 

6.37 

6.58 

6.74 

6.37 

Broke  in  middle 

at 

at 

at 

at 

at 

during  tension  test 

1,00 

5.00 

10.00 

47.00 

500.00 

7-23-76 

9-3 

71.63 

802.00 

6.70 

6.70 

6.51 

6.51 

6.51 

Broke  in  middle 

at 

at 

at 

at 

during  tension  test 

1.00 

5.00 

10.00 

802,00 

8-5-76 

4-5 

76.80 

500.00 

6.65 

6.65 

6.65 

6.54 

6.54 

6.54 

Broke  in  top  tab  • 

at 

at 

at 

at 

at 

during  tension  test 

1.00 

5.CU 

10.00 

50,  OC 

SCO. CO 

8-3-76 

6-6 

77.16 

784.00 

7.02 

6.72 

6.72 

6.50 

6.71 

Broke  near  lov/er  tab 

at 

at 

at 

at 

during  tension  test 

1.00 

5.00 

10.  CO 

734.0-.} 

TABLE  12  RANDOM  FATIGUE  TEST  RESULTS 

G/E  0,90, +^45/  Panel  No.  2  Group  No.  11  Wide  Band 
Mean  Stress  -  26.80  ksi 
RMS  Stress  =  5.96  ksi 
Clipping  at  Mean  +  4  RMS 


Remarks  (FAILURE  MODE) 


\o 

t 

1 

7-4P2 

79.78 

500 

6.65 

6.65 

6.96 

6.96 

6.5/ 

6.49 

6.49 

Broken  in  middle 

at 

at 

at 

at 

at 

at 

1 

5 

10 

50 

100 

500 

9-13-76 

10-6P2 

85.03 

500 

6.93 

6.93 

6.93 

6.74 

6.74 

Broken  ‘in  middle 

at 

at 

at- 

at 

1 

5 

10 

500 

9-9-76 

4-1P2 

78,23 

500 

6.84 

6.65 

6.65 

6.65 

6.65 

6.65 

6.65 

Broken  in  middle 

at 

at 

at 

at 

at 

at 

1 

5 

10 

50 

100 

SCO 

9-7-76 

12-6P? 

79.34 

510 

6.67 

6.64 

6.64 

6.45 

6.82 

6.42 

Broken  in  middle 

at 

at 

at 

at 

at 

1 

5 

10 

50 

510 
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TABLE  13  RANDOM  FATIGUE  TEST  RESULTS 
G/E  0,90, +45/  Panel  No.  2  Group  No.  14  (Wide  Band) 
Mean  Stress  =  26.30  ksi 
RMS  Stress  =  7.40  ksi 
Clipping  at  Mean  +  3  RMS 


8-10-76  6-*lP2  76.40  501 

813-76  6-5P2  77.92  501 

8-13-76  4-4P2  74.00  500 

8-19-76  5-3P2  78.27  461 

8-24-76  13-6P2  65.44  803 

8-26-76  6-2P2  75.21  501 


6.79 

6.79 

at 

6.94 

at 

6.32 

1 

5 

10 

6.58 

6.51 

7,03 

6.66 

at 

at 

at 

1 

5 

10 

6.63 

6.55 

G.63 

6.48 

at 

at 

at 

1 

5 

50 

6.69 

6.58 

6.61 

6.69 

at 

at 

at 

1 

5 

11 

6.77 

6.39 

6.39 

6.39 

at 

at 

at 

1 

5 

10 

7.03 

6.91 

6.91 

6,84 

at 

at 

at 

1 

5 

104 

6.52  Broken  in  middle 

501 

C.58  6.51  i^'oken  in  middle 

at  at 
7.19  501 

6.26  Broken  In  middle 

at 

500 

6.73  Broker,  in  middle 

at 
461 

6.39  Broken  in  middle 

at 
303 

6.46  Broken  in  middle 

at 

501 
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TABLf;  14  .  CONSTANT  AMPLITUDE:  FATIGUE  TEST  RESULTS 


G/E  [U,90+45]j^  Panel  1 

Average  Ultimate  Strength  =  59.566  KSI 
Stress  Ratio  >  0.1,  Cyclic  Frequency  =  10  Hz 


Mean 

Maximum 

Res  idual 

'"Applied 

Iq  E, 

(Initial)  , 

E,  E. 

Specimen 

Stress 

Stress 

Stress 

Cycles 

^  ^  Failure 

Number 

(KSI) 

(KSI) 

(KSI) 

(103) 

(10^  psi )  {10^  psi ) 

(10°  psi)  (10°  psi)  Mode 

End  tab 


11-11-75 

0 

5-3 

29.03 

52.70 

0 

3.01 

11-20-75 

9 

5-1 

27.85 

50.63 

0 

0.22 

11-18-75 

8 

3-9 

27.85 

50.63 

0 

0.48 

12-18-75 

37 

5-13 

27.85 

50.63 

0 

0.16 

12-17-75 

36 

6-4 

27.85 

50.63 

0 

0.20 

11-20-75 

9 

6-1 

26.21 

47.65 

0 

0.26 

1-12-75 

61 

2-4 

26.  n 

47.65 

0 

0.54 

11-20-75 

c 

6-2 

26.21 

47.65 

0 

0.40 

1-17-76 

61 

3-10 

26.21 

47.65 

0 

0.46 

n-20-75 

9 

4-2 

24.57 

44.67 

0 

93.63 

n-21-75 

10 

1-6 

24.5/ 

44.67 

0 

0.92 

1-2-76 

51 

2-9 

24.57 

44.67 

0 

1.69 

1-6-76 

55 

3-3 

24.57 

44.67 

0 

27.21 

1-5-76 

54 

2-8 

24.57 

44.67 

0 

527.84 

1-5-76 

54 

5-9 

24.57 

44.67 

0 

403.17 

12-21-75 

100 

2-12 

22.93 

41.70 

0 

2.04 

12-1-7S 

20 

2-11 

22.93 

41.70 

0 

2.74 

1-2-75 

51 

3-11 

22.93 

41 .70 

58.75 

945.10 

1-2-76 

51 

1-9 

22.93 

41  .70 

65.46 

1124.88 

12-30-75 

49 

1-7 

21.30 

38.72 

59.78 

110.00 

12-30-75 

49 

2-2 

21.30 

38.72 

53.48 

1686.29 

12-22-75 

41 

3-12 

21.30 

38.72 

43.76 

1000.00 

12-23-75 

42 

5-12 

21.30 

38.72 

53.25 

1000.00 

11-8-75 

57 

C-8 

21.30 

38.72 

0 

39.39 

11-8-75 

57 

6-9 

21.30 

38.72 

62.40 

1003.50 

8.23 

End  tab 

End  tab 

6.35 

6.32 

n=100 

End  tab 

End  tab 

7.89 

End  tab 

6.66 

End  tab 

7.07 

End  tab 

6.27 

End  tab 

7.30 

7.06 

n=5,000 

8.28 

n=24,10n 

End  tab 

7.41 

End  tab 

6.10 

End  tab 

6.51 

6.43 

n=l,000 

6.25 
n--- 4,000 

End  tab 

6.18 

6.07 
n=2 ,000 

5,22 

259.240 

End  tab 

5.11 

4.93 
at  5000 

End  tab 

6.91 

End  tab 

7.14 

Fnd  tab 

6.45 

6.32 

n=2,b00 

6.26  6.26 
n=5,300  n=945,lCG 

End  tab 

Near  middle 

6.18 

6.02 

n=5,000 

6.02  5.99  End  tab 

n=50,000  n=110,000 

6.04 

5.72 

n=5,000 

5.71  5.31 

n=l 0,000  n -66, 000 

End  tab 

5.36 

4.71 

n’^5,000 

4.57  5.61 

n-:85.000  n=l,000,000 

End  tab 

4.92 

4.91 

n-5,000 

6.34 

1,000, 000 

End  tab 

6.26 

6.00 

n=5,000 

End  tab 

6.49 

6.45 

n^5,000 

6.42 

n-1  ,003,500 

Near  middle 
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TABLE  15  .  CONSTANT  AMPLITUDE.  FATIGUE  TEST 
Type  of  Specimen  [0,+^45,90]2,  Stress  Ratio=0.1 


Date 

(Days) 

Specimen 

Number 

Mean 

Stress 

(KSI) 

Maximum 

Stress 

(KSI) 

Residual 

Stress 

(KSI) 

"Applied 

Cycles 

(103) 

Eo 

(Initial ) 
(10^  psi) 

(10®  psi) 

4 

(10®  psi) 

®3 

(10®  psi) 

Failure 

Mode 

1-26-76 

2325#! 2 

22 

40 

62.68 

1,000 

7.29 

6.95 

n=103 

6.56 

n=10'' 

5.34 

n=105 

Middle 

1-27-76 

2325=5 

22 

40 

69.52 

1,000 

7.90 

7.43 

n=103 

7.02 

0=10"^ 

6.84 

n=117,570 

Middle 

1-29-76 

2335=6 

22 

40 

69.17 

1,000 

7.52 

8.59 

n=103 

7.39 

n-lO'* 

8.20, 

n=105 

End  tab 

2-3-76 

2335=21 

22 

40 

61.36 

1,000 

7,76 

7.33 

n=ic3 

7.35 

n-lO'* 

7.44 

n=105 

End  tab 

2-10-76 

2337=3 

22 

40 

0 

810 

8. 87 

6.32, 

n=10- 

5.67. 

n^-lO'’ 

6.64 

n=10'^ 

End  tab 

2-20-76 

2339=18 

22 

40 

71.22 

1,000 

7.05 

6.34 

n-103 

5.96 

n=10'' 

5.76, 

n=105 

Middle 

2-24-76 

2339=22 

22 

40 

69.98 

1,000 

6.36 

6.88 

n=103 

6.58, 

n=10'' 

6.39 

n=105 

End  tab 

2-25-76 

2340.^8 

22 

40 

78.91 

1,000 

7.05 

6.75 

n=10-> 

6.68, 

n=lQ‘' 

6.64^ 
n  =  10^ 

End  tab 

3-23-76 

2375=7 

22 

40 

72.27 

1,000 

6.71 

6.17 

n--103 

5.52 

li-io'' 

n=105 

Middle 

3-25-76 

2360=20 

22 

40 

70.24 

1,000 

6.04 

6.88 

n=105 

5.95 

n=5xl05 

5.81 

n=10<^ 

End  tab 

3-21-76 

2374#4 

22 

40 

63.96 

1,000 

7.11 

7.57 

n=103 

7.11 

n=10'' 

6.52 

r, =1.1x10- 

End  tab 

I 

3-22-76 

2374fl8 

22 

40 

60.56 

1,000 

6.67 

6.31 

n=103 

6.41 

n“10^ 

6.49^ 

n=105 

End  tab 

3-20-76 

2379^4 

oo 

40 

64.21 

'  ,000 

6.83 

6.78 

n=103 

6 .35 
n=10‘' 

7.52 

11=10-^ 

Middl- 

3-19-76 

2379=14 

22 

40 

74.05 

1 ,000 

6.33 

5.82 

n=103 

5.65 

n=10'' 

6.21 

n=108 

End  tab 
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TABLE  16.  CONSTANT  AMPLITUDE  FATIGUE  TEST 


Type  of  Specimen  -  q/e  [0, +^45,035 
Average  Ultimate  Strength  =  100.4  K$I 
Stress  Ratio  =  0.1  Cyclic  Frequency  =  10  H/ 


Date 

.  (DdvsJ_  _ 

Specimen 

Number 

Mean 

Stress 

(KSI) 

Maximum 

Stress 

(KSI) 

Residual 

Stress 

(KSI) 

Applied 

Eo 

{Ini tial ) 
(106  psi) 

(10®  psi) 

®2 

(10®  psi) 

Failur'j 
(10®  psi  Mode 

9-3-75 

838 

G5I26 

46.93 

85.33 

0 

0.23 

End  tab 

9-9-75 

844 

G5I129 

42.76 

77.74 

0 

18.79 

11.26 

End  tab 

11-10-75 

905 

GSIIj 

43.12 

78.40 

0 

1.21 

7.81 

End  tab 

11-6-75 

901 

651122 

43.12 

78.40 

108.71 

1099.00 

11.47 

10.73 

n’25,000 

9.62 

n^l,099.000 

End  tab 

12-4-75 

929 

65131 

43.12 

78.40 

122.92 

1020.00 

14.90 

14.25 

n-100.000 

12.987 

n=850,000  n= 

13,03 

=1,020,000 

11-24-75 

919 

G5I30 

43.12 

78.40 

144.15 

i071.60 

11  ,63 

10.31 

n=231,500 

9.85 

n=l, 071,600 

10-27-75 

^'^92 

65181 

43.12 

78.40 

127.00 

2283.24 

15.31 

13,51 

nv5,000 

13.54 

n^2.283,?40 

End  tab 

10-23-75 

878 

G5I85 

43.12 

78.40 

123.01 

11-24-75 

879 

65159 

43.12 

78.40 

108.81 

11^1.90 

9.36 

10.37 

n=5,000 

11.91 

n= 789, 000  n 

10.38 

-1.141 ,900 

10-31-75 

8S6 

65 1 70 

43.12 

78.40 

102. 9e 

1065.83 

10.59 

10.59 

n==20,000 

10.37 

1,065.330 

10-24-75 

889 

G5I19 

43.12 

78.40 

0 

.14 

10.51 

End  tab 

9-10-75 

845 

65113 

40.03 

72.33 

81.76 

1.57 

9-3-75 

838 

G5I77 

38.65 

70.03 

97.46 

1926.60 
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ANALYSIS  OF  FATIGUE  DATA 


Jann  N.  Yang 

George  Washington  University 
(Formerly  Virginia  Polytechnic  Institute 
and  State  University) 
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OBJECTIVES 


Strength  Degradation  Model  for  Unnotched 
Composite  Laminates 

(a)  Theoretical  Investigation 

(b)  Experimental  Verification 

Theory  of  Periodic  Proof  Tests  for  Composites 

(a)  Theoretical  Investigation 

(b)  Experimental  Verification 

Reliability  Prediction  and  Cost  Optimization 
for  Composites  Under  Periodic  Proof  Tests  in 
Service 


PRINCIPAL  FINDINGS  AND  CONCLUSIONS 


The  theoretical  residual  strength  degradation 

MODEL  CORRELATES  REASONABLY  WELL  WITH  THE 
EXPERIMENTAL  DATA  FOR  G/E  UNNOTCHED  LAMINATES 

Correlation  between  the  theory  of  periodic 

PROOF  TESTS  FOR  COMPOSITES  AND  THE  EXPERI¬ 
MENTAL  DATA  IS  OUTSTANDING 

Damage  to  G/E  unnotched  laminates  due  to 

REPETITIVE  PROOF  TESTS  APPEARS  TO  BE  NEGLIGIBLE 
FOR  THE  PROOF  LOAD  LEVEL  UP  TO  90%  OF  THE 
MEAN  ULTIMATE  STRENGTH 

Methodologies  established  for  the  reliability 

PREDICTION  OF  COMPOSITES  UNDER  PERIODIC  PROOF 
TESTS  IN  SERVICE  CAN  BE  USED  WITH  CONFIDENCE 


NUMBER  OF  CYCLES 
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RESIDUAL  STRENGTH  AFTER  26,000 
FATIGUE  CYCLE,  KSI 


RESIDUAL  STRENGTH  AFTER  26,000 
FATIGUE  CYCLES  FOR  SURVIVING 
SPECIMENS,  KSI 
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NO.  OF  STRESS  CYCLE  rf 
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DISTRIBUTION  FUNCTION  *"  R  (2T-) 
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RESIDUAL 


NO.  OF  CYCLES 
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RESIDUAL  STRENGTH  IN  H  RESIDUAL  STRENGTH  IN  fl 


DYNAMICS  OF  COMPOSITE  MATERIALS 


by 

George  Herrmann 
Stanford  University 

Figure  Captions 


Fig.  1 

Fig.  2 
Fig.  3 

Fig.  4 
Fig.  5 

Fig.  6 
Fig.  7 
Fig.  8 
Fig.  9 
Fig.  10 
Fig.  11 

Fig.  12 

Fig.  13 

Fig.  14 


Spectral  Surface  for  Harmonic  Waves  in  a  Homogeneous  Body  (i.e. 
frequency  w  as  a  function  of  wave  numbers  ,  K^) 

Laminated  Composite 

Qualitative  Representation  of  Spectral  Surface  for  Laminated  Com¬ 
posite  for  SH-Wave.  Frequency  Q  as  a  Function  of  Wave  Number 
n  (in  periodic  direction)  and  ^  (singular  direction).  B. 
are  Ends  of  Brillouin  zones. 

Spectral  Lines  for  Constant  c  of  Laminated  Composite 

Spectral  Lines  and  Associated  Mode  Shapes  along  Ends  of  Brillouin 
Zones  for  Laminated  Composite 

Curves  of  Constant  Frequency  (Fermi  Lines) 

Spectral  Lines  of  a  Good  Conductor 

Spectral  Lines  of  a  Good  Insulator 

Qualitative  Dependence  of  Phase  Velocity  on  Wave  Number 

Qualitative  Dependence  of  Group  Velocity  on  Wave  Number 

Spectral  Lines  of  Six-Parameter  MES  Theory  as  an  Extension  of 
Effective  Stiffness  Theory 

Spectral  Lines  of  Six-Parameter  Approximate  Theory  as  an  Extension 
of  Kunin  Theory 

Comparison  of  Mode  Shapes  at  End  of  First  Brillouin  Zone.  Layer 
Thicknesses  are  normalized.  Interfaces  are  at  ~1  Primes  refer 
to  Stiffer  Layers 

Equations  of  Motion  and  Frequency  Equation  of  Two  Approximate 
Theories 
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=  Const,  (projection) 


Cg  =  Vk  CO 


1  Spectral  Surface  for  Harmonic  V.'aves  In  a  Homogeneous  Body 
(I.e.  frequency  oi  as  a  function  of  wave  numbers  K  ,  K 


2  Lam I na  ted  Compos  I te 


SYM-SYM’ 

ASYM-ASYM‘ 


ASYM-SYM 

SYM-ASYM 


Fig.  6  CurvRs  of  Con*itant  Frequency  (Fermi  Lines) 
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PERFECT  INSULATOR 


.  10  Qualitative  Dependeilce  of  Group  Velocity  on  Wave  Number 


- EXACT  SOLUTION 

- MATCHED  EFFECTIVE  STIFFNESS  THEORY 

- 'OLD'  EFFECTIVE  STIFFNESS  THEORY 


Fig.  n  Spectral  Lines  of  Six-Parameter  MES  Theory  as  an  Extension  of 
Effective  Stiffness  Theory 


EXACT 

APPROX. 


EXACT 


—  APPROX. 


FJg.  13  Comparison  of  Mode  Shapes  at  End  of  First  Brillouin  Zone,  Layer 
Thicknesses  are  normalized.  Interfaces  are  at  ±1  »  Primes 
Refer  to  Stiffer  Layers 
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Fig.  Equations  of  Motion  and  Frequency  Equation  of  Two  Approximate 
Theories 


CONTINUUM  THEORY  OF  FRACTURE 

Morton  E.  Curtin 
Carnegie-Mellon  University 
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LINEAR  AND  NONLINEAR  EFFECTS 


IN  THE  VIBRATION  OF 
ELASTIC  STRUCTURES 


Grant  AFOSR-73~2479 


Lawrence  W.  Rehfield 
Professor 

School  of  Aerospace  Engineering 
Georgia  Institute  of  Technology 
Atlanta,  Georgia  30332 


NONLINEAR  EFFECTS  IN 

STRUCTURAL  VIBRATIONS 

1)  Drift  of  resonant  frequencies  with 
increasing  amplitude  of  the  motion. 

2)  Alteration  of  time  history  of  transient 
response. 

3)  Jump  phenomena. 

4)  Dynamic  buckling  (oil  canning)  in  the 
presence  of  prestress  or  geometric 
curvature. 

5)  Subharmonic  or  superharmonic  resonances. 

6)  Bounded  limit  cycle  response. 
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NONLINEAR  MATERIAL  EFFECTS  IN 
RESIN  MATRIX  ADVANCED  COMPOSITES 

1)  Different  stiffness  in  tension  and  compression. 

2)  Nonlinear  matrix-dominated  shear  behavior. 

3)  Accentuating  influence  of  moisture  absorption 
and  elevated  temperature  on  the  above. 


OBJECTIVES 


1)  To  determine  the  nature  and  extent  of  the 
influence  of  nonlinear  material  characteristics 
on  vibration  behavior  of  resin  matrix  composite 
structures. 

2)  Assess  the  impact  of  moisture  absorption  and 
elevated  temperature  on  the  above. 
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SCOPE 


1)  Restrict  attention  to  flexural  behavior 
of  simple  beams. 

2)  Consider  [O]  and  lay-ups. 

3)  Perform  trend  analysis. 

4)  Define  and  perform  definitive  experiments. 


REQUIRED  DEVELOPMENTS 

1)  Nonlinear  theory  and  method  of  analysis. 

2)  Uniaxial  nonlinear  material  model. 

a  =  E(8  +  +  pc.^) 

3)  Compressive  test  method. 

4)  Test  method  for  producing  large  amplitude 
structural  vibrations. 
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ANALYTICAL  RESULTS  FOR  FREE 
VIBRATION  OF  SIMPLY  SUPPORTED  BEAMS 


r  -  mode  number 


If  >  0,  structural  response  is  *‘hard." 
If  B  <0,  structural  response  is  *‘soft." 


if  [i  + 


+  (3P|J.  -  2o(^) 


ENVIRONMENTAL  EFFECTS  DUE  TO 

MOISTURE  ABSORPTION  AND  ELEVATED 

TEMPERATURE  ON  FIBER  DOMINATED 
BEAM  VIBRATIONS 


Consider  [0]  laminates  of  T300  Narmco  5208  graphite/ 
epoxy  with  the  following  room  temperature  dry  properties: 

Ej,  =  21.3  X  10^  psi 

E  =  18.8  X  10^  psi 
c 

e  ==  0.007 
a 

Under  the  least  favorable  environmental  conditions, 

E  will  be  reduced  by  about  25  percent.  With  this  assumption 
c 


(U)  )d 

or 


0.94 
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C  O  N  C.L  U  S  I  O  N  S  F  R  O  M  TREND  ANALYSIS 


1)  Nonlinear  material  effects  are  increasingly  important  as 

a.  the  slenderness  ratio  decreases 

b.  the  mode  number  increases 


2)  Graphite/epoxy  |^oJ  beams  exhibit  hardening  behavior. 

3)  Boron/epoxy  |^oJ  beams  exhibit  a  tendency  to  softening  behavior, 

4)  r+^^l  laminates  exhibit  a  pronounced  tendency  to  softening 
behavior. 

5)  Effect  of  moisture/elevated  temperature  appears  to  pose  no 
major  problems  for  fiber  dominated  flexural  vibrations.  Matrix 
dominated  vibrations  cannot  be  rationally  evaluated  at  this  time. 


FUTURE  THEOR  Y/A  NALYSIS  WORK 


1)  Material  behavioral  models  for 

a.  two  and  three  dimensional  characterization  of 
lamina  and  laminates. 

b.  realistic  description  of  degradation  due  to  moisture 
and  elevated  temperature. 

2)  Flexural  vibration  of  laminated  beams  accounting  for 
transferse  shear  and  other  secondary  nonlinear  effects, 
including  boundary  restraint  variation. 

3)  Flexural  vibration  of  angle  ply  panels. 
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20  22  24  26  28  30  32  34 


EXCITATION  FREQUENCY,  Hz. 

STIFFENING  OVERHANG  BEHAVIOR  FOR  MEAN  FORCE, 
Pg»30  POUNDS,  AND  ALTERNATING  FORCE,  P,  14  POUNDS. 

DATA  FOR  INCREASING  FREQUENCY  -  • 

DATA  FOR  DECREASING  FREQUENCY  -Q 


FUTURE  EXPERIMENTAL  WORK 

1)  Develop  experimental  means  of  quantifying 
nonlinear  effects. 

2)  Environmental  simulation  of  moisture/elevated 
temperature  effects 

a.  material  characterization  studies. 

b.  large  amplitude  vibration  experiments. 

c.  damping  experiments. 


3)  Vibration  of  angle  ply  panels. 
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MECHANICS  OF  COMPOSITE  MATERIALS 


G.  A.  Hegemier 

University  of  California,  San  Diego 


RESEARCH  AREAS 


1,  WAVE  PROPAGATION  IN  COMPOSITE  MATERIALS: 

PREDICTION  OF  STRESS,  DEFORMATION. 
THERMAL  FIELDS 

2.  MICROMECHANICS  OF  FAILURE: 

PREDICTION  OF  FAILURE  MODES.  DAMAGE, 
DEGRADED  PROPERTIES 
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GEOMETRIES  CONSIDERED 


LAMINATES: 

PROPAGATION  AT  OBLIQUE  INCIDENCE 


FIBROUS  COMPOSITES: 


ID  ARRAYS 
2D  ARRAYS 
COMPLEX  3D  ARRAYS 


PROPAGATION  AT 
0®  AND  90® 
INCIDENCE 


CONSTITUENT  BEHAVIOR  CONSIDERED 

LINEAR  ELASTIC 

LINEAR  VISCOELASTIC 

ELASTIC-PERFECTLY  PLASTIC 

(FINITE  DEFORMATIONS,  THERMODYNAMICS) 


VALIDITY  TESTS  OF  THEORIES 


1.  COMPARISON  OF  PHASE  VELOCITY  SPECTRA 

2,  COMPARISON  OF  TRANSIENT  PULSE  RESULTS 

COMPARISON  DATA 

1.  exact  SPECTRA  FOR  ELASTIC  LAMINATES  AT 
OBLIQUE  INCIDENCE  AND  ELASTIC  FIBERS  AT  90® 

2.  EXACT  TRANSIENT  PULSE  RESULTS  FOR  ELASTIC 
LAMINATES  AT  0®  INCIDENCE 

3.  DETAILED  NUMERICAL  TRANSIENT  PULSE  RESULTS 
FROM  2D  NONLINEAR  FINITE  DIFFERENCE  CODES  FOR 
FINITE  AMPLITUDE  90®  PROPAGATION  IN  LAMINATES 
AND  UNIDIRECTIONAL  FIBER  ARRAYS 

.  EXPERIMENTAL  RESULTS  ON  LAMINATES  AND  FIBROUS 
COMPOSITES  (PHASE  VELOCITY  AND  TRANSIENT  PULSE) 


2  83 


4. 


LINEAR  AND  NONLINEAR  WAVE  PROPAGATION  IN  LAMINATES 


Fig.  1.  Geometry  and  Coordinate  System  for 

Wave  Guide- type  Propagation  in  Laminates* 


TYPICAL  MIXTURE  FORM  FOR  LAMINATES 


THORNEL  -  CARBON  PHENOLIC  90*  LAMINATE 

2(h"'+ , 

^CXACT 

X 

.MIXTURE  THEORY 

f  3"*  ORDER 

THEORY  \ 

experiment  ^ 

(WHITTIER  A  PECK) 

7^ 

1**  ORDER  THEORY*^ 

/ 

MODIFIED  1**  ORDER 

THEORY  ^ 

/  \ 

12  3  4 

fREOUEKCt  (  10*  J  CPS 


I 


BORON- CARBON  PHCNOLIC  90*  LAMINATE 
^"^^^MIXTURE  THECRy  X 

ORDER 

THEORY 

EXPERIMENT 
(WHITTIER  A  PECK)  j 

I*’ ORDER  THEORY  \ 


MOOlPlEO  ("ORDER  THEORY 


2  9  4  9 

FREOUENCY  (  tO^  1  CPS 


EXPERIMENT  /  f  N. 

(WHITTIER  A  PECK)  /  \ 


Fig.  2a,  b.  Lowest  symmetric  mode  of  the  phase  velocity  spectrum  for  typical  layered  composites. 


\  "x 


EFFECTIVE  STIFFNESS 
4^*^  THEORY 


\ 

MODIFIED 
1*^  ORDER  THEORY 


t**ORDER  THEORY^, 


2  ORDER  THEORY 


> 

ui  0 
tn 


1234567 

NONOIMENSION/.L  WAVE  NUMBER  ,  K*  *  2Kh*** 


MODIFIED  r'  ORDER  THEORY  - 

.EFFECTIVE  STIFFNESS  THEORY 


,2 '“ORDER  THEORY 


r  ORDER  THEORY, 
rvArT  / 


01234  3676 

NONOIMENSIONAL  WAVE  NUMBER  ,  K*  •  ZKN**' 

Fig.  5a,  b.  Phase  velocity  ^ectrum  comparisons  for  several  theories. 


0 


Fig.  9. 


1.0  2.0  3.0  4.0  5.0  6.0 

Time  (usee) 


Comparison  of  theory  with  experiment  No.  7, 


289 


Phase  Velocity  vs*  Wave  Number  for  Several  Propagation  Angle* 


LINEAR  AND  NONLINEAR  WAVE  PROPAGATION  IN 
UNIDIRECTIONAL  FIBROUS  COMPOSITES 


Fig.  12.  Unidirectional  fiber- reinforced  composite  with 
hexagonal  aeray. 
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Fig.  13.  Comparison  of  mixture  and  experimental  rear  surface 
velocity  results  from  quartz  fiber  reinforced  phenolic. 


Phenolic 


Ar  =  0.01  cm 
Ax  =  0.02  cm 


Problem  #1 

Vk  =  0.05  _ 

cm//i.sec 

. -j _ i — 

i.O  2.0 
time  (/xsec) 


Problem  #2 

Vx  =  0.05 
cm//tsec 

t.O 

time  (/isec) 


Fig.  H. 


Geometry  for  2D  CRAM  calculation. 
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Velocity  in  Nicnulic  SliealU  (tUtcrial  2)  at 


ComparLaon  of  TINC  and  average  CRAM  particle  Comparison  of  TINC  and  average  CRAM  particle 

velocity  histories  for  radial  fiber*  velocity  histories  for  sheath* 


Radius  (cm) 


Comparison  of  TINC  and  CRAM  predictions  of 
axial  stress  (x  =  0*  2  cm,  t  =  1.0  psec). 


LINEAR  AND  NONLINEAR  WAVE- PROPAGATION  IN  3D  COMPOSITES 


Fig.  18.  AVCO  3DQP  "Black"  material. 


ed  Rear  Surface  Velocity,  V/Y, 


Input^StrcsSp 


Fig*  21#  Input  streBS  profiles  for  AFWL  shot  326, 


Fig,  22,  Stress  history  at  rear  of  0.  254  in,  3DQP 
plus  0,060  in.  plexiglass  (shot  326), 
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Avertge  Stress,  kb 


Fig#  23.  Stress  history  at  rear  of  0.  154  in.  3DQP 
plus  0.060  in.  plexiglass  (shot  314). 
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1/1 


Fig,  24,  Temperature  distribution  at  T  =  1,  1.5  and  2  for  similar 
material  properties  and  small  duration  of  input  pulse 
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T/X 


0  .2  4  6  8 

y/A 

Fig.  25,  Temperature  distribution  at  T  =  i  for  widely 
differing  material  properties  and  small  t^. 
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RESPONSE  OF  CANTILEVER  PLATE  (FAN  BLADE) 
DUE  TO  IMPACT  LOADING 

P.  C.  CHOU 
DREXEL  UNIVERSITY 


Objectives 

I .  Design  Curve 

To  obtain  a  design  curve  for  estimation  of  peak  stress  in 
structural  response  (in  milliseconds). 

2*  Finite-difference  Code 

To  develop  a  finite-difference  computer  code  for 
anisotropic  material,  for  local  response 
(in  micro-seconds). 


1.  Beam  -  Review  of  general  approach. 
Analytical  Approach  -  5  methods 

*  Lumped  spring  energy  conserved 

*  lumped  mass  momentum  conserved 

*  Two-degree-of -freedom 
•McQuillen  solution 
•Timoshenko  solution 

Experimental  Approach 

'  Maximum  strain  measurement 
'  Failure  experiment 

2.  Simply  Supported  Plate 

3.  Cantilever  Plate  (Fan  blade) 
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Kj  =  beam  spring  constant 
2.  Dynamic  deflection  curve  assumed  same  as  static. 
K  .  iO. 

^x=c  "  (^•'d)'^lmax 


3.  Define  generalized  strain: 
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Imoact  Design  Curve  for  Cantilever  Plate 


lo 


I.  Conservation  of  energy: 


I  >"2 ''^-2^1 


2.  Perform  experiments  to  determine  spring  constant  and  relationship 
between  strain  and  deflection  by  applying  static  load  at  impact  point 

P  =  Ki  Wi 


P  =  K2  € 


•*.  €  =  J- 


K2  I 


3.  Define  generalized  strain: 


€  =  £ 
where 


yD|i//>h 


Final  result: 


where  C^  = 


D||  ab  K| 
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Finite-Difference  Code  (ANEL) 


(For  Local  Response) 


Approach  :  Based  on  HEMP 

HEMP  -  Lagrangian,  elastic-plastic-hydrodynamic 
isotropic,  plane  strain  or  cylindrical 
symmetry. 


ANEL  -  Lagrangian,  elastic,  monoclinic  (plane 
strain  or  plane  stress),  or  transversely 
isotropic  (axisymmetric  problem) 


Results 


1.  Comparison  with  other  calculations. 

2.  Comparison  with  experiments. 


ANEL 

A  Finite-Difference  Lagrangian  Code 
Governing  Equations 
Cauchy's  Equation  of  Motion 


D^x 


r  ^  yij.j  X2  X 


Conservation  of  Mass 


Conservation  of  Energy 


^33 


Strain  Displacement 


•33 


1  ^  •  •  V 


^2 


Constitutive  Relations 

^ij  "I  ‘  "^ij] 

^  3  ^ijkk  ^qqkk  ^ 

'  3  f  ^iikjt,  ^kj  *  ^iikk  ^mj 
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Loading  Function  Used  in  ANEL  Calculation 
(To  be  Compared  with  the  Static  Solution  by  Greszczuk) 


Q Present  Calculation, 
^  ^  Dynamic 

X - Static,  Finite  Element 

(Greszczuk) 


Experiment 


Comparison  for  Edge  Impact  of  Composite  Plate 


Conclusions 


I.  Design  Curve 

♦  A  method  is  deveioped  for  constructing  design  curves 
for  piate  under  impact.  For  a  given  type  of  structure, 
calculation,  or  testing,  of  a  few  cases  is  sufficient; 
it  can  be  used  to  estimate  other  sizes  and  other  impact 
conditions. 


•  Maximum  strain  proportional  to  velocity  and  to  square 
root  of  impactor  mass. 

Olrax  OC  V 


2.  Anisotropic  Code 

♦  An  anisotropic  finite-difference  computer  code  is 
developed  for  local  impact  response  calculation. 


•  dynamic  stress  higher  than  quasi-static  stress 


•  May  be  extended  to  include  plastic  and 
hydrodynamic  constitutive  equations. 
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MECHANICS  OF  COMPOSITE  MATERIALS 
WITH  DIFFERENT  MODULI  IN 
TENSION  AND  COMPRESSION 

ROBERT  M.  JONES 
SOUTHERN  METHODIST  UNIVERSITY 
DALLAS,  TEXAS 


OUTLINE 


•  INTRODUCTION 

•  NONLINEAR  MULTIMODULUS  MATERIAL  MODEL 

•  MATERIAL  MODELING  ACCOMPLISHMENTS 

•  MATERIAL  MODELING  IN  PROGRESS 

•  SUMMARY 
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INTRODUCTION 


BILINEAR  MODEL  FOR  MULTIMODULUS  MATERIALS 
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TABLE  1 


TENSION  AND  COMPRESSION  MODULI  RELATIONSHIPS 
FOR  SEVERAL  COMMON  COMPOSITE  MATERIALS 


MATERIAL 

FIBROUS 

OR 

GRANULAR 

REPRESENTATIVE 

MODULI 

RELATIONSHIP 

GLASS/EPOXY 

FIBROUS 

BORON/EPOXY 

FIBROUS 

GRAPHITE/EPOXY 

FIBROUS 

CARBON/CARBON 

FIBROUS 

ZTA  GRAPHITE 

GRANULAR 

ATJ-S  GRAPHITE 

GRANULAR 

EFFECT  OF  MOISTURE  AND  TEMPERATURE 
ON  STRESS  “STRAIN  CURVE  NONLINEARITY 


r 
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TYPICAL  STRESS -STRAIN  BEHAVIOR 
OF  FIBER -REINFORCED  COMPOSITE  AAATERIALS 


OBJECTIVES 

MECHANICS  OF  COMPOSITE  MATERIALS  WITH  DIFFERENT  MODULI  IN  TENSION  AND  COMPRESSION 

•  DEVELOP  MATERIAL  MODEL  FOR  NONLINEAR  MULTIMODULUS  BEHAVIOR 
OF  VARIOUS  CLASSES  OF  COMPOSITE  MATERIALS 

•  ANALYZE  STRESS  EQUILIBRIUM  BEHAVIOR  OF  SOLID  BODIES 

(E.G.,  THERMAL  STRESSES  SIMULTANEOUS  TENSILE  AND  COMPRESSIVE  STRESSES) 

0  LINEAR  ELASTIC  MULTIMODULUS 
0  NONLINEAR  ELASTIC 
0  NONLINEAR  ELASTIC  MULTIMODULUS 

0  ANALYZE  BENDING,  BUCKLING,  AND  VIBRATION  BEHAVIOR  OF  LAMINATED  PLATES  AND  SHELLS 
0  LINEAR  ELASTIC  (LAMINATION  ASYMMETRIES) 

0  LINEAR  ELASTIC  MULTIMODULUS 
0  NONLINEAR  ELASTIC 
0  NONLINEAR  ELASTIC  MULTIMODULUS 
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JONES  -  NELSON 
NONLINEAR  MATERIAL  MODEL 


ORTHOTROPIC  STRESS-STRAIN  RELATIONS 

AXISYMMETRIC  BEHAVIOR 
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MATERIAL  PROPERTY  -  ENERGY  RELATIONS 


MATERIAL  PROPERTY:  =  A:  1  -  B 


(t) 


=  INITIAL  ELASTIC  VALUE 
=  INITIAL  CURVATURE  OF  Cr-6  CURVE 


RATE  OF  CHANGE  OF  CURVATURE 


=  STRAIN  ENERGY 
=  UNIT  STRAIN  ENERGY 


STRAIN  ENERGY  FUNCTION 


U  “  ^9  Xrz) 
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UNIAXIAL  INTERPRETATION  OF  THE  STRAIN  ENERGY 


REPRESENTATION  OF  STRESS-STRAIN  BEHAVIOR 

DIRECT  MODULI 
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ITERATION 

PROCEDURE 

PMD  =  PRINCIPAL 
MATERIAL 
DIRECTIONS 

PSD  =  PRINCIPAL 
STRESS 
DIRECTIONS 


CONVERGENCE  OF  ITERATION  PROCEDURE 


STRAIN 


322 


ACTUAL  VS.  EXTENDED  APPROXIMATE 
MECHANICAL  PROPERTY  VS.  U  CURVE 


MECHANICAL 

PROPERTY 


STRAIN  ENERGY,  U 


LIMIT  OF  MEASURED  DATA  FOR 
MECHANICAL  PROPERTY  VS.  U  CURVE 
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LIMITS  OF  STRESS  -  STRAIN  DATA  FOR 
DIFFERENT  PRINCIPAL  MATERIAL  DIRECTIONS 


EXTENDED  STRESS -STRAIN  CURVE  APPROACH 

JONES-NELSON -MORGAN  A,  B,  AND  C  PLUS  INPUT  SLOPE 


€ 
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EXTRAPOLATION  OF  MECHANICAL  PROPERTY  BEHAVIOR 

JONES- NELSON -MORGAN  A,  B,  AND  C  PLUS  INPUT  SLOPE 


DEFINITION  OF  A  MATERIAL  MODEL 


•  OFTEN  - 

THE  STRESS -STRAIN  RELATIONS  FOR  THE  MATERIAL 


•  HERE  - 

THE  SET  OF  MECHANICAL  PROPERTY  VERSUS  ENERGY  EQUATIONS 
PLUS  THE  ASSOCIATED  ITERATION  PROCEDURE  FOR  THE  MATERIAL 
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STEPS  TO  DEFINE  AND  VALIDATE  MATERIAL  MODELS 

(1)  FIT  MECHANICAL  PROPERTY  VERSUS  ENERGY  EQUATIONS 

TO  UNIAXIAL  STRESS  -  STRAIN  DATA  IN  PMD* 

=1^  MATERIAL  MODEL  IS  DEFINED 

(2)  USE  MATERIAL  MODEL  TO  PREDICTED  STRAINS 
UNDER  UNIAXIAL  LOADING  IN  OTHER  THAN  PMD* 

AND  COMPARE  PREDICTED  WITH  MEASURED  STRAINS 

(3)  USE  MATERIAL  MODEL  TO  PREDICT  STRAINS 

UNDER  BIAXIAL  LOADING  IN  PMD*  AND  OTHER  THAN  PMD* 
AND  COMPARE  PREDICTED  WITH  MEASURED  STRAINS 

(4)  IF  ALL  COMPARISONS  ARE  SATISFACTORY, 

THEN  MATERIAL  MODEL  IS  REGARDED  AS  VALIDATED 

*PMD  =  PRINCIPAL  MATERIAL  DIRECTIONS 
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MATERIAL  MODELING  ACCOMPLISHMENTS 


MATERIAL 

TYPE 

OF 

MATERIAL 

ATJ-S  GRAPHITE 

TRANSVERSELY 

ISOTROPIC 

BORON/EPOXY 

ORTHOTROPIC 

LAMINA 

GRAPHITE/EPOXY 

ORTHOTROPIC 

LAMINA 

BORON/ALUMINUM 

ORTHOTROPIC 

LAMINA 

CARBON -CARBON 

ORTHOTROPIC 

LAMINA 

AND 

ANISOTROPIC 
3-D  WEAVE 

MATERIAL  MODEL 


MULTIMODULUS  NONLINEAR  NONLINEAR 
TENSION  I  COMPRESSION  I  ELASTIC  MULTIMODULUS 


NUMBER  OF 
NONLINEARITIES 


MATERIAL  MODELING  ACCOMPLISHMENTS 


COMPARISON  WITH  EXPERIMENT 


MATERIAL 


UNIAXIAL 


BIAXIAL 


PMD*  OFF-AXIS  TUBE  DISK  LAMINATE 


ATJ-S  GRAPHITE  / 
BORON/ EPOXY  / 
GRAPHITE/EPOXY  / 
BORON/ALUMINUM  / 
CARBON-CARBON 


*PHD  =  PRINCIPAL  MATERIAL  DIRECTIONS 


OFF-AXIS  Cr-€  BEHAVIOR 
TENSION,  a=70° 


OFF-AXIS  Cr-6  BEHAVIOR 
COMPRESSION,  a=70°  WCM 


HOLLOW  GRAPHITE  BIAXIAL  TEST  SPECIMEN 


WEDGE-SHAPED  DISK  AND  IDEALIZATION 
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DISK  INMER  DIAMETER  CHANGE  VS.  TIME 


INNER 

DIAMETER 

CHANGE 

inches 


TIME,  seconds 


OFF-AXIS 
STRESS -STRAIN 
BEHAVIOR 


2 


4 


10 


MATERIAL  MODELING  IN  PROGRESS 


MATERIAL  MODELING  IN  PROGRESS 

•  NONLINEAR  LAMINATE  DEFORMATION  AND  BUCKLING 

•  BORON/EPOXY 

•  GRAPHITE/EPOXY 

•  BORON/ALUMINUM 

•  EXTENSION  OF  MODEL  TO  CARBON-CARBON 

•  BEHAVIORAL  CHARACTERISTICS  OF  CARBON-CARBON 

•  BASIC  CHARACTERISTICS  OF  MODEL 

•  APPARENT  FLEXURAL  MODULUS  AND  STRENGTH 
OF  MULTIMODULUS  MATERIALS 

•  COMPARISON  OF  PREDICTED  AND  MEASURED  OFF-AXIS  STRAINS 
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SEQUENCE  OF  POINTS  IN  AN  ITERATION  PROCEDURE 
TO  FIND  STRAINS  CORRESPONDING  TO  A  GIVEN  LOAD 


PREDICTED  AND  MEASURED 
LAMINATE  LOAD-STRAIN  BEHAVIOR 


r 
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BUCKLING  LOAD  SEARCH  PROCEDURE 


BUCKLING  OF  UNIDIRECTIONAL  0°  LAMINATES 


t,  INCHES 
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MODELING  OF  CARBON-CARBON 


•  NONLINEAR  ORTHOTROPIC  OR  O' 

ANISOTROPIC  MATERIAL 

•  MODELED  AS  A  NONLINEAR 

STRESS-STRAIN  CURVE  WITH 
A  LINEAR  PORTION  ABOVE 
A  TRANSITION  STRESS  £ 


MODELING  OF  CARBON-CARBON,  continued 

•  JNMDATA  COMPUTER  PROGRAM  (AUTOMATED  MODELING  AND  INPUT) 

•  NOW  FOR  TRANSVERSELY  ISOTROPIC  MATERIALS 

•  EXTENDING  TO  ORTHOTROPIC  MATERIALS 

•  APPARENT  FLEXURAL  MODULUS  AND  STRENGTH 

•  NOT  INDEPENDENT  PROPERTIES  OF  MULTIMOCULUS  MATERIALS 

•  CAN  BE  PREDICTED  WITH  AVAILABLE  DATA 

•  SHOULD  NOT  BE  MEASURED 

•  COMPARISON  OF  PREDICTED  AND  MEASURED  OFF-AXIS  STRAINS 

•  USE  SoRI  DATA 

•  22  1/2°  and  45°  OFF-AXIS  DATA 
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PRINCIPAL  FINDINGS  AND  CONCLUSIONS 

MECHANICS  OF  COMPOSITE  MATERIALS  WITH  DIFFERENT  MODULI  IN  TENSION  AND  COMPRESSION 
•  MATERIAL  MODEL 

•  STRESS-STRAIN  RELATIONS  FOR  MATERIALS  WITH  DIFFERENT  MODULI  IN  TENSION  AND 
COMPRESSION,  AIAA  JOURNAL.  FALL  1976 

•  A  NEW  MATERIAL  MODEL  FOR  THE  NONLINEAR  BIAXIAL  BEHAVIOR  OF  ATJ-S  GRAPHITE 
JOURNAL  OF  COMPOSITE  MATERIALS.  JANUARY  1975. 

•  FURTHER  CHARACTERISTICS  OF  A  NONLINEAR  MATERIAL  MODEL  FOR  ATJ-S  GRAPHITE. 
JOURNAL  OF  COMPOSITE  MATERIALS.  JULY  1975. 

•  MATERIAL  MODELS  FOR  NONLINEAR  DEFORMATION  OF  GRAPHITE.  AIAA  JOURNAL.  JUNE  1976. 

•  THEORETICAL-EXPERIMENTAL  CORRELATION  OF  MATERIAL  MODELS  FOR  NONLINEAR 
DEFORMATION  OF  GRAPHITE,  AIAA  JOURNAL,  SEPTEMBER  1976. 

SOLID  BODIES 

•  ABOVE  REFERENCES  PLUS 

•  NONLINEAR  DEFORMATION  OF  A  THERMALLY  STRESSED  GRAPHITE  ANNULAR  DISK 
(SUBMITTED  TO  AIAA  JOURNAL) 

•  APPARENT  FLEXURAL  MODULUS  AND  STRENGTH  OF  MULTIMODULUS  MATERIALS 
(SUBMITTED  TO  JOURNAL  OF  COMPOSITE  MATERIALS) 


O  LAMINATED  PLATES  AND  SHELLS 

•  LAMINATION  ASYMMETRIES 

•  BUCKLING  AND  VIBRATION  OF  UNSYMMETRICALLY  LAMINATED  CROSS-PLY  RECTANGULAR  PLATES, 

AIAA  JOURNAL,  DECEMBER  1973. 

•  BUCKLING  AND  VIBRATION  OF  ANTISYMMETRICALLY  LAMINATED  ANGLE-PLY  RECTANGULAR  PLATES, 
JOURNAL  OF  APPLIED  MECHANICS.  DECEMBER  1973. 

•  BUCKLING  AND  VIBRATION  OF  CROSS-PLY  LAMINATED  CIRCULAR  CYLINDRICAL  SHELLS 
AIAA  JOURNAL,  MAY  1975. 

•  DEFLECTION  OF  UNSYMMETRICALLY  LAMINATED  CROSS-PLY  RECTANGULAR  PLATES, 

PROCEEDINGS  OF  12th  ANNUAL  MEETING  OF  THE  SOCIETY  OF  ENGINEERING  SCIENCE,  OCTOBER  1975. 
O  LINEAR  ELASTIC  MULTIMODULUS 

•  BUCKLING  OF  STIFFENED  LAMINATED  COMPOSITE  CIRCULAR  CYLINDRICAL  SHELLS  WITH 
DIFFERENT  MODULI  IN  TENSION  AND  COMPRESSION,  AFOSR-TR-75-0547 ,  FEBRUARY  1975. 

•  BUCKLING  OF  LAMINATED  COMPOSITE  CIRCULAR  CYLINDRICAL  SHELLS  WITH  DIFFERENT  MODULI 
IN  TENSION  AND  COMPRESSION,  PROCEEDINGS  OF  1975  INTERNATIONAL  CONFERENCE  ON 
COMPOSITE  MATERIALS.  APRIL  1975. 

•  BENDING  AND  EXTENSION  OF  CROSS-PLY  LAMINATES  WITH  DIFFERENT  MODULI  IN  TENSION 
AND  COMPRESSION,  PROCEEDINGS  OF  17th  AIAA/ASME  SDM  CONFERENCE,  MAY  1976. 

•  NONLINEAR  ELASTIC 

•  ANALYSIS  OF  NONLINEAR  STRESS-STRAIN  BEHAVIOR  OF  FIBER-REINFORCED  COMPOSITE  MATERIALS, 
PROCEEDINGS  OF  17th  AIAA/ASME  SDM  CONFERENCE,  MAY  1976. 
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SUMMARY 


•  SIGNIFICANT  MULTIMODULUS  EFFECTS  FOR  GRANULAR  COMPOSITE  MATERIALS 

•  MIXED  MULTIMODULUS  EFFECTS  FOR  FIBER-REINFORCED  COMPOSITE  MATERIALS 

•  LOW  FOR  BORON/EPOXY  AND  GRAPHITE/EPOXY  AT  ROOM  TEMPERATURE  AND  LOW  HUMIDITY 

•  (EXPECTED)  HIGHER  FOR  BORON/EPOXY  AND  GRAPHITE/EPOXY  AT  ELEVATED  TEMPERATURE  AND  HUMIDITY 

•  (EXPECTED)  HIGH  FOR  CARBON-CARBON  UNDER  ALL  CONDITIONS 

•  MULTIMODULUS  AND  NONLINEAR  EFFECTS  ARE  SUSCEPTIBLE  TO  RATIONAL  ANALYSIS 

•  FUTURE  WORK  WILL  BE  CONCENTRATED  ON  CARBON-CARBON 
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